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Abstract 



This thesis examines the non-leptonic 5-decays within QCD factorisation 
and beyond, to challenge the assumptions and limitations of the method. 
We analyse the treatment of the distribution amplitudes of light mesons 
and present a new model described by simple physical parameters. The 
leading twist distribution amplitudes of light mesons describe the leading 
non-perturbative hadronic contributions to exclusive QCD reactions at large 
energy transfer, for instance electromagnetic form factors. Importantly, they 
also enter into the two-body B decay amplitudes described by QCD factori- 
sation. They cannot be calculated from first principles and are described 
by models based on a fixed-order conformal expansion, which is not always 
sufficient in phenomenological applications. We derive new models that are 
valid to all orders in the conformal expansion and characterised by a small 
number of parameters related to experimental observables. 

Motivated by the marginal agreement between the QCD factorisation results 
with the experimental data, in particular for B — > tttt, we scrutinise the in- 
calculable non-factorisable corrections to charmless non-leptonic decays. We 
use the available results on B — > to extract information about the size 
and nature of the required non-factorisable corrections that are needed to 
reconcile the predictions and data. We find that the best-fit scenarios do 
not give reasonable agreement to 2er until at least a 40% non-factorisable 
contribution is added. Finally we consider the exclusive B — > V^y decays, 
where we analyse the recently updated experimental data within QCD fac- 
torisation and present constraints on generic supersymmetric models using 
the mass insertion approximation. 
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Introduction 



The quest for comprehension of the physical world around us has dominated 
philosophy and science since the beginning of history. The ancient peoples 
looked up to the stars and the Gods for answers. The philosophers called 
upon the elements of nature - the Ancient Greeks, such as Leucippus, Dem- 
ocritus or Epicurus were the first to analyse and categorise the nature of all 
things. The belief that matter was composed of the four fundamental ele- 
ments: earth, air, fire and water survived for over two thousand years. It was 
they who introduced the first "elementary particles" , the indivisible atomos. 
As mysticism and natural philosophy morphed into scientific principle, the 
foundations of modern particle physics were established; beginning with the 
physical theories of Newtonian mechanics and gravity to the work of Thomp- 
son and Rutherford which allowed physics to leave the 19th century knowing 
the "indivisible atoms" were not in fact, indivisible at all. 

The birth of the 20th century and the pioneering study of quantum mechanics 
was followed in the 1950s and 60s by the discovery of a bewildering variety 
and number of particles, once physicists had discovered that protons and 
neutrons were themselves composite particles. Yet, as with Mendeleev's pe- 
riodic table of 1869, these particles could again be categorised and classified 
revealing a simple underlying structure in terms of a few elementary parti- 
cles. The ideal to capture the beauty and complexity of nature in a simple 
mathematical form culminated in the formulation of the Standard Model of 
particle physics. This is based upon quantum field theory and gauge princi- 
ples - combining the electromagnetic, weak and strong interactions into one 
unified theory. To date the Standard Model is the most successful math- 
ematical theory of particle physics ever created. Particle accelerators have 
probed scales down to 10~ 16 cm, and this theory is consistent with virtually all 
physics down to this scale - indeed the Standard Model has passed indirect 
tests which probe even shorter distances. 

The Standard Model has 25 elementary particles, so perhaps one may ask 
"how deep the rabbit hole goes"? The search for another layer of structure 
beneath the quarks and leptons has however so far been fruitless. More im- 
portantly, the Higgs boson which is required to generate mass for all of the 
Standard Model particles, has still not been discovered in collider experi- 
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ments. There are a number of other unresolved issues with the Standard 
Model, that have led to recent conclusions that it cannot be complete, and 
is more likely to be an "effective theory" of some more encompassing theory 
at higher energy. These issues include some fundamental questions - such as 
the large number of arbitrary parameters that appear in the Standard Model 
Lagrangian. Why are there copies of quarks and leptons in three generations? 
Why are the masses split in the way they are? Why is the weak scale so dif- 
ferent from the expected Planck scale? These questions cannot be answered 
with our current understanding. On top of all these questions, the prob- 
lem of including General Relativity with the Standard Model still remains 
- combining gravity and quantum mechanics produces a non-renormalisable 
quantum field theory - this again suggests the presence of another theory, 
and hence some new physics at some higher energy. 

Many ideas have been put forward to solve these problems, the most popular 
being supersymmetry, string theory or extra dimensions. Yet to date there 
is no theory that is truly simpler or less arbitrary than the Standard Model. 
More importantly, there is little experimental evidence for the existence of 
the holy grail of "new physics" and to guide the directions of the model 
builders. There are perhaps a handful of discrepancies between the theoreti- 
cal predictions and experimental results, none of which have been significant 
enough to claim as new physics. 

Within this thesis I work predominantly within the Standard Model, taking 
some brief forays into generic supersymmetric models. To summarise the 
motivation in a line, I say the following: 'before any claims of new physics 
can be made, we must be absolutely sure that our theory predictions are as 
accurate as we can make them'. The general aim being to better understand 
some elements of the Standard Model - specifically relating to the weak 
decays of 5-mesons. 

Beauty decays are a rich and powerful playground for studying flavour physics 
and CP violation - we can test predictions and constrain Standard Model 
parameters. From experimental observables (decay rates, parameters etc.) 
we can perform indirect searches for new physics via any measured deviation 
from the Standard Model expectation. Due to the rapid decay of the top 
quark, the b quark is the heaviest quark to bind into mesons which can 
be observed. Theoretical techniques based upon an expansion in the heavy 
b quark mass can provide vast simplifications of calculations and model- 
independent predictions of observables. 

The theoretical framework within which any 5-physics analysis is based is 
that of effective field theory, where the operator product expansion and renor- 
malisation group evolution are invaluable techniques which we exploit. These 
encompass the concepts of factorisation which allow the calculation of de- 
cay amplitudes to be separated into perturbatively calculable short-distance 
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Wilson coefficients, and long-distance matrix elements. We need some non- 
perturbative method such as QCD sum rules or Lattice QCD to calculate 
these matrix elements fully. For the exclusive two-body decays of the B- 
meson we can make use of the powerful technique of QCD factorisation which 
makes use of the hierarchy between the heavy b quark mass mj, and the in- 
trinsic scale of QCD, Aqcd- This gives a factorisation formula for the evalua- 
tion of hadronic matrix elements, and is a further separation of long-distance 
contributions from a set of perturbative short- distance contributions. The 
long-distance part must still be calculated via some non-perturbative tech- 
nique, but is of much less complexity than the original matrix element. 

QCD factorisation was developed by Beneke, Buchalla, Neubert and Sachra- 
jda (BBNS) [1-3] for exclusive decays of the B into two mesons. The method 
is however not without its uncertainties, nor gaps where calculations cannot 
be completed without some model-dependent assumptions. It is these un- 
certainties, in the decays of the 5-meson into non leptonic final states, that 
we discuss in this work. 

The main new result of this work involves examining the uncertainty from 
one of the important non-perturbative inputs in the factorisation formulae, 
namely the light-cone distribution amplitudes of the light mesons. We de- 
velop and introduce a set of new models for the distribution amplitude which 
improve upon the truncated conformal expansion that is widely used in the 
literature. We quantitatively take into account contributions from higher or- 
der moments in the conformal expansion by a resummation to all orders. We 
show how these models can be expressed in terms of a few simple physical 
parameters that are directly related to experimental observables. 

We then present an analysis of the second major source of uncertainty to 
QCD factorisation, namely corrections of order (Aq CD /771;,) to the leading 
calculation performed in the heavy quark limit (mj — > oo), that cannot be 
calculated in a model-independent manner. Evidence suggests that QCD fac- 
torisation may considerably underestimate these corrections in some decay 
channels, specifically B — > irn. We construct various scenarios of additional 
non-factorisable contributions and use the wealth of experimental data cur- 
rently available from the B factories and accelerator experiments to quantify 
how good the agreement is with the experimental measurements. 

Finally, we present an analysis using recently released (August 2005) results 
on the B — > (p, uj) 7 decays, based on the extension of QCD factorisation to 
radiative decays B — > V7 by Bosch and Buchalla [4]. We also consider the 
possibility of new physics in the b — > dj and b — * 37 decays, and use the new 
experimental results to constrain contributions from generic minimal super- 
symmetric models using the mass insertion approximation. 

The structure of the thesis can be outlined as follows: 
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We begin with a whistlestop tour of the basic concepts of the Standard Model 
and the theoretical techniques required for an analysis in effective field theory; 
renormalisation group perturbation theory and a discussion of the AB = 1 
effective Hamiltonian. We conclude this introductory chapter with a brief 
introduction to two important non-perturbative techniques: QCD sum rules 
and Lattice QCD. We make use of many results from these methods and so 
it is important to understand their basis and limitations. 

Chapter 2 provides an in-depth discussion on the formulation of QCD fac- 
torisation. We begin by placing it in context with a discussion of naive 
factorisation, and then follow with details of the structure of the factorisa- 
tion formula and its input. We introduce the "non-factorisable corrections" 
and the problems with calculation of power-suppressed diagrams in a model- 
independent manner. We also include the isospin decomposition of B — > nn 
amplitudes which is used extensively in our Chapter 4 analysis, and finally 
conclude with the limitations of the QCD factorisation framework. 

Chapter 3 presents the main result of this work: the development of a new 
resummed model of the light-cone distribution amplitude for light mesons. 
We present a detailed discussion of the treatment of the DA, beginning with 
the application of conformal symmetry techniques, and the expression of the 
DA as a partial wave expansion in conformal spin. In order to place our 
new models in context we give a brief review of the literature and previous 
constructions with specific reference to the pion wavefunction, before going 
on to present full details of our new models. These sections are primarily 
based on the work presented in [5]. We discuss the full implications of the 
new models, and give numerical results to show how variation of the size of 
contribution from the higher-order moments can affect the branching ratios 
and CP asymmetry predictions, using the example of B — > nn. 

We then go on in Chapter 4, to present our analysis on the non-factorisable 
corrections to the charmless 5-decays. We investigate how the extensive 
experimental information available can be used to extract information about 
the size and nature of non-factorisable corrections that are needed to provide 
agreement between the prediction and measurement of B — > nn branching 
fractions and CP asymmetries. We split the B — > tttt isospin amplitudes 
into factorisable and non-factorisable parts - the former being calculated via 
QCD factorisation and the later fitted to the experimental data. We show 
that there is evidence for sizable non-factorisable corrections in the B — > nn 
system, and discuss the application of these results to B — > nK. We also 
discuss the possibility and likelihood of a charming penguin contribution to 
the charmless decays. 

Our final analysis is presented in Chapter 5, and discusses the exclusive 
radiative decays B — > (p, oo)^ and B — > K*^y. We examine the predictions 
from QCD factorisation in the context of the recent experimental results for 
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this system. These decays are rare decays occurring with branching fractions 
of 10~ 5 or less, as they occur only at the loop level within the Standard Model. 
They £1X6 clS such, very sensitive to the possibility of new physics. After 
introducing the basic principles and motivation for one of the most popular 
extensions of the Standard Model - the Minimal Supersymmetric Standard 
Model (MSSM) - we discuss how the mass insertion approximation can be 
utilised to constrain the parameter space for new physics in the b — > (d, s) 
transitions. We give graphical constraints on insertions Sf 3 and 6$ 3 . 

Finally, we present our conclusions and outlook for the future in Chapter 6. 

We consign various technical details to the Appendices: We give the explicit 
formulae for the AB = 1 Wilson coefficients; the analytic evolution of the 
LCDA; a set of useful expressions from QCD factorisation, specifically the 
decay amplitudes for B — > tttt and the annihilation contributions to the B — > 
V'-y decays. We also present a full summary of numerical input parameters 
that have been used in our analyses. 
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Chapter 1 

Basic concepts of B-physics 



It's a job that's never started that takes longest to finish 

J.R.R. Tolkien 

The Standard Model is the cornerstone of particle physics and this chapter 
introduces its most pertinent features. We begin with a brief summary of 
the structure of the Standard Model and then introduce the concepts of 
effective theories that underpin calculations of B-decays. Via a discussion on 
the renormalisation group we give the full effective Hamiltonian for AB = 
1 decays. Finally, we discuss two important methods of calculating non- 
perturbative information - QCD sum rules and Lattice QCD. 

1.1 The Standard Model 

The Standard Model is the most successful and comprehensive theory of par- 
ticle interactions to be developed in modern times. The model conjoins the 
theories of strong and electroweak forces into a unified framework based upon 
gauge symmetries. The dynamics can be described by a single fundamental 
Lagrangian, constructed of contributions from the three sectors: Quantum 
Chromodynamics (QCD), electroweak interactions and the Higgs sector. The 
gauge structure of the Standard Model is summarised as 

SU(3)c®SU(2) L ®U(l) Y 

QCD is the theory of the strong interaction that describes the gauge inter- 
actions between quarks and gluons, and was first introduced in the 1970's 
[6-10]. This force acts upon "colour charge" and is based upon the gauge 
group SU(3)c- The eight generators of this group represent eight force carri- 
ers (the gluons) which communicate the force between coloured objects. The 
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quarks carry colour charge and most importantly so do the gluons due to 
the Yang-Mills or non-Abelian nature of the QCD gauge theory; this allows 
the gluons to interact with each other and is an essential ingredient in the 
asymptotic freedom of QCD. This property allows the perturbative treatment 
of the strong interactions at short distances - where the coupling constant a s 
becomes small. At long distances, i.e small energies, the coupling becomes 
large and there is a total confinement of the quarks into colourless hadrons. 
Figure 1.1 shows the running of the strong coupling with the energy scale. 




Figure 1.1: Running of a s {^) with energy scale [i and summary of experimental 
measurements. Data points are (in increasing fi) from r widths, T 
decays, deep inelastic scattering, JADE data, TRISTAN data and 
e + e~ event shapes [11]. 



The next sector of the Standard Model is that of the electroweak interactions, 
described by the Glashow-Salam- Weinberg model [12-14] and which is based 
on the gauge group SU (2) l (E)U(1)y- This theory is a unified treatment of the 
weak and electromagnetic interactions, where the full group is spontaneously 
broken to QED U(1)q. This is the Abelian theory of pure electromagnetic 
interactions, which has the photon as its propagator. The weak interaction 
controls all particle decays which do not proceed via the dominant strong 
or electromagnetic interactions. The weak decays of the 5-meson involve 
a change of flavour of the b quark, and are controlled by a charged current 
interaction. We discuss the dynamics of flavour in detail in the next section. 
The spontaneous symmetry breaking of the electroweak theory occurs via 
the non-zero vacuum expectation value of the Higgs field - a scalar isospin 
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doublet 




(1.1) 



The Higgs field has four degrees of freedom, three of which provide masses to 
the electroweak gauge bosons W + , W~ , Z°. The remaining degree of freedom 
is theorised to manifest itself as a massive scalar boson, but to date there has 
been no experimental detection to verify this claim. Data from electroweak 
precision tests suggest the Higgs mass ran, should be light, and a lower bound 
from LEP exists of ~ 114GeV [15]. 

It would seem that the Higgs is the last missing piece in the particle puzzle. 
However even if the Higgs is found, there are still a number of issues that 
need to be addressed: for example, there is yet to be a consistent method 
of incorporating gravity into the model, nor is there an inclusion of masses 
for the neutrinos and their recently discovered oscillations. We also must 
address the hierarchy problem that exists if the Higgs is discovered as light, 
instead of having a mass ~ 10 16 GeV as would be expected from calculation 
of quantum corrections to ra^. 

Many different scenarios have been developed to explain these issues and 
many more which are not elaborated upon here. However, before one can 
study extensions and revisions of the Standard Model, we must understand 
properly the Standard Model itself! This thesis is concerned with the theory 
behind decays of the 5-meson, and so we now introduce and discuss in detail 
the most relevant sector of the Standard Model - the dynamics of flavour. 

1.1.1 The flavour sector 

The fermions appear in the Standard Model in three generations, each gen- 
eration identical in all its quantum numbers differing only by the masses of 
the particles. Concerning the electroweak interactions, the quarks and lep- 
tons are spilt into left-handed doublets and a right-handed singlet under the 
gauge group SU(2) L . These contain the leptons: 




and the quarks: 




t R 
b R 



d R 



The existence and nature of the right-handed neutrino is still under question 
and is currently subject to a great amount of theoretical and experimental 
probing [16,17]. This does not affect our subset of non-leptonic B-decays so 
we retain the "standard", massless neutrino version of the Standard Model. 

The electroweak interactions are described by the following Lagrangian, which 
is made up of a charged current and a neutral current. 

Miit — L CC + L-NC 

= -4= Ut W+ ^ + J^W-"] - e J™A» 9 — [J°Z»] 

\/2 M M cos9 w M 

The neutral current part of the Lagrangian is made up of the neutral elec- 
tromagnetic and weak currents J" n and J° which are given in terms of the 
(electric) charge and isospin of the fermions. We have: 

j; m = QfhJ 

Jl = H[(H-2Q f sm 2 9 w )-lh 5 ]f 

summing over all flavours. The charged current in the quark sector is given 
by 



where the L subscript again represents the left-handed projector |(1 — 75) 
which reflects the vector - axial- vector (V — A) structure of the weak interac- 
tion. Vckm is the Cabbibo-Kobayashi-Maskawa (CKM) matrix [18, 19] and 
is a 3 x 3 unitary mixing matrix which rotates the mass eigenstates (d, s, b) 
into their weak eigenstates (<f, s', b'), and allows for transitions between the 
quark generations. The leptonic sector is described by an analogous mixing 
matrix which (in the absence of neutrino masses) is given by the unit matrix. 

Symbolically the CKM matrix is written as 

Vckm = V cd V cs V cb I (1.2) 




v ud 


v us 


v ub 


v cd 


Vcs 


v cb 


v td 


V ts 


v tb 



Unitarity ensures that there are no flavour changing neutral currents at tree 
level. The elements V^- are in general complex numbers which are restricted 
only by the unitarity condition - they are free parameters of the Standard 
Model and can only be determined by experiment. In general, an n x n 
unitary matrix is described by n 2 parameters. If this were to correspond to 
n quark doublets say, then the phases of each of the 2n quark states can be 
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redefined whilst leaving the Lagrangian invariant. Hence, V should contain 
n 2 — (2n — 1) real parameters. As an orthogonal n x n matrix can have only 
jn(n — 1) real parameters, then we will be left with n 2 — (2n— 1) — \n(n — 1) = 
\(n— 1) (rz — 2) independent residual phases in the quark mixing matrix. Thus, 
we see that Vckm must be parameterised by three independent angles and 
one complex phase; it is this phase which leads to a non-zero imaginary part 
of V and which is essential to describing CP violation in the Standard Model. 

The standard parameterisation of Vckm is written in terms of three angles 
% (h3 = 1)2,3) and a CP violating phase 5 [11]. The most useful pa- 
rameterisation we use in this work is the Wolfenstein parametrisation [20], 
where each element is expanded as a power series in the small parameter 
A = \V US \ 0.22. This reads to £>(A 3 ) as 

/ i - f A A\ 3 ( P - i V ) \ 

Vckm = -A 1 - f AX 2 (1.3) 

V AX 3 (1 -p-irj) -AX 2 1 J 

We can extend this to higher orders in A, and corrections up to order 0(X 5 ) 
alters only the element V td , transforming p ^ p = p(l — A 2 /2) and rj — > fj — 
77(1 — A 2 /2). The latest experimental determination of these parameters is 
included in a full summary of input parameters which is given in Appendix 
D. 



1.1.2 Unitarity triangle 

The unitarity of the CKM matrix gives six independent relations, each of 
which can be geometrically represented as triangles in an Argand diagram. 
The relation that is normally employed to represent unitarity in the B-system 
is 

v ud v: b + v cd v: b + v td v* = (1.4) 

A phase convention is chosen where V cd Y* b is real and we rescale the triangle 
by |V^V^| = AX 3 . Graphing this in the complex plane (p,r)) leads to a the 
unitarity triangle with co-ordinates of the vertices at (0,0) and (1,0) and the 
apex at (p,fj) as shown in Figure 1.2. 

The sides are expressed as 

\/p 2 + fj 2 



R u = 



Rt 



\ V udVub\ 

\VcdV; b \ 

\v*vs>\ 
\v cd vi\ 




Figure 1.2: The unitarity triangle. 



We introduce the shorthand X p = V p bV* d which is used extensively in phe- 
nomenological applications. The unitarity relation (1.4) is invariant under 
phase transformations, which means that the sides and angles of the tri- 
angle remain unchanged with a change of phase - hence they are physical 
observables. These, and the elements of Vckm have been subject to as many 
experimental determinations as possible in an attempt to over-constrain the 
parameters of the triangle and to test the Standard Model. The constraints 
can come from a wide variety of different decays and parameters: for example 
B° — ► J/i/jKs for sin 2/5; B — ► pir, B — > Tin for a... These constraints are 
neatly summarised graphically in Figure 1.3 1 . 




-1 -0.5 0.5 1 1.5 2 

P 



Figure 1.3: Constraints on the unitarity triangle from global CKM fit as per- 
formed by CKM Fitter Group [21]. 



Updated results and plots available at: http://ckmfitter.in2p3.fr. 



1.2 B decays and effective field theory 



There are three types of decays of the _B-meson, categorised by the final state 
decay products. Firstly there are the leptonic decays, such as B + — > £ + U£ and 
semi-leptonic decays such as B° — > D~l + vp. The decays B — > £ + i>e + X, with 
X as anything, accounts for around 10% of the total B decay rate. Finally, 
and most important in the context of this work, are the fully hadronic (non- 
leptonic) decays. The complication of calculating weak decays in QCD is 
illustrated in Figure 1.4. This demonstrates the non-trivial interplay between 
the strong and electroweak forces which determine the dynamics of the decay. 



d 




d 



Figure 1.4: Illustration of QCD effects in B — ► irn decay. 

The decay is characterised by several very different energy scales, from the 
mass of the PF-boson and the heavy quarks to the intrinsic QCD scale and 
the masses of the light quarks. We have the ordering 

m t , M w > m h)C > Aqcd > rn uAs 

QCD effects at short distances (i.e. that at higher energies) can be calculated 
perturbatively thanks to the asymptotic freedom of the theory. However, 
there are unavoidable long-distance dynamics from the confining of light- 
quarks into bound states - the hadrons. These are characterised by a typical 
hadronic scale of // ~ lGeV, a 8 (fx) is no longer small at this scale so we 
cannot use perturbation theory. This means that there are non-perturbative 
QCD interactions which enter into the calculation of the decay. 

To systematically disentangle the long and short distance contributions we 
make use of the operator product expansion (OPE) [22-24]. The basic idea is 
that any decay amplitude can be expanded schematically in terms of 1/M W , 
since My/ is much heavier than all of the other relevant momentum scales. 
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The products of the quark current operators that interact (via the W- ex- 
change) are expanded into a series of local operators Qi multiplied by a 
Wilson coefficient C{. These coefficients represent the strength that a given 
operator enters into the amplitude. Schematically we have 

A = C t (Mw/fjL, a s ) ■ (Qi) + O { P 2 /M^) (1.6) 

In this way we can define an effective weak Hamiltonian to describe weak 
interactions at low energies. In this effective theory the W boson and the top 
quark are removed as explicit degrees of freedom - i.e. they are "integrated 
out". We describe this as an effective field theory with rif "active" quarks, 
i.e. at the scale of mj,, we have 5 active flavours. As a practical example of 
this we can consider the basic (tree-level) W- exchange process of b — > duu, 
for which the OPE gives the amplitude: 

C M 2 

with C — 1 and the local operator Q = (duj y _ A (ub) v _ A This is represented 
diagrammatically in Figure 1.5. 




Figure 1.5: Tree-level diagram for b — > duu. 



The full weak (effective) Hamiltonian, including QCD and electroweak cor- 
rections has the following structure 

v i 

where the factor Vq KM denotes the CKM structure of the particular operator. 
The decay rate for a two-body non-leptonic decay B — > MiM 2 , is then simply 

r = — ^ — \(M 1 M 2 \'H c r\B)\ 2 (1.8) 
167rm B 

with S — 1/2 ii Mi and M 2 are identical or S — 1 otherwise. The Qi denote 
the relevant local operators which govern the particular decay in question; 
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they can be considered as effective point-like vertices. The Wilson coefficients 
are then seen as "coupling constants" of these effective vertices, summarising 
the contributions from physics at scales higher than fi. Operators of higher 
dimensions corresponding to the terms of O (p 2 /M^) can be neglected. 

In short, the OPE gives essentially a factorisation of short and long distance 
physics. The Wilson coefficients contain all the information about the short 
distance dynamics of the theory, i.e. that at energy scales greater or equal 
to ji. They depend intrinsically on the properties of the particles that have 
been integrated out of the effective theory, but not on the properties of the 
external particles. The factorisation implies that the coefficients are entirely 
independent of the external states, i.e, the Cj are the same for all amplitudes. 

The long-distance physics, that at energy scales lower than /i, is parame- 
terised purely by the process-dependent matrix elements of the local opera- 
tors. The renormalisation scale \i can be thought of as a "factorisation scale" 
at which the full contribution splits into the high and low energy parts. 

The matrix elements of the local operators are not easily calculated and as 
we have discussed, must contain a degree of non-perturbative information. 
Methods of determining these matrix elements using a systematic formalism 
are introduced in the next chapter. The Wilson coefficients however, can 
be fully calculated perturbatively by matching the full theory (with the W 
propagators) onto the effective theory; this ensures that the effective theory 
reproduces the corresponding amplitudes in the full theory. The steps to 
compute the C, are summarised as: 

• Compute full amplitude Af u ii with arbitrary external states 

• Compute the matrix elements (Qi) with same external states 

• Extract the Cj using expression (1.6) 

Both ultraviolet and infrared divergences occur in calculation of the ampli- 
tude Afuif, we discuss the removal of the UV-divergences in the next section. 
In the matching procedure, the IR-divergences are regulated by setting the 
momenta of the external quarks to p 2 ^ 0. From the point of view of the 
effective theory, all of the dependence on p 2 (representing the long-distance 
structure of the full amplitude) is contained in the matrix elements (Qi). 
Hence, the Wilson coefficients are free from this dependence - and so inde- 
pendent of the external states. For convenience, in the matching the external 
states are chosen to be all on-shell quarks (or all off-shell), but in general any 
arbitrary momentum configuration will work. 

This procedure will give the initial conditions for the Wilson coefficients at 
the matching scale, in our case \i = M w . We can then use the equations of 
Renormalisation Group Evolution (RGE) [25] to find the value at any desired 
scale. 



Q 



1.2.1 Renormalisation group evolution 



Put simply, renormalisation removes infinities from a theory. Feynman di- 
agrams with internal loops often give ultraviolet divergences, as the virtual 
particle running through the loop is integrated over all possible momenta 
(from zero to infinity). Renormalisation allows the isolation and removal of 
all of these infinities from any physical quantity [26]. We relate the bare 
(unphysical) parameters with a set of renormalised (physical) parameters - 
such as masses or coupling constants - and rewrite all of the observables we 
need in terms of the new physical quantities. We can then "hide" all of the 
divergences in redefinitions of the parameters in the theory Lagrangian. 

The procedure of renormalisation introduces a dependence on a dimension- 
ful parameter known as the renormalisation scale, /i. We can subsequently 
obtain the scale dependence of the renormalised parameters from the \i- 
independence of the bare ones. If we choose a set of parameters at a certain 
scale q (which gives g{q),m(q) etc.), then the set of all transformations that 
relate parameter sets with different values of q is known as the renormalisa- 
tion group. 

An an example, we can consider the coupling constant of QCD: a s (fj) = 
g 2 / (47r). The renormalisation group equation (RGE) for the running coupling 
is given by 

where the beta function, (3(g), is related to the renormalisation constant for 
the coupling. In QCD this is given by 

9 \ 2 Q , / 9 



with 



M = -9\[^) Po+{^) /?! + ..•> (1.10) 
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A) = y N c - -n f 

A = ^N 2 -™N c n f -2C F n f 

N 2 ~l 

N c is the number of colours, rif the number of active flavours and Cp 



The leading order solution for the coupling a s (fx) is then found via 



2N C 



a sW = " 7— \" (1-11) 



l + g^o)m(V) 



This equation can be re-expressed in terms of a mass scale A, which is the 
momentum scale at which the coupling becomes strong as q 2 is increased. 
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We have: 

A) Ma) 

From (1.11) we can see that the renormalisation group equations allow for 
the resummation of large logarithms, which could otherwise be problematic. 
The large logarithms can spoil the validity of the perturbative expansion, 
even if the value of a s is still small. The form of the correction terms to 
higher orders can be summarised in the following diagram [27], denoting for 
example, the "large log" at scale A = M w : L = \n/j,/M w . 

LL NLL 

a s — > a s L a s 

a 2 — > a 2 L 2 a 2 s L a 2 

a 3 . — > a 3 L 3 a 3 L 2 a 3 L a 3 

O(l) 0(a s ) 

The rows of this table correspond to the expansion in powers of a s from ordi- 
nary perturbation theory. This is no longer true in the presence of the large 
logarithms, but is resolved by resumming the terms (a s L) n to all orders in n. 
This re-organisation is obtained by solving the RGE equation. Expanding 
the leading order terms in (1.11) we have 

°° / '8 2 \ m 

a s (n) = a s {no)Y^ (j^(^o)hi^f ) (1.13) 

m=0 ^ 47F ^ / 

This sums logs of the form In /i 2 which can become large if /i /i. We 
then speak of "leading logarithmic order" (LL) and "next-to-leading loga- 
rithmic order" (NLL), although we carelessly use these synonymously with 
the terms LO and NLO. 



1.2.2 RGE for Wilson coefficients 

As discussed, the Wilson coefficients can be interpreted as effective cou- 
plings for the operators Qi of the effective Hamiltonian. These operators 
still have to be renormalised (or more specifically the operator matrix el- 
ements (Qi)^). This is done using renormalisation constants for each of 

1/2 

the four external quark fields Z q and a matrix Z^- which allows operators 
with equivalent quantum numbers to mix under renormalisation: we have 
(Qi)^ = Z~ 2 Zij(Qj) . Since the operators are always accompanied by their 
corresponding Wilson coefficients the operator renormalisation is in general 
entirely equivalent to the renormalisation of their "coupling constants" Cj. 
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We can therefore write the RGE for the Wilson coefficients as 



(1.14) 



where 7^ = 7 is the anomalous dimension matrix for the operators, defined 
via 

7 _i dZ kj 



1.15) 
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The anomalous dimension matrix (ADM) is itself given in terms of a pertur- 
bative expansion in a s 



lij 
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(0) , 

Tij + 



(1.16) 



We can then give the solution for the Wilson coefficients in terms of an 
evolution matrix Uij(/i, fj, ) 

The evolution matrix is given generally by 



At leading order, this reduces simply to 
t> ,0 Wo) = 



9O0) 



(1.17) 



a s {jj) ) 



= V 



a s (n) 



7 (Q) 

2/3 



v- 1 



(1.18) 
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where V is the matrix that diagonalises and 7^ is the vector (with 

elements 7^) containing all of the eigenvalues of the leading order ADM 



7 (o)_ 



7 g) = v-^v 



It is also possible to calculate the exponentiated matrix directly without 
diagonalisation. At the next-to-leading order, we find that the evolution 
matrix becomes a little more involved, and includes dependence on the NLO 
ADM 7 (1) . The solution is 



An 



t> (0 Wo) 



' _ a a (/x ) j 
An 



(1.19) 
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with 



j = vsv- 1 

= 5.V 0) — ^ 

3 " jJl IB 2 ofl i (°) (°) 

G = V- l ^ T V 

In the weak decays of 5-mesons the matching is performed at the scale Mw, 
so that both the top-quark and the W-boson are integrated out (removed as 
explicit degrees of freedom). This gives matching conditions and hence values 
for the Wilson coefficients at this scale. We can then use the procedures just 
outlined to evolve these down to the appropriate scale, e.g m b . We write 

C i { l i) = U{ l i,M w )C i {M w ) (1.20) 

with an expansion of the coefficients given to the same accuracy as the evo- 
lution, i.e. to NLO 

Ci(M w ) = Cf\M w ) + ^Cf»(M w ) (1.21) 

The scale can be any required, e.g fi = mj for i?-decay amplitudes or /j, ~ 
lGeV for discussions on the wavefunctions of light mesons. Care must be 
taken to include the effects of the flavour thresholds at lower energies, which 
can be done simply, by applying the evolution equations in two stages with 
the correct number of active quark flavours in each stage. 

1.2.3 The AB = 1 effective Hamiltonian 

We consider a basis for computing non-leptonic B decays with change in 
beauty of AB = 1. The expressions below are for those decays with un- 
changing strangeness and charm, AS = AC = 0, but can easily be adapted 
to decays with AS = 1 by replacing d — > s. 

We begin with the tree-level process without QCD corrections, which is de- 
scribed by one dimension 6 operator. When we include QCD corrections, 
another current- current operator is generated: these are labelled as Q\ and 
Q2, although sometimes in the literature their definitions are interchanged. 
The tree-level diagram and the 0{ot s ) corrections to it are shown in Figure 
1.6. 

QCD corrections also produce four new gluonic penguin operators Q3 to 
Qq. If we include terms from the electroweak sector up to 0{a) we obtain 
an additional set of electroweak penguin operators Q 7 to Q w . These are 
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Figure 1.6: Tree-level exchange and 0(a s ) corrections; p = u,c. 



considered as a next-to-leading order effect due to their proportionality to 
the electroweak gauge coupling a. There are, in principle, QED corrections to 
the matrix elements of the QCD operators Q\ . . . Qq, but these are suppressed 
and usually neglected. We also have additional terms which contribute in 
some AB = 1 processes, namely the magnetic penguin operator Qj^ (which 
is important for the radiative decays) and the chromo magnetic penguin Q§ g . 
Examples of all these diagrams are shown below in Figures 1.7 and 1.8. 






Figure 1.7: Gluonic and electroweak penguin diagrams. 





Figure 1.8: Magnetic photon and chromomagnetic penguin diagrams. 



This leads us to a AB = 1 effective Hamiltonian of 



G 



p=u,c 



G\Q\ + C2Q2 + CiQi + Cj^Qj-f + C$gQ£ 



i=3,...,10 



+ h.c. 
(1.22) 
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with the full operator basis given as 

Qi = (diPj) v _ A (Pjbi) v _ A 
Ql = {dp) v _ A {pb) v _ A 

Qs = {db) v _ A J2(^)v-A 
i 

Q 4 = {d i b j ) v _ A ^{q j q i ) v _ A 
q 

Qs = {db) v _ A J2^v+A 
q 

Q 6 = {d i b j ) v _ A ^{q j q i ) v+A 
q 

q 

Qs = (sibj)y_ A J2hq(mi)v+A 
q 

Q9 = (sb) v _ A ^\e q (qq)y„ A 

q 

Q10 = (sjbj) v _ A ^ \e q (q~jqi) v _ A 
q 

= ^m 6 scr^(l + 75)F^fe 
Q 8g = JL mb s i<T i» (1 + ^)^0% 

We use the usual notation that (qiq2)y± A = ?i7^(l ± 75)^2; h j are colour 
indicies, F^ u and are the photonic and gluonic field strength tensors re- 
spectively. There is also implied summation over all flavours q. The matching 
conditions for the Wilson coefficients and the anomalous dimension matrices 
required for scale evolution are given in detail in Appendix A. 



1.3 QCD sum rules 

In later chapters we make use of a number of results from the technique 
to study hadronic structure known as QCD sum rules [28]. Without any 
detailed technical discussion, we take a moment here to briefly describe the 
purpose and form of the sum rule approach. They were originally used to 
calculate simple characteristics of hadrons, such as masses, but are also ap- 
plicable to more complicated parameters such as form factors or hadronic 
wave functions. The sum rule approach is used to calculate non-perturbative 
effects in QCD. 
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The hadrons are represented by interpolating quark currents, which are 
formed into correlation functions. These are treated in an operator prod- 
uct expansion to separate the short and long distance contributions from 
the quark-gluon interactions. The short distance parts are again calculable 
in normal perturbative QCD. The long-distance contributions are parame- 
terised by universal vacuum condensates or by light-cone distribution ampli- 
tudes. (See later chapters). The results of this calculation are then matched 
via rigorous dispersion relations to a sum over the hadronic states. This gives 
us a sum rule which allows us to calculate the observable characteristics of 
hadronic ground states. Together with experimental data this can be used 
to determine quark masses and universal non-perturbative parameters. 

There are some limitations in the accuracy of the sum rules which come 
from the approximation of the correlation functions. Additionally there is 
uncertainty from the dispersion integrals, which have complicated and often 
unknown structure. The method is however very successful and all uncer- 
tainties can be traced and estimated well. 

The classical sum rule approach is based on two-point correlators, but this 
is not the only approach. Combining sum rule techniques with a light-cone 
expansion gives the very successful technique of light- cone sum rules [29,30]. 
This procedure expands the products of currents near the light-cone and 
involves a partial resummation of local operators. 

1.4 Lattice QCD 

Ideally we would like to be able to estimate all of the fundamental parameters 
of QCD from first principles. To do this, we need complete quantitative 
control over both the perturbative and non-perturbative aspects of QCD. 
The numerical simulations of Lattice QCD in principle are able to do this 
- answering almost any question about QCD and confinement, from Yang- 
Mills theories to the running coupling constant to the topology of the QCD 
vacuum. In our work, the relevant parameters that can be determined via 
lattice techniques include decay constants [31,32], form factors [33,34] and 
hadronic matrix elements [32]. As with the results we frequently rely upon 
from QCD sum rules, it is important to understand the importance and 
limitations of lattice results that we use; a full review can be found in [35]. 

QCD can be expressed in Feynman path integrals, which need to be cal- 
culated in the continuous space-time in which we exist. This is very diffi- 
cult but we can make an approximation by discretising the continuous four- 
dimensional space-time onto a 4-D lattice as represented in Figure 1.9. The 
Lagrangian consisting of fields and derivatives of fields is also discretised by 
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replacing the continuum fields with fields at the lattice sites, and replac- 
ing the derivatives with finite differences of these fields. The path integrals 
over different field configurations become mult i- dimensional integrals over 
the values of the fields at the lattice sites and on the links. These can be 
evaluated numerically, for example using Monte Carlo simulation; this is the 
basic tenet behind lattice QCD. The problem of solving a non-perturbative 
relativistic quantum field theory in transferred into a "simple" matter of 
numerical integration. 

node 



link 



a 

Figure 1.9: Nodes are separated by a lattice spacing of a. Quark fields live on the 
lattice sites (nodes) and gluons on the links connecting these sites. 

Large amounts of computing power are required for simulations, and the 
computation time depends on a number of factors - the overall volume of the 
spacetime L 4 where L is often taken as L ~ 1 — 2fm, and the finiteness of the 
grid, a ~ 0.05 — O.lfm. Ideally the grid spacing should be small enough that 
results are within a few percent of the continuum results. Errors associated 
with the discretisation are an important source of systematic errors in a 
lattice calculation that need to be controlled. 

The purely gluonic part of QCD can be simply discretised onto the lattice 
but complications arise for the quarks due to their fermionic nature. There 
are however a number of different methods of including quarks on the lattice 
such as Wilson quarks or staggered quarks, each with its own improvements 
and disadvantages. Implementing an algorithm for quarks gives an additional 
problem as fermionic variables in path integrals require the calculation of a 
fermionic determinant - requiring hundreds of matrix inversions for every 
vacuum configuration considered. This calculation amounts to determining 
the non-local interaction between the gluons, i.e the computation of dy- 
namical quark loops. There is no technical obstruction to performing this 
calculation and a number of methods exist for the discretisation of fermions, 
however the computation time and difficulty of the calculation increases by 
several orders of magnitude! 

This has led to many simulations being performed in the quenched approx- 
imation, where vacuum configurations that include only gluons are consid- 
ered. This neglects the additional term which gives rise to the sea quarks 
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(and hence the ability to produce qq pairs from the vacuum). The quenched 
approximation performs satisfactorily for processes dominated by valence 
quarks, but for the converse dynamical quarks are required to produce correct 
answers. Unquenching can have other side effects such as allowing processes 
which require decays into qq pairs e.g. g — > qq — > g which will enable the 
strong coupling to run at the correct rate. This also means that decays such 
as p — > 7T7T are now visible on the lattice, but it will be difficult to determine 
m p . The cost of including the dynamical quarks increases as we move to 
lighter and lighter quarks, and it is the lightest three u, d, s that have the 
most significant effect on most quantities. 

Heavy quarks, such as the b quark could in principle be treated the same as 
the lighter quarks, however in order to gain sufficient accuracy an extremely 
fine lattice would be needed. Instead, exploiting the non-relativistic nature 
of the b quark bound states proves much more efficient. This can be done 
using several techniques including static quarks - using the heavy-quark limit 
of m h = oo and NRQCD - a non-relativistic formulation of QCD. 
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Chapter 2 

QCD factorisation 



Still round the corner there may wait, 
A new road or a secret gate. 

"The Lord of the Rings", J.R.R. Tolkien 

The concept of factorisation in its various forms is key to many of the ap- 
plications of perturbative QCD, including deep inelastic scattering or jet 
production in hadron colliders. Perturbative QCD only describes quarks and 
gluons, so in observed, hadronic physics there will always be some relevant 
non-perturbative dynamics that must be isolated and dealt with in a system- 
atic fashion. In the following sections, we describe the ideas of "naive fac- 
torisation" for the hadronic decays of heavy mesons, and then introduce the 
"improved QCD factorisation" formalism - a systematic method for treating 
exclusive two-body B decays. There have been a number of other approaches 
to try to describe non-leptonic B decays, although none have been nearly as 
successful as QCD factorisation. These methods include the perturbative 
QCD (pQCD) methods (which treats the form factors as perturbatively cal- 
culable quantities) or the soft-collinear effective theory (SCET), which is 
based upon factorisation of two energy scales: m\ ^> Aqcd^& ^ Aqcd- 
QCD factorisation is currently the accepted and most widely used methodol- 
ogy for treating B decays, and we work exclusively in this framework in this 
thesis. Finally, we conclude this chapter with a discussion of the limitations 
of the QCD factorisation method and ask how well predictions match up 
with the experimental data. 

2.1 Naive factorisation 

A complete understanding of the theoretical framework for non-leptonic B 
decays is a continuing challenge which is yet to be resolved. For two-body 
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decays of B — > MiM 2 , the major difficulty involves the evaluation of the 
hadronic matrix elements (M 1 M 2 \Qi\B). 

We introduced the concept of factorisation in relation to the OPE where 
we separate the long and short-distance contributions to decay amplitudes. 
Similarly, the hadronic matrix elements for B decays can also be factorised 
by disentangling the long and short-distance contributions. This seems in- 
tuitive for leptonic and semi-leptonic decays, where we can factorise the 
amplitude into a leptonic current and the matrix element of a quark current. 
Gluons do not interact with the leptonic current and so the factorisation 
is exact. For non-leptonic decays, the full matrix element is separated into 
a product of matrix elements of two quark currents; so we must also have 
"non-factorisable" contributions from gluons which can connect these two 
currents. 

The first general method proposed is that of naive factorisation [36,37]. 
Here the full matrix element is assumed to factorise into a product of matrix 
elements of simple, colour singlet, bilinear currents. The first is the transition 
matrix element between the P-meson and one of the final state mesons, 
and the second is the matrix element of the other final state meson being 
"created" from the vacuum. For example 



(n+n-\(ub) v _ A (du) v _ A \B°) — (n-\(du) v „ A \0)(n + \(ub) v _ A \B°) (2.1) 



These matrix elements can then be expressed in terms of a form factor 
F B ^ Ml and decay constant fu 2 ■ The Hamiltonian is re-expressed in terms of 
effective coefficients C^ s . These are constructed from the Wilson coefficients 
and all of the scheme and scale dependent contributions from the hadronic 
matrix elements. This ensures the effective coefficients are free from these 
dependencies. The complete matrix element is then made up of the fac- 
torised matrix elements and the coefficients Oj, which are combinations of 
the effective Wilson coefficients 



where Pj involves the QCD and electroweak penguins (but not taking into 
account final state flavour) and are zero for % — 3, 5. 

A qualitative justification for this approach comes from the concept of colour 
transparency [38, 39] which implies the decoupling of soft (low-energy) glu- 
ons from colour-singlet pairs of quarks, e.g. the final state pion in example 
(2.1). This occurs as the heavy-quark decays are highly energetic, so the 
hadronisation of the qq pair occurs far away from the remaining quarks. 

The main source of uncertainty in this approach is the neglect of the "non- 
factorisable" contributions, exacerbated by the assumption that all final state 
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interactions are absent. Corrections arise from the exchange of gluons be- 
tween the two quark currents, which allows a dynamical mechanism (namely 
re-scattering in the final state) to create a strong phase between different am- 
plitudes. The second major issue is that the cancellation of renormalisation 
scale dependence in the amplitude is destroyed, as the scale independence of 
the form factor and decay constants in the factorised matrix element is in 
conflict with the scale dependence of the original matrix element. Unphys- 
ical dependencies are also caused by the Wilson coefficients at NLL level 
which develop scheme dependence in addition to the scale dependence, and 
for which there is no cancellation from the (scheme independent) factorised 
matrix elements. 

2.2 QCD factorisation 

QCD factorisation was introduced by Beneke, Buchalla, Neubert and Sachra- 
jda (BBNS) [1,2,40], and is based upon the important simplifications that 
occur in the limit where the b quark mass is large as compared with the 
strong interaction scale, m& ^> Aqcd- The factorisation that occurs in this 
approach is the separation of the long-distance dynamics (the matrix ele- 
ments) and the short-distance interactions which depend only on the large 
scale m b . The short-distance contributions are calculated perturbatively to 
order a s , and the long-distance information is encoded in various process 
independent non-perturbative parameters, or is obtained directly from ex- 
periment. However, it is structurally much simpler than the original matrix 
elements. 

In the limit mj ^> Aqcd, the underlying physics involves the decoupling 
of fast-moving light mesons, produced from point-like interactions from the 
weak effective Hamiltonian, from the soft QCD interactions. This is a result 
of the concept of colour transparency as discussed above, the systematic 
implementation of which is provided by QCD factorisation. It also provides 
a soft factorisation which allows the matrix elements to be given at leading 
order in the Aqcd I expansion. The full matrix elements can then be 
represented in the form 

{M l M 2 \Q,\B) = (M 1 \j 1 \B)(M 2 \j 2 \0) 

where r n denote the radiative corrections in a s , and ji are bilinear quark 
currents. If the order a s corrections are neglected, we see that at leading 
order in Aqcd/^6 we recover the naive factorisation results. 

In the QCD factorisation approach, the non-factorisable power-suppressed 
corrections are in general neglected, with two non-trivial exceptions: the 



1 + r «< + ° ( A QCD/m b ) (2.3) 
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hard-scattering spectator interactions and annihilation contributions which 
are chirally enhanced and cannot be ignored. In the context of this thesis, 
it is important to note that these cannot be calculated within the actual 
framework of QCD factorisation; they are included in the BBNS approach 
via model-dependent assumptions. This will be discussed in detail in the 
following sections. 

2.2.1 Structure of the QCD factorisation formula 

The QCD factorisation framework is represented by a "master formula" 
for the matrix element (MiM 2 \Qi\B) , where the final state mesons can be 
"heavy-light" (e.g. B — > Du) or "light-light" (such as B — > tttt). For exclusive 
non-leptonic decays to two light mesons, we have the following expression for 
(M 1 M 2 \Qi\B) , to leading order in the Kq C D/ m b expansion: 

(M x M 2 \Qi\B) = F?^{ml) f duT^u) <S> M2 (u) + [M ± M 2 ] 

Jo 

+ [ d£ du dv T^{i, u, v) $ B (0 $Af! (v) $m 2 (u) (2.4) 
Jo 

$m denote the light cone distribution amplitudes (LCDA) for the valence 
quark states. Both the LCDA $m and the B — > M form factor, are much 
simpler than the original non-leptonic matrix element, and can be calculated 
by some non-perturbative technique such as QCD sum rules, on the lattice, 
or taken directly from experiment. 

The perturbative information is contained in the hard-scattering kernels 
T[j(u), Tf k (v) and T/ / (^,-u,f), which are calculable functions dependent on 
the light-cone momentum fractions of the constituent quarks in the light final 
state mesons (u, v) and the B meson (£). All of the hard-scattering kernels 
and the LCDA have a factorisation scale dependence. The expression is split 
into two categories of contributions - "type I" and "type-II" . The type-I con- 
tributions consist of the leading terms which give the tree level contributions, 
and hard vertex and penguin corrections at 0(a s ). The type II contributions 
originate from hard interactions between the spectator quark and the emit- 
ted meson, which enter at 0(a s ). If the spectator quark in the interaction 
can only form into one of the final state mesons (such as in B° — > tt + K~), 
then the second form factor term is absent. 

We can see that from the tree level terms in T 1 we reproduce the naive 
factorisation results, and the convolution integral in (2.4) will reduce to a 
meson decay constant. We also get a large simplification in the case where 
one of the mesons is heavy (B — > HiM 2 ), as the spectator interactions will 
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be power-suppressed in the heavy quark limit and the type II kernel T is 
absent. 



2.2.2 Non-perturbative parameters 

We will now elaborate further on the non-perturbative input into the factori- 
sation master formula, and introduce the concept of light-cone distribution 
amplitudes which will be discussed in detail in Chapter 3. 



Light-cone distribution amplitudes for light mesons 

The distribution amplitude §M(u,q 2 ) is in essence a probability amplitude 
for a meson to be found in a particular state, i.e for finding a valence quark 
with a light- cone longitudinal momentum fraction u in the meson at a mo- 
mentum q 2 , independent of the process. The light-cone frame for some vector 
Pn = ipo, Pi, Pi-, P3), is a choice of co-ordinates which naturally distinguishes 
between the transverse and longitudinal degrees of freedom. Two of the spa- 
tial dimensions (transverse) remain unchanged and the other (longitudinal) 
spatial dimension is combined with the temporal co-ordinate via 

Po ± Ps fo _v 

P± = — (2.5) 

so that p^ = (p + , p-, p±), and p 2 = for a light-like vector. The distribution 
amplitude is defined in terms of the expectation value of non-local operators, 
near the light-cone. For example, considering the pion distribution amplitude 
<f) n , we express the LCD A in terms of a matrix element of a gauge-invariant 
non-local operator defined between the vacuum and the pion [41] 

(0\d(z)^ l5 [z,0}u(0)\7r(P)) (2.6) 

[x, y] is a path-ordered gauge factor along the straight line connecting the 
points x, y in order to preserve gauge invariance and is known as a Wilson 
line 



[x, y] = Pexp 



ig 



[ dt(x-y)^(tx + (l-t)y) 
Jo 



(2.7) 



These matrix elements give a set of integral equations of two and three par- 
ticle light-cone distribution amplitudes. These distribution amplitudes can 
be classified in terms of twist. The twist of an operator is defined as its 
"dimension minus spin", i.e. 

t = d-s (2.8) 

This is related to power counting of 1/Q, where Q denotes momentum trans- 
fer, and controls the relative size of contributions from the operator product 
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expansion. In general, an operator of dimension d in the OPE has a coeffi- 
cient function with dimension (mass) 6 ~ d [42], corresponding to a suppression 
factor in the Fourier transform of the OPE of 



Q 



[2.9) 



If the operator has spin s, the operator matrix element gains s contributions 
from the momentum P M , so the total contribution is of the order of 



21' -qY ( 1 N '' * 2 



Q 2 J \Q 



(2.10) 



with q 2 = —Q 2 . The leading twist is then t — 2, allowing us to identify a 
distribution amplitude of twist- 2 and also of higher twists. The higher-twist 
distribution amplitudes contain contributions from lower twists as well. For 
example the twist-3 DA has contributions from both twist-3 and twist-2 
operators. In the QCD factorisation formalism, and in Chapter 3, where we 
study models of LCDAs, we consider the simplest case of the leading-twist 
distribution amplitudes. From the matrix element in (2.6) the twist-2 DA is 
defined as 

(0|J(z) 7 M 75MM0)7r(P))| 22=0 = if n pn I' due iu ^ P Un(u,q 2 ) (2.11) 



A detailed discussion of light-cone distribution amplitudes is continued in 
Chapter 3. 



Light-cone distribution amplitude for B-meson 

Unlike the distribution amplitudes for the light mesons, there is very limited 
knowledge of the parameters determining the .B-meson wavefunction. In the 
context of QCD factorisation, it is the inverse moment of the distribution 
amplitude that is most relevant (see below). For scales much larger than 
nib, (ft b should tend to a symmetric form, as with a light meson distribution 
amplitude. However, at or below m^, the distribution is expected to be very 
asymmetric in momentum fraction £, with £ ~ 0(AqcD/ m b)- 

The .B-meson light-cone distribution amplitude only appears in the type-II 
(hard spectator interaction) term in the main QCD factorisation formula. 
The amplitudes for these interactions, at order a s , depend only on the prod- 
uct of the momenta of the light meson which absorbs the spectator quark, 
p', and that of the spectator quark itself, /. Converting to light-cone co- 
ordinates we can arrange the momenta so that the hard-spectator amplitude 
is dependent only on l + . The wavefunction can then be integrated over the 



other momenta l± and Z_. Following the method of [2], the hard-spectator 
integrals are re-expressed in terms of the longitudinal momentum fraction £, 
and the 5-meson LCDA can be decomposed (at leading order in 1/ mj) into 
two scalar wavef unctions. These wavefunctions describe the distribution of 
the longitudinal momentum fraction £ = l+/p+ of the spectator quark inside 
the meson and are defined via: 

(Q\q tt (Q)b p (z)\B(p)) = l -lf[(p'+rn b ) l5 }^J 1 d{;e- i ^ z - [$ B1 (0+^- ^(flU 

(2.12) 

where n is an arbitrary light-like vector, which is chosen in the direction of 
one of the final state momenta, n_ = (1,0,0,-1). The wavefunctions are 
normalised 

/ de^Bi(o = i / ^ 2 (o = o (2.i3) 

In the calculations involving the _B-meson distribution amplitude, we need 
only the first inverse moment of the wavefunction $bi(£), parameterised as 

f d ^ = <M (2 . 14) 

There is little information as to the actual value of the parameter X B . There 
is a known upper bound which implies 3\b < 4A [43], where A = tub — mj, 
corresponding to \b ~ 600MeV. Numerical estimates have been suggested 
from a number of different models [43-45] which we can combine to give an 
estimate of around \b = 350 ± 150MeV. 



Form factor 

Hadronic form factors describe the inner structure of the hadron, and are 
functions of scalar variables arising from the decomposition of some matrix 
element. In the coupling of a particle to a photon, we have an electromagnetic 
form factor, which is a momentum dependent function reflecting the charge 
and magnetic moment distribution, and hence the internal structure. In 
this case we are interested in the transition form factor, which describes 
the overlap of the B meson and the decay product (i.e. some pseudoscalar 
meson) during the actual decay. Decays of B — > ix are fully described by 
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three form factors F + , F and F T , which are defined by [46] 



(n(p)\u^b\B(p B )) = F + {q 2 ) j(p fl + p), - + ™*~ ^ F (q 2 )q, 

(2.15) 

(n(p)\da^(l + l5 )b\B(p B )) = i {( PB +P),q 2 - q,(m 2 B - ml)} ' iW) 



m B + m^ 
(2.16) 



where q = Pb — P, q 2 = rn B — 2mBE 7T . In QCD factorisation the vector 
current arises most often. We should also note that at the scale q 2 = the 
two form factors coincide, F + (0) = -Fo(O). The asymptotic scaling behaviour 
of the form factors at q 2 = is given as [47] 

F+(0) = F (0) ~ (2.17) 
m b 



The form factors can be calculated using various non-perturbative techniques, 
such as QCD sum rules (as we discussed in Section 1.3), and on the lattice 
[33,34]. The method of light-cone sum rules [46] allows the calculation of the 
form factor within a controlled approximation, relying on the factorisation 
of an unphysical correlation function whose imaginary part is related to the 
form factor in question. This correlation function is that of a weak current 
and a current with the same quantum numbers as the 5-meson, evaluated 
between the vacuum and the n; it is related to the form factor via 

U,(q,p B ) = i f d*xe^(7t(p)\TV,(x)j B (0)\0) (2.18) 
= Tl + (q 2 ,p 2 B)(p + PB)» + n(q 2 ,p 2 B )q, (2.19) 

where j B = im h d^b. When p 2 B m 2 , these correlation functions can be 
expanded on the light-cone 

ni c (g 2 ,p|) = / 1 ^Ti" ) ( M ,g 2 ,j 9 |, / i)0(«)K/i) (2.20) 

n J ° 

This sum runs over contributions from the pion distribution amplitude in 
increasing twist, where n = 2 is the leading- twist contribution. The functions 
T± are hard-scattering amplitudes, determined by a perturbative series in a s , 
known to 0(a s ) at twist-2 and twist-3. The form factor can then be extracted 
via a light-cone sum rule which depends on the spectral density p^ c of the 
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correlation 11+ 

e-</ M2 m 2 B f B F^(q 2 )= dse~ s ' M2 pl C {s,q 2 ) (2.21) 

Jml 

where M 2 and s are specific parameters from the sum rule, related to a 
Borel transformation [46]. 

2.2.3 Contributions to hard-scattering kernels 

It is useful to discuss qualitatively the diagrams that contribute to the hard- 
scattering kernels T^{u) and T^^^u, v) at the leading order and with ex- 
change of one gluon. We use the decay of B^ — > 7r + 7r~ as an illustrative 
example. The factorisation formula of equation (2.4) has been proved to 
one-loop for decays into two light mesons and to two loops for decays into 
heavy- light final states. There is also proof to all orders for B — > Dir [48]. 

The leading order diagram represents the quark level process b — > uud, and 
there is only one contributing diagram with no hard gluon interactions, as 
shown in Figure 2.1. Since the spectator quark is soft and does not undergo 



Figure 2.1: Leading order contribution to the hard-scattering kernel T 1 . 

a hard interaction, it is absorbed by the recoiling meson, and is described 
by the B — > n form factor. In the heavy-quark limit, we can represent the 
scaling of the decay amplitudes as [40] 

A(B° -> 7T+7T-) ~ G F m 2 b F B ^(0)f n (2.22) 

Radiative corrections to this diagram are suppressed by a power of a s or 
A-QCD/mb, or are already accounted for by the definition of F 3 ^ 71 or /„-; 
gluon exchange diagrams that do not fall into the above category are those 
that are "non-factorisable" with respect to the naive factorisation formalism. 
The hard-scattering kernel T 1 contains contributions from diagrams of this 
type, including vertex corrections, penguin contractions and contributions 
from the chromomagnetic dipole operator. These are illustrated in Figures 
2.2 and 2.3. 

The type-II hard-scattering kernel T 11 contains the hard spectator interac- 
tions, of the type shown in Figure 2.4. These would violate factorisation if 
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A/ V >>C \A 



Figure 2.2: The "non-factorisable" vertex corrections. 



Figure 2.3: Contributions from penguin dipole operator (left) and chromomag- 
netic dipole operator (right). 



there was soft-gluon exchange at leading order, however they are suppressed 
because of the endpoint suppression of the light-cone distribution amplitude 
for the recoiling ir. 




Figure 2.4: Hard-spectator contributions to kernel T . 



2.3 Basic formulae for charmless B-decays 

It is not necessary to reproduce all of the numerous expressions of which the 
QCD factorisation framework is composed, as they are discussed at length 
in the literature. On the other hand, an exposition of the most important 
of these formulae which make up the backbone to the QCD factorisation 
method is illuminating, and will be outlined here. We highlight the places 
where model- dependent contributions arise, and discuss the methods used by 
BBNS to parameterise them (Section 2.3.2). Our in-depth analysis of non- 
factorisable contributions to non-leptonic decays is found in Chapter 4. The 
full decay amplitude expressions as expressed in [2] for B — > -kit are given for 
completeness in Appendix B. 

2.3.1 Factorisable contributions 

We have introduced the non-perturbative input for the factorisation formula 
(2.4); what remains is to discuss the perturbative part. We will concentrate 
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on the factorisable part - that which is fully calculable in the factorisation 
framework. This is done by translating the effective weak Hamiltonian into 
a transition operator so that the matrix element is expressed, for example, 

as 

(7m\n eff \B) = % E M^WB) (2.23) 

p=u,c 

The transition operator T p is constructed from operators labelled by the 
flavour composition of the final state, corresponding to the operators of the 
weak Hamiltonian. These are multiplied by a QCD factorisation coefficient, 
af(7T7r), which contains all of the information from evaluation of diagrams 
corresponding to that topology. 

For example, the coefficients ai(7T7r) and a2(7T7r) are related to the current- 
current operators and a^im) . . . a§{itit) are related to the QCD penguin op- 
erators. The coefficients a 7 {jm) . . . ai (7r7r) are partially induced by the elec- 
troweak penguin operators Q 7 . . . Q 10 and are of order 0(a). These coeffi- 
cients do not include any annihilation topologies - they are model-dependent, 
power suppressed corrections and as such are evaluated separately, as dis- 
cussed in the next section. 

The general form of the factorisation coefficients at the next-to-leading order 
in a s is 



+ 



Ci±\ Cf 



N c 



+ P?(M 2 ) 
(2.24) 



where the upper signs apply for i odd and lower for % even, and the superscript 
p again runs over p = u,c (except for i — 1,2 where there is no flavour 
dependence and p is omitted). The Wilson coefficients Cj = Cj(/x) are taken 
at next-to-leading order in a s . The quantities that enter these coefficients 
are split up as follows: 

Leading order - the coefficient Ni(M%) describes the normalisation of the 
relevant LCDA, and is unity in all cases except % — 6, 8. 

Vertex contributions - V^(M 2 ), originate from processes as in Figure 2.2. 
These are calculated via convolution integrals of the LCDA for the meson 
M 2 with a scalar hard-scattering function (this function takes different forms 
depending on i). The pion LCDA 0^ is used for B — > nir decays, while 
both <f) n and <pK are used in B — > -kK decays. For all % except % = 6, 8, 
the leading twist distribution amplitude is used; the i — 6,8 terms use the 
twist-3 distribution amplitude, which is $ 3 (x) = 1 for pseudoscalar mesons. 
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Penguin contributions - Pj(M 2 ), as in Figure 2.3. These are present at 
order a s only for % = 4, 6. Again, these are convolutions of a light-cone 
distribution amplitude with a hard-scattering function, where the function is 
dependent on the internal quark mass in the penguin diagrams. This enters 
as z p , where z u = and z c = m 2 c jm\. 

Hard spectator terms - Hi(M\M2), as in Figure 2.4. These terms originate 
from the hard gluon exchange between the final state meson (M 2 ) and the 
spectator quark, and make up the kernel T n . They are calculated at a lower 
scale than the other contributions which are evaluated at fj, — m h . For these 
terms we associate a lower "hard" scale //^ ~ (AgcD^b) 1 ^ 2 , where we use 
Aqcl> ~ Ah = 0.5GeV. We also evaluate the relevant Wilson coefficients in 
(2.24) at the same scale. The contributions are again convolution integrals 
of the distribution amplitudes of the light, final state mesons and that of the 
£>-mcson. The light LCDA enter both at leading twist, and in the chirally 
enhanced twist-3 contributions. Only the twist-2 contribution is dominated 
by the hard gluon exchange and so is calculable, the twist-3 terms having 
logarithmic endpoint singularities. 

The coefficients a6(MiM 2 ) and as(MiM 2 ) are power suppressed by a ratio 
r^ 1 which is proportional to the QCD quark condensate - these are therefore 
referred to as chirally enhanced corrections. The suppression ratio is 

2m 2 

x m b (n)(m u (n) + m d {fi)) 

for B — > 7T7T decays. Although this factor is formally suppressed by AQCD/m>b, 
it always appears in conjunction with the coefficients a^, as, and can cause a 
large enhancement of power-suppressed corrections, as we consider next. 



2.3.2 Power-suppressed corrections 



The first main source of non-factorisable corrections to QCD factorisation is 
through the chirally-enhanced power corrections that are identified with the 
endpoint divergences arising in the hard spectator terms. These terms are 
calculated via 



H i (M 1 M 2 ) 



/b/tt 



771 



1^(0) I 



[ dx [ 
Jo Jo 



dy 



$M 2 {x)$ Ml (y) 



xy 



_l_ r 



<£> M2 (x)$ 3 mi 



xy 



(2.26) 



where x = (1 — x). This simplifies, for example for B — > nn, to 



<i>n(x)<i>M ! r n M^)^iy) 

xy x xy 

(2.27) 



since the asymptotic twist-3 amplitude 0^ = 1 (for all pseudoscalar mesons), 
and the pion distribution amplitude is symmetric in exchange x <-> x. The 
first integral in the bracket in (2.27) is finite. However the second, involving 
the twist-3 amplitude (where the explicit power suppression can be seen), is 
logarithmically divergent. BBNS introduce a model-dependent parameteri- 
sation of this divergence by defining 

-m _ f 1 d y is t„ .\ f 1 d y 



This parameter represents the soft-gluon interactions with the spectator 
quark and is expected to be of the order Xff ~ In fj, b / Aq C d- It is also 
treated as universal, i.e. the same value for all cases, so there is no depen- 
dence on M or the index % of the contribution in which it arises. It is written 
in the form 

X H = {l + p H e^)\^ (2.29) 

Aft 

The second source of non-factorisable contributions to QCD factorisation is 
the evaluation of the annihilation topologies as in Figure 2.5. 



Figure 2.5: Power suppressed annihilation contributions. 



Their contribution can be numerically significant, especially in B — > ttK. 
The calculational difficulty arises in the appearance of endpoint singulari- 
ties even for the leading twist contributions, which cannot be dealt with in 
the hard-scattering formalism. Since there is no consistent way of treating 
these contributions, an approximated model ignoring these soft endpoint di- 
vergences is used on top of the QCD factorisation setup in order to estimate 
the size of the annihilation contributions. The contributions are added at 
the amplitude level, via a new transition operator T® nn giving a new set of 
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coefficients, bi(MiM 2 ). 



{™\Keff\B) = J= J2 M^\T P + T£ nn \B) (2.30) 

* p=u,c 

The divergences that appear are treated in a similar way to the power cor- 
rections for the hard-spectator diagrams, and are written in terms of an 
analogous complex parameter Xa, introduced via 




(2.31) 



where again Xa can be simplified under the assumption of universality, so the 
value is independent of the identity of the meson or the weak decay vertex, 
and is parameterised by 

X A = (l+p A e^)\n^ (2.32) 



2.3.3 Isospin decompositions for B — > tttt 

For our study of the non-factorisable corrections to charmless B-decays in 
Chapter 4, we decompose the B — > nn decay amplitudes in terms of their 
isospin amplitudes. The factorisable contributions to these are calculated 
within the QCD factorisation framework, so it is useful at this point to discuss 
the decomposition of the decay amplitudes into their isospin components. 

In B — > 7T7T decays, the final state (nn) can only have total isospin of I = or 
1 = 2. This implies that the three decay amplitudes, B° — > 7r + 7r~, B — > 7r°7r° 
and B + — > 7r + 7r°, have only two decay paths, with change in total isospin 
between the initial and final state of AI = 1/2 or AI = 3/2. They must 
therefore obey a triangle relation [49] 

V2A(B + -> tt+tt ) - A(B° -> 7T+7T-) - y/2 A(B° -> = (2.33) 

and similarly for the CP-conjugate decays. Assuming isospin invariance in 
the matrix elements A(B — > nn) = (7T7r|7i e //|S), we can decompose the 
decay amplitudes into two isospin amplitudes A3/2 and A1/2 (where the index 
corresponds to the change in isospin). Explicitly we have: 

V2A(B + -> tt+tt ) = 3A 3/2 

^(S°^7r + 7r-) = A 3/2 -A 1/2 
V2A(B° -> ttV) = 2A 3/2 + A 1/2 (2.34) 

with equivalent expressions in terms of Aj for the charge-conjugate decays. 
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We can break the isospin amplitudes down further into contributions from 
two weak amplitudes, 



m = KA u t + a:a c 

A = K^i + \ C A\ (2.35) 

with \ q = V q bV* d . The phase difference between \ u and A c is related to the 
angle 7 of the unitarity triangle: 

hL = -^P^f e ^ (2.36) 



As was first reported in [50], it is straightforward to use the explicit expres- 
sions for the decay amplitudes (B.l) to obtain the following expressions for 
the factorisable part of the isospin amplitudes: 



AU 

/1 l/2 



AU 

^3/2 



^3/2 



A ~ 



6 



A c 



2 

A, 



[4ai - 2a 2 + 6« + rial) + 3(a 7 - a 9 ) + 3« + r£a£)] 
[2(a c A + r£<) + (a 7 - a 9 ) + K + rjajj)] 
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[2ai + 2a 2 - 3(a 7 - a 9 ) + 3« + r*a%)] 
[_( a? _ a9 ) + ( a « Q + r £ a «)] 



(2.37) 



with 



4nr = i^ml^CO)/, 



2ml 



r = 

x 



2m\ ) m q 



(2.38) 



2.4 Limitations to QCD factorisation 

The predominant problem with the factorisation framework is that it is only 
explicitly valid in the heavy quark limit - the master formula (2.4) is exact 
only as mj — > 00. However, as we know is ~ 5 GeV it is justifiable to 
examine the validity of this limit. The power corrections in the Aqcd/^6 
expansion are often difficult to calculate, and include non-factorisable cor- 
rections that cannot yet be determined in a model-independent manner. 

The expected size of the non-factorisable corrections is C(Aq C £>/to 6 ) ~ 10%, 
yet as we will show in detail in Chapter 4, the power corrections can be and 
are often required to be, considerably larger than this estimate. 

The power-suppressed terms can be large with respect to the leading order 
contribution thanks to various different sources of suppression - such as small 
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values of Wilson coefficients, CKM factors or colour suppression. Numerical 
enhancement of the power-suppressed terms themselves can also occur, where 
the coefficients can be larger than suggested by power counting, e.g., the 
chirally enhanced terms in coefficients aQ^nn) or as(7rn) discussed in Sections 
2.3.1 and 2.3.2. 

Comparing the predictions from QCD factorisation to the available exper- 
imental data gives very good agreement in many cases. For charmless B 
decays however, the correspondence between prediction and data gives what 
has been called "The B — > nn and B — > ixK puzzles" . We can consider two 
independent ratios of CP-averaged branching fractions: 



BR(i? ± -> 71" 71" ) 
BR(5° -> 7T + 7T-) 



T~B" 



T3TX-K 

-"-00 



BR(-B° -> 7T 7T ) 
BR( J B° -> 7T + 7T-) 



The branching ratio for B — > 7r + 7r~ is unexpectedly small compared to the 
prediction, and 5 — > 7r°7r°, unexpectedly large. 5 ± — > 7r ± 7r is however in 
agreement with the theory estimates. The experimental values, using the 
HFAG averages, are K%L ~ 2.2, K$ ~ 0.32. Using the default parame- 
ters the QCD factorisation estimates are somewhat different: FF+_ ~ 1.25, 
Rqq ~ 0.03. Even using the "favoured scenario" for the parameters and with 
enhanced annihilation contributions we find: R™_ ~ 1.8, -Rqq ~ 0.13 [3]. An 
analogous situation exists for B — > nK where moderate discrepancies occur 
in some branching ratios, such as the ratio of B — > n T K ± and B ± — > ^K . 

There is a reasonable agreement for the B — > 7T7t CP asymmetry predictions 
from QCD factorisation. There is however difficulty in making conclusive 
statements regarding this due to the discrepancies between the experimental 
measurements from BaBar and Belle. This can be seen graphically in the 
presentation of our results in Chapter 4, where we clarify and quantify the 
discrepancies between the predictions and measurements. 
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Chapter 3 

Light-cone meson distribution 
amplitudes 



This chapter is devoted to the study of distribution amplitudes of light 
mesons, where we derive new models characterised by a small number of 
parameters directly related to experimental observables - this is the first 
new result of this thesis. We discuss how these new models are constructed 
and how they are more descriptive than the current expressions of light me- 
son distribution amplitudes. We also show how we can numerically evolve 
the new model distribution amplitudes in the scale //, to exactly reproduce 
the leading order scaling behaviour of the usual approximated DA. Finally 
within the framework of QCD factorisation, we study the effect of these 
model DAs on important 2-body non-leptonic B-decays such as B — > nn. 

3.1 General framework 

The basic equation to describe bound states in relativistic quantum field 
theory is known as the Bethe-Salpeter equation, which was first formulated 
in the 1950's [51,52]. The leading twist distribution amplitude of a meson 
is related to this wave function by integrating out the dependence on the 
transverse momentum p±, 



A light from the shadows shall spring 



J.R.R. Tolkien 




(3.1) 



The light-cone distribution amplitudes were originally introduced in the con- 
text of hadronic electromagnetic form factors and the pseudoscalar-photon 
transition form factor [53-57], but have attracted interest in B physics due 
to their appearance in QCD sum rules on the light-cone [46,58-66], form 
factors [67,68], and their use in factorisation for B decay amplitudes [1,2]. 

Using the formalism introduced in Section 2.2.2, we can define the distribu- 
tion amplitude for a general pseudoscalar meson P in terms of the matrix 
element of a non-local operator near the light-cone. We have 

r 1 T 1 

(0\q 1 (x)^ 5 [x : -x]q 2 (-x)\P) = tf P pJ due^ x <j> P {u) + - m 2 P x 2 A P (u) + 0(x 4 ) 

Jo L 4 

+ \ fpmp — xj du e^ x [M P (u) + 0(x 2 )] 
I px J 



(3.2) 



where £ = 2u — 1 and [x,y] is the Wilson line defined in (2.7). In equa- 
tion (3.2), (j) P is the leading twist-2 distribution amplitude, whereas A P and 
Bp contain contributions from higher-twist operators. The corresponding 
definitions of vector meson DAs can be found in [69,70]. 

The theory of meson DAs is well understood [41,69-72], and suggests their 
parameterisation in terms of a partial wave expansion in conformal spin (see 
Section 3.2), allowed due to the conformal invariance of QCD. This expansion 
is in terms of contributions with different conformal spin, which do not mix 
with each other under a change of scale. This is true to leading logarithmic 
accuracy, but is no longer the case at higher order, as the underlying sym- 
metry is anomalous. The leading- twist distribution amplitude <f>(u), for both 
pseudoscalar and vector mesons has a conformal expansion given in terms of 
G eg enbauer polynomials Cf/ 2 , 

oo 

<P(u, /i 2 ) = 6«(1 - u) a n (fi 2 )C 3 J 2 (2u - 1) (3.3) 

?1=0 

The Gegenbauer polynomials are generalisations of the associated Legendre 
polynomials; the first few are 

Cf(x) = 1 
Cf /2 {x) = 3x 

Cl /2 {x) = 3 -(-l + 5x 2 ) 
The coefficients a n in (3.3) are Gegenbauer moments, and renormalise mul- 
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tiplicatively to leading logarithmic accuracy 



a n (Q 2 ) = a n (^) (^y) (3-4) 

where @ = 11 — (2/3)n/. For pseudoscalar and longitudinally polarised 
vector mesons, the one- loop anomalous dimension ^\ given by [73] 

^ = 8C '(*(" + 2 ) + ^-!- 2(» + l)(» + 2) ) (3 - 5) 
and for transversely polarised vector mesons 

fn) - ~ ( . - ~ 3 



7^ = 8C F ^P(n + 2)+ lE --j (3.6) 

For 7r, p, a;, 77, 77' and 0, G-parity ensures that a 0( jd = and that the DA 
is symmetric under u <-> 1 — u, whereas for X and fT* the nonzero values 
of a Q dd induce an antisymmetric component of the DA. ao = 1 is fixed by 
normalisation 



-1 

.2\ 
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and all the other a n are intrinsically nonperturbative quantities. As they do 
not mix under renormalisation to leading order accuracy, equation (3.3) is 
well suited to construct models for <fi: truncating the series after the first few 
terms yields a parameterisation of the DA that is "stable" under a change 
of scale, except for the numerical values of a n . Despite there being no small 
expansion parameter, such a truncated conformal expansion is often a mean- 
ingful approximation to the full distribution amplitude. An advantage of 
this expansion is that the contributions of higher conformal spins to the 
convolution integrals involving the LCDA (and therefore also the physical 
amplitudes) are suppressed by the highly oscillating behaviour of the par- 
tial waves. This suggests a construction of models for the DAs based on a 
truncated conformal expansion, where only the first few waves are included. 

We now move on to discuss the construction of the distribution amplitude 
via the conformal expansion using the conformal invariance of QCD. 



3.2 Conformal symmetry 
3.2.1 The conformal group 

The conformal group is defined as the group of general co-ordinate transfor- 
mations of 4-D Minkowski space that conserve the interval ds 2 = g iiv dx il dx v 
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up to a change of scale - i.e. transformations that preserve angles and leave 
the light-cone invariant. Transformations of this type are: scale transfor- 
mations (x^ — > Arc M ), inversions (x^ — > x^/x 2 ), translations (x^ x^ + c) 
and Lorentz rotations. The conformal group is the maximal extension of the 
Poincare group that leaves the light cone invariant [74]. In four dimensions 
the full conformal group has 15 generators, which are denoted by 



P M 4 translations 

M M!y 6 Lorentz rotations 

D dilatation 

K M 4 special conformal transformations 



An important subgroup of the full conformal group is the collinear subgroup, 
which (as we will see later) is the most relevant in the study of QCD. Denoted 
SL(2,~R), this subgroup is made up of a special case of the special conformal 
transformation 

... x^ + a^x 2 , „. 

x* _> = 3.7 

1 + 2a ■ x + a 2 x 2 v ; 

for a M light-like; this can be reduced to 

x_ -> x' = — (3.8) 

" l + 2ax_ V ; 

which maps the light-ray onto itself in the x_ direction. Transformations of 
this type together with the dilations and translations along the x_ direction 
make up the collinear subgroup. 

Transformations of SL(2,M) are governed by four generators P+, M_ + , D 
and K_. These are most often presented in the linear combinations [75]: 

L + = Li + iL 2 = -iP+ 
L = L x - iL 2 = iK_ 

L = f(D + M_ + ) 

E = i(D-M_+) (3.9) 

so that the algebra of SL(2, M) is written 

[L , L T ] = =fL t [L_, L+] = -2L (3.10) 

These generators act on fundamental fields $(x) which are equivalent to 
primary fields in conformal field theory. If we consider the parton model of 
hadron states, the hadron is replaced by partons moving collinearly in some 
direction n^. We can then consider only the fields confined to the light-cone 

— > $(an) 
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with a a real number. Additionally, we choose the field $ as an eigenstate of 
the spin operator so that is has a fixed spin projection s in +direction. The 
action of the generators can be described in terms of differential operators 
acting on the primary field. <p(a) = <f>(an) 

[L+,$(a)] = -d«$(a) = L+$(a), 
[L_,$(a)] = (a 2 d a + 2ja)<f>(a) = L_$(a), 
[L ,$(a)] = {ad a +j)$(a) = L $(a) 

[E,$(a)] = ±(Z -*)<&(«) (3.11) 

The generator E measures the collinear twist of the field, defined as the 
dimension minus spin projection in + direction. The field &(ax), with a fixed 
spin projection on the lightcone (x 2 = 0) is an eigenstate of the quadratic 
Casimir operator L 2 , defined by 

L 2 = L 2 + L 2 + L 2 . = + ^+ [L 2 ,L 4 ]=0 (3.12) 

so that 

L 2 $(«) = J2 [ L i, [ L i, *(«)]] = JU ~ W«) (3-13) 

i 

Here j = |(Z + s) and is referred to as the conformal spin, where I is the 
canonical mass dimension of the field and s the (Lorentz) spin. Consequently, 
we see that the field $(«) is transformed according to representations of the 
collinear conformal group SL(2,R), specified by the conformal spin j. 

3.2.2 Conformal symmetry in QCD 

QCD is conformally invariant at leading order, but is broken at next-to- 
leading-order by the inclusion of quantum corrections. However, the confor- 
mal spin is seen as a good quantum number in hard-processes up to small 
corrections of order a 2 . The non-local operators of QCD can be written 
in terms of conformal operators. We can construct a complete basis of lo- 
cal operators on SX(2,R) from $(0) by applying the "raising" operator L + 
/c-times. 

O = $(0) 

O k = [L + ,...,[L + ,[L + ,$(0)]]] = (-9 + ) fc $(a)| Q=0 
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The primary fields $(a) can be formally expressed as a Taylor series expan- 
sion of local conformal operators 



OO / \ JU 



fc=0 



(3.14) 



This is a simple example of a conformal tower. It is possible to construct 
conformal towers for general operators built of a number of primary fields and 
derivatives [75] . We can define a local conformal operator O n that transforms 
under the collinear conformal subgroup. This is equivalent to demanding the 
operator satisfies 

[L 2 ,O n ] = j(j -l)O n 
[L ,O n ] = j O n 
[L_,O n ] = 

The tower of operators can be built by repeatedly applying the raising oper- 
ator, as above 



''n,n+k 



[L + , . . . , [L + , [L. 



-d + y 



It is possible to construct explicitly a general local conformal operator as [74] 



Q>^ n (x) 



-1) 



d + + d 



(3.15) 



where P% ,b {x) are Jacobi polynomials and the superscripts indicate the con- 
formal spins of the constituent fields. 

Returning to QCD, we can extract different conformal operators related to 
the different spin and dimension of the constituent quark or gluon fields. For 
a quark ~ anti-quark operator at the leading twist we have the relevant local 
conformal operator O^' 1 , given by 



(x) = (id + ) n \i,{x) 1+ CT (*D + /d + ) ij(x) 



(3.16) 



where the Gegenbauer polynomials are a special case of the Jacobi polyno- 



mials py = C\ 



3/2 
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3.2.3 Conformal partial wave expansion 



The aim of the partial wave expansion is to make full use of the underlying 
symmetry of the theory in order to simplify the dynamics of the problem. In 
the case of QCD, we make use of the conformal invariance to create a partial 
wave expansion in conformal spin. The conformal expansion of the LCDA is 
analogous to the partial wave expansion for a spherically symmetric potential 
in quantum mechanics. Here, the rotational symmetry allows the separation 
of angular and radial degrees of freedom. The angular dependence is encoded 
in the spherical harmonics Y l m (0, <p) which form an irreducible representation 
of 0(3), and the radial information is governed by an evolution equation, 
namely the 1-d Schrodinger equation in R(r). 

In QCD, the partial-wave expansion decomposes the distribution amplitude 
into longitudinal and transverse degrees of freedom. The dependence on 
the longitudinal momentum is expressed in terms of orthogonal polynomials 
which form irreducible representations of SL(2,M). For a general two-body 
wavefunction, these polynomials are the Jacobi polynomials. For conformal 
spins ji = J2, we obtain the Gegenbauer polynomials, defined by the gener- 
ating function 

(1 _ 2xf + t .y =EW 

v ' n=0 

A multiple-particle state that is built of primary fields can be expanded in 
terms of irreducible representations of SL(2, M.) with increasing conformal 
spin. The lowest possible conformal spin for the multi-particle state equals 
the sum of spins of the constituents, and its wavefunction is simply a product 
of non-degenerate one particle states. This lowest state defines what is known 
as the asymptotic distribution amplitude, as it is reached at the formal limit 
of q 2 — > oo. This state is expressed, in the general 

*„(«!,«,...«„,) = ^^'''r^j ^ 1 ' 1 ^' 1 ■ • • M - m " 1 (3 - 18) 

where Ui is the momentum fraction of the i th constituent, and the normali- 
sation is chosen so that f[dui](f) as (ui) = 1. 



Light-cone meson distribution amplitudes 

The distribution amplitude for the pion <f> n (u) is an expansion over an infinite 
series of Gegenbauer polynomials, with multiplicatively renormalisable co- 
efficients 

oo 

faiu, fi 2 ) = 6«(1 - «) *n(v 2 )C 3 J 2 (2u - 1) (3.19) 

n=0 
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The Gegenbauer polynomials form an orthonormal set on the interval < 
u < 1 with weight function u(l — u) and satisfy the orthogonality condition 

f duu{\ - u) Cf/ 2 (2u - \)C% 2 {2u - 1) = 5 mn N n (3.20) 
Jo 

with 

_ (n+l)(n + 2) 
7Vn " 4(2n + 3) 

The renormalisation group equation for <p w (u, fi 2 ) is known as the Efremov- 
Radyushkin-Brodsky- Lepage (ER-BL) evolution equation [55,57] 

d f 1 

fi 2 -^(frn(u,fi 2 ) = J dvV(u,v,a s )(p w (v,fi 2 ) (3.22) 

The evolution kernel V(u, v,a s ) has been calculated perturbatively to next- 
to-leading order in a s [76,77] 

rv (y^ 

V(u, v, a s ) = ^V (u, v) + -^V^u, v) (3.23) 

The pion distribution amplitude can be expressed as matrix elements of 
renormalised local operators 

(0|J(0) 7+ 75(^ + )^(0)|7T + (p))=2/ 7r (p + )" +1 f du(2u-irMu,fi 2 ) 

Jo 

(3.24) 

so the relevant conformal operators (3.16) are then [74] 

O n (x) = (id + ) n [d(x) 1+l5 C 3 J 2 (jD + /d + ) u(x)] (3.25) 

Due to the flavour structure of this operator, there is no mixing with opera- 
tors involving two-gluon fields, nor with operators of three or more fields, as 
these appear at higher twist. Thus the operators of (3.25) do not mix at the 
leading order, and are multiplicatively renormalised. The evolution equation 
(3.22) can then be solved by a conformal spin expansion 

Mu,fi 2 ) = ^p^-C 3 J 2 (2u - l)((O n (fi 2 ))) (3.26) 



N, 

n=0 1 



where 



((O n (fi 2 ))) = f duC 3 J 2 (2u-l)Mu,V 2 ) (3.27) 
Jo 

The sum runs over even values of n to ensure symmetry of the pion DA, 
due to charge conjugation invariance: <f> n (u) = <f> n (l — u). We can construct 
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the Gegenbauer moments, a„(/i 2 ), using the reduced matrix elements of the 
conformal operators given in (3.27), such that 



a„(/i 2 ) = 



2(2n + 3) 
3(n + l)(n + 2) 

2(2n + 3) 
3(ra + l)(ra + 2) 



f duC^ 2 {2u-l)<t>^u^ 2 ) 
Jo 



(3.28) 



3.3 Non perturbative input 

In general, not much is known about the amplitudes of these partial waves 
and hence the numerical values of the moments a n . There is some informa- 
tion available on the lowest moments of the 7r DA, but much less for the 
other mesons. For the ir and r] mesons, there is some experimental informa- 
tion available on the ^(77)77* transition form factor [78], which can be used 
to extract the values of the first two amplitudes. The 77* — > re form factor 
factorises into a partonic hard scattering amplitude and a soft hadronic ma- 
trix element parameterised by the pion DA <j) v . The leading twist expression 
for F ni is given by [79,80]: 

4 f 1 

F^{q 2 ) = -j=J dx <f> v {x, q 2 ) T^ n (x, q 2 ) 

^ vfef*^ (1 + OK,) (3 ' 29) 

using <f> n (x) = <f> n (l—x) and since the amplitude for the subprocess 77* — > n 
(labelled T^J is 

This means that the 7 — n transition form factor approximately probes the 
sum of the Gegenbauer moments via 

f 1 dx ^M = i + y an (3.30) 

By making some assumption (such as truncating this series after n — 4) it is 
possible to extract information about the moments. 

On the theory side, there exist a few lattice calculations for the second mo- 
ment of the 7r DA; some quite dated [81-83], and a recent retry [84,85]. 
Unfortunately, these results are still preliminary and cannot yet be used 
in phenomenological applications. Other theoretical calculations have been 
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done using QCD sum rules for both pseudoscalar and vector mesons; these 
are reviewed below. 



3.4 Previous constructions of ^ 

This section briefly reviews and consolidates the previous calculations (and 
constructions) of the pion DA and its moments a n . 

The earliest calculation of one of the moments of the pion DA was that of 
a 2 by Chernyak and Zhitnitsky (CZ) [41]. They developed a model of the 
pion DA based on local QCD sum rules whereby (p n is expanded in terms 
of local operators from which the a n can be extracted. The form of the 
DA is strongly peaked at the endpoints, due to the approximation of the 
vacuum quark distribution by a delta function and derivatives of, the non- 
perturbative contributions to the sum rule forcing the DA to be end-point 
concentrated. The DA is shown graphically later in this section in Figure 
3.1, and is described by a model wavefunction of 

<j)C z (u) =30u(l-u)(l-2u) 2 (3.31) 

which corresponds to the numerical result 

a 2 (0.5GeV) = 2/3 (3.32) 

The drawback of extracting the a n in this way is that the expansion of an 
intrinsically non-local quantity, namely the pion DA, in terms of contribu- 
tions from the local operators causes an increased sensitivity to the non- 
perturbative effects. The coefficients of the condensates in the sum rule 
increase with n and can in fact dominate over the perturbative contribu- 
tions for a large enough value of n. This implies that this method should be 
reliable for at least n = 2 but not for large n. They also obtained the result 

0(1/2, 1/2 GeV) = (3.33) 

neglecting all a n > 4 . This however is an artifact of the neglect of higher order 
moments. 

This method was improved by Braun and Filyanov [29], who found that 
0(1/2) was required to be non-zero; they found a constraint on 0(1/2) at the 
scale of 1 GeV of 

3 9 45 
(1/2, 1 GeV) = 1.2 ± 0.3 = - - -a 2 (l) + — a 4 (l) (3.34) 
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They also re-determined the first two moments using the CZ procedure: 



a 2 (lGeV) = 0.44 a 4 (lGeV) = 0.25 

Other calculational methods have also been used to determine the moments of 
the DA. For example, the most recent calculation is from Ball and Zwicky [46, 
58] using light-cone sum rules, who have determined the first two moments 

as 

a 2 (lGeV) = 0.115 a 4 (lGeV) = -0.015 

Another alternative method is to use sum rules with non-local condensates, 
as developed in [86-90]. This method constructs the DA by connecting the 
dynamic properties of the pion with the QCD vacuum structure. This allows 
for the fact that quarks and gluons can travel through the QCD vacuum 
with non-zero momentum k q , which implies a non-zero average virtuality, 
(k^) = \g, using the ratio 

K = {qa f 9) (3-35) 



The most recent calculations using this technique was by Bakulev, Mikhailov 
and Stefanis (BMS) [87], who used a virtuality of \ 2 q = (0.4 ± 0.1)GeV 2 
and extracted the values of a 2 to a w from the non-local condensate sum 
rules. The non-local quark condensate represents a partial resummation 
of the operator product expansion to all orders, for the vacuum expectation 
value of the (non-local) operator {q(0)E[0, z]q(z)). This is expressed in terms 
of analytic functions Fsy, whose derivatives are related to condensates of the 
corresponding dimension 



F s (z-) - 1 ±F V (*) 



(3.36) 



E[0, z] is an appropriate Wilson line operator, taken in an appropriate gauge 
so that the path-ordered exponential equals unity. The lowest order conden- 
sates are found via standard QCD sum rules [91], and are 

Q 3 = (qq) 4* = ig.(qG' a 'w) Q % = <(??) 2 > ( 3 -37) 

The functions F$y are given in terms of a virtuality distribution which de- 
scribes the distribution of the vacuum fields in terms of \ 2 q 

f oo 



where 



/•oo 

F s ,v(z 2 ) = / e az2 / A f sy (a)da (3.38) 
Jo 

f°° f 1 S 

I fs ^ a) da= \0 V and chiral limit (3 ' 39) 
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In the BMS approach, this distribution is modelled by performing an expan- 
sion where the large average virtuality (compared to the relevant hadronic 
scale) is included in only the first term of the series, instead of the usual 
expansion of the non-local condensate in terms of local condensates (qq), 
(qD 2 q) .... This corresponds to taking the form of the virtuality distribu- 
tion as f s (a) = 5(a - X 2 q /2) and f v ~ a s Q 3 5'(a - X 2 q /2) [90]. 

The pion distribution amplitude can then be connected directly to the non- 
local condensates by means of a sum rule 

i=S,V,Tj 

The s and M are Borel parameters and the are local condensates, with 
i running over all scalar, vector and tensor types. Using this sum rule pro- 
vides values for the first ten Gegenbauer moments of the pion distribution 
amplitude; for n > 4, the a n were found to be ~ 10~ 3 and so were neglected 
in their calculations. The non-zero moments were found to be 

a 2 (1.16GeV) = 0.19 a 4 (1.16GeV) = -0.13 (3.41) 



Finally, the extraction of the moments a-i and a 4 from the photon-pion tran- 
sition form factor (as discussed above), has been developed and refined by a 
succession of authors [92-95], culminating in the constraint 

a 2 (lGeV) + a 4 (lGeV) = 0.1 ±0.1 (3.42) 



The moments of the DA which are obtained in all of these models are sum- 
marised in the table below, where the moments are all normalised to a scale 
of q 2 = 1 GeV 2 



Model 


a 2 


a 4 


AS 








CZ 


0.56 





BF 


0.44 


0.25 


BMS 


0.20 


-0.14 


BZ 


0.115 


-0.015 



The models of the DAs are shown for comparison in Figure 3.1. 
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Figure 3.1: Model pion distribution amplitudes for comparison between asymp- 
totic (dotted), BZ (blue), BMS (red) and CZ (green), all taken as 
scale q 2 = 1 GeV 2 . 



Here we discuss the main new result of this chapter, where we introduce new 
models for the leading-twist distribution amplitudes of the it and the K based 
on the fall-off behaviour of the nth Gegenbauer moment of these amplitudes, 
a n in n, which is assumed to be power-like [5]. The models are formulated 
in terms of a few parameters, notably the first inverse moment of the DA 
(which for the 7r is directly related to experimental data) and the strength 
of the fall-off of the a n , and can be summed to all orders in the conformal 
expansion, parameterising the full DA at a certain energy scale. Although 
we focus on the pseudoscalar mesons n and K, the models we propose are 
equally well applicable to vector mesons p, u>, K* and <f>. 

3.5.1 Motivation 

For many processes involving distribution amplitudes it is usually argued 
that a truncated conformal expansion would be sufficient for the calculation 
of physical amplitudes so long as the perturbative scattering amplitude is 
"smooth" — the reason being the highly oscillatory behaviour of higher order 
Gegenbauer polynomials. In order to quantify this statement, consider the 
simplest case of one meson in the initial or final state, so that the convolution 
integral reads 



3.5 New models for LCDA 




(3.43) 



where T is the perturbative scattering amplitude. 
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This can be separated into three cases: 



(i) T is nonsingular for u G [0,1], 

(ii) T has an integrable singularity at one of the endpoints, 

(iii) T contains a non-integrable singularity at one of the endpoints. 

As a typical example for case (i), consider T(u) = y/u, which yields 

f\ a ( \ r V- (-l)"36(n+l)(n + 2) 

Jo dU<P[U)Vu ~^ (2n- l)(2n + l)(2n + 3)(2n + 5)(2n + 7) ° n 

This result implies a strong fall-off ~ 1/n 3 of the coefficients of higher Gegen- 
bauer moments a n : assuming a« = 1 for all i, already the first three terms in 
the sum account for 98.8% of the full amplitude. In reality the convergence 
will be even better as all existing evidence points at \a n \ <C 1 for n > 1. 

As an example for case (ii) consider T(u) = lnw, giving 

The singularity at u — evidently worsens the convergence of the series; 
again assuming Oj = 1, the first three terms now overshoot the true result by 
35%. In order to approximate the full amplitude to within 5% one now has 
to include nine terms, but the convergence will again be better in practise, 
thanks to the fall-off of a n with n. 

Case (iii) is more complicated and depends on the asymptotic behaviour of 
the a n . For T(u) — \/u, for instance, one obtains 



oo 

d M 0( M )- = 3^(-l) n a n (3.44) 



Here the amplitude is finite only if the a n fall off sufficiently fast in n. For 
stronger endpoint divergences the coefficients multiplying the a n start to 
grow in n and for T ~ 1/u 2 the integral diverges, even for the asymptotic 
DA, which would indicate a breakdown of factorisation for that process. 

This discussion suggests that models of <p, based on a conformal expansion 
truncated after the first few terms, can be used appropriately for cases (i) 
and (ii), but are less reliable for case (iii). Convolutions with T ~ 1/u are 
very relevant both in hard perturbative QCD e.g. 77* — > n (Section 3.3), and 
in decays such B — > tytt. The a n can have different weighting in the different 
convoluted amplitudes, and have the highest impact for the convolutions of 
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type (iii). Instead of basing a parameterisation of <fi on a<i and (taking all 
On>4 = 0), it is not unreasonable to use the case in (3.44), where all a n enter 
with the highest possible weight factor. This is the basic idea behind our 
models of leading-twist DAs. 



3.5.2 Model construction 



In beginning to construct a model of the (leading-twist) pion DA we consider 
the two properties that any viable model must fulfil. Firstly, the large-g 2 
behaviour must be satisfied for all light-meson DAs, in that the distribution 
amplitude must tend to the asymptotic form <f) n (u,fi 2 — > 00) = 6u(l — u). 
Also, in order for QCD factorisation in B decays to make sense the first 
inverse moment, j Q du(f) n (u)/u, which as we have seen is related to the 7177* 
transition form factor, must exist. Both these conditions are fulfilled for 
models based on a truncated conformal expansion, which has in addition the 
prediction that for u — > 0, 1, (f) ~ u(l — u), independently of the factorisation 
scale. This prediction is in general, not fulfilled for models that are not 
truncated at fixed order in the conformal expansion, as shown below. 

As was introduced in equation (3.30), the first inverse moment of the pion 
DA <p n is related to the sum of all the Gegenbauer moments. We define this 
as our first "important parameter" , A, via: 

dU 3^ = AG " 2) = 1 + D" 1 )™ ^ 2 ) ( 3 - 45 ) 

n=l 



If this sum is to be convergent then the d2 n must fall-off in n sufficiently fast 
(recalling that for the ir, all odd-numbered a n vanish due to G-parity). If 
we assume an asymptotic equal-sign behaviour for the moments with large 
n, then the slowest possible fall-off is power-like a 2n ~ l/n p , with p slightly 
larger than 1. When the distribution amplitude is defined with a power-like 
fall-off it is possible to explicitly sum over the Gegenbauer moments. This is 
done using the generating function for the Gegenbauer polynomials: 

(it^t?^ = 5 c " /2K)fl (3 - 46) 



To show this, we begin from the full conformal expansion for the pion distri- 
bution amplitude, suppressing notation for the /^-dependence for compact- 
ness and defining u = (1 — u). We define the moments with a power-like 
fall-off as 

an = WbTW (3 ' 47) 
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for n even. This gives us an expression for the DA (for even n) of 

oo 

<j,+ (u) = Quu C n /2 ( 2 u - 1) a n 

n=0 

oo 

= Quu J2 C n /2 (2u-l) 



n=0 



(n/b+ l) a 



where + denotes the same-sign behaviour of the Gegenbauer moments. We 
can quite easily extend this across all values of n, giving 



oo , 

<P + (u) = 3uu C"J 2 {2u - 1) (1 + (-1)") 

n=0 ^ 



b + n 



We can express the (a, b) dependent part of this equation in terms of an 
analytic integral equation 

f _ 5— V = -J_ f 1 dti-kity- 1 ^" 

\b + nj T(a)J V ; 



giving 



<f> + (u) 



3uu 



r(o) Jo 



dt(-\nt) a ^ ^C n 3/2 (2 M -1) (l + (-l) n ) t n/b 

n=0 

The sum can be split into two parts 

oo oo 

cli\2u - 1) t n ' h + Cn 3/2 (2« - 1) HT /b 



n=0 



n=0 



which can be directly related to the generating function given in equation 
(3.46), and gives a final expression for the DA of 



3uu 



J 1 dt(- In ty- 1 (j (2u - 1, t l ' b ) + f (2« - 1, -t 1/6 ) ) 

(3.48) 



Using an analogous method we can obtain a re-summed DA with alternating- 
sign behaviour of the Gegenbauer moments, 



ra/2 



(-1) 
(n/6 + l) a 



for n even 



of 



= fjfi [ *t(- lnO"" 1 (/ (2« - 1, it 1 /') + / (2u - 1, ) 

(3.49) 

From equation (3.45) we can find the corresponding values of A 

= ( & /2) a C(«, 6/2) A- b = (6/4) a (C(a, 6/4) - C(«, 1/2 + 6/4)) (3.50) 

where ((a, s) = X)*~=o + s ) a * s a generalisation of the Riemann zeta 
function known as the Hurwitz zeta function. In order to obtain models for 
arbitrary values of A we split off the asymptotic DA and write 

0a, 6 (A) = 6u«+^^(^(u)-6uti) (3.51) 

a, ft 

valid for a > 1 and 6 > 0. This equation implies that the asymptotic DA is 
recovered for A = 1 and also from <p^ b in the limit a — > 1. Examples of these 
DAs for various values of a, 6 are shown in Figures 3.2 and 3.3 in the next 
section. 

The fall off of the a n in inverse powers of n gives a compact, closed expression 
for the DA using equation (3.51). This sort of fall-off behaviour is in fact in- 
trinsic to QCD, which we can see if we consider the behaviour of the moments 
a n (fx 2 ) under a change of scale. Taking the models defined at the hadronic 
scale /i ~ 1.2GeV, we know the a n scale with jj, according to equation (3.4). 
For large n we have 

7 J n) "~° O 8C F lnn + 0(l) 

and 

a n(Q 2 ) ~ n 4C F /A, ln(l/L) M/^) 

with L = a s (Q 2 )/a s (fi 2 ). This shows that the leading-order scaling induces 
a power-like fall-off of the a n , at least for large n. Another consequence of 
this is that as Q 2 — > oo (L — > 0) the suppression of the higher order moments 
is power-like so that 

^(Q 2 -> oo) = 4^ (a -> oo) = 6u(l - u) 

Hence the DAs (p^ b defined in (3.51) approach the asymptotic DA in this 
limit, so both of the necessary conditions for the DA construction are satis- 
fied. 
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3.5.3 Properties of model DAs 



Once the parameters A, a and b are specified at a certain reference scale, 
the model DA can be evolved to a different scale via the evolution equation 
(3.22), using the method which we describe in Section 3.6 - which shows 
how the evolution equation can be solved numerically to leading logarithmic 
accuracy. The procedure for the full analytic evolution of a DA expressed as 
a conformal expansion is included for completeness in Appendix C. 

After evolution to higher scales, the DA is no longer described by the function 
0^ b (A) with a simple suitably chosen set of A, a and b. This lack of "form- 
invariance" under changes of scale is not however disadvantageous to our 
model. We determine the value of A at some low scale around \i = 1.2GeV 
(chosen according to the analysis of [87]) and ideally the other parameters 
a, b could then be fixed from some experimental or theoretical determination. 
There is an analogy to this elsewhere in QCD, for the parton distribution 
function, which has a parameterisation fixed at a low scale (typically ~ 2 
GeV), but has parameters which are fitted to experimental data obtained at 
a large variety of scales [96-98]. 

Figures 3.2 and 3.3 show that the two models + and 0~ have quite a dis- 
similar functional dependence on the momentum fraction u, notably that (\T 
becomes non-analytic at u — 1/2 for values of a < 3. The spike at u — 1/2 
characteristic for <p~ causes these models to significantly deviate from the 
asymptotic distribution amplitude even for A close to 1. 
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Figure 3.2: Left: Examples for model DAs (j)^ b as functions of u, for a = 
1.5,2,3,4, and constant b = 2 and A = 1.2 (solid curves). For 
a — > 1, (f)^ b approaches the asymptotic DA. Right: the same for (f>~ 2 . 
The asymptotic DA is also shown for comparison (dashed curve). 

What we can also discern from the form of (f> a> b is that it is much more sen- 
sitive to variation in the fall-off parameter a, than b. Since the experimental 
data available are too scarce to constrain all of the parameters of the model 
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Figure 3.3: Left: Examples for model DAs 4>^ b as functions of u, for b = 0.1, 1, 10, 
and constant a = 3 and A = 1.2 (solid curves). Right: the same 
for 0g b . The asymptotic DA is also shown for comparison (dashed 
curve) . 



we fix b = 2. Using the set of model DAs <f) a ,b we can find the possible 
values of the lowest non-zero Gegenbauer moments a 2 and 04. For a — 1, 
approaches the asymptotic DA with = 0; at a — > 00 the models 

approach the standard NLO fixed order conformal expansion with a 2 = A — 1 
and a 4 = 0. 

As discussed above, for models based on the truncated conformal expansion, 
the endpoint behaviour of the DA (as u — > 0, 1) is assumed to follow ~ 
u(l — u). This is no longer true for our new models for all possible values of 
a. For example, for the model + , we have 

a/m(1 — u) for a = 2 

u{\ — u) \au{\ — u) for a = 3 
u(l — u) for a > 3 

We see that the linear endpoint behaviour is only obtained for models with 
higher values of a. This is not however a drawback to the model. The original 
argument in the literature for this endpoint behaviour was put forward in [41] 
from the calculation of DA moments via QCD sum rules: 

Jo n 

This is a result of a leading order calculation of the perturbative contributions 
to the QCD sum rule. The NLO result [99, 100] does not have the same 
dependence, instead following 

0~u(l-u) In 2 (u/(l-u)) 
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0+ ~ 
(f) + ~ 



which is equivalent to a n ~ (1/n 3 ). The large-n behaviour of the non- 
perturbative terms cannot be determined from the sum rules, but there is no 
reason why it should not follow the behaviour of NLO perturbation theory 
or some other scaling. 

3.5.4 Constraints on model parameters 

The asymptotic distribution amplitude is recovered for all models with A = 1 
and with + in the limit of a — > 1. We must impose the constraint that 
a > 1, otherwise the models 0^ b will not vanish at the endpoints u — 0, 1. 
In order to constrain our parameters, we can now introduce the experimental 
restrictions on our parameter set A, a, b. As discussed above, we fix b = 2 and 
we also fix our reference scale as /i = 1.2 GeV. The available experimental 
data for the n, as summarised in [87], points to a value of A around 1.1 at 
the scale fx w 1.2 GeV. 

If we require a 2 to be positive, which as seen in Section 3.4 is a reasonable 
conclusion from all of the previous determinations, then we require A > 1. 
We can infer an upper bound from experimental data or by imposing the 
requirement that a 2 < 0.2 - the upper bound as found by the most recent 
analyses. This corresponding bound on A is A < 1.2 for + , and is smaller 
again for <f>~. The final constraint we must consider is the allowed range of 
7r (l/2) from light-cone sum rule analyses, which give [29]: 

0.9 < 7r (l/2, 1 GeV) < 1.5 

For A > 1, 0(1/2) is always smaller than 1.5, so only the lower bound is 
relevant. The upper bound on A* 1 for various values of a (as implied by the 
lower bound 0^(1.2) > 0.9) is 



a 


2 3 4 5 6 00 


A+ 

max 
^max 


2.04 1.58 1.43 1.36 1.33 1.27 
1.04 1.11 1.16 1.19 1.22 1.27 



The constraint on A + is weaker than those discussed above, but for 0~ the 
minimum value of 0^(1/2) poses a nontrivial constraint on A. Concentrating 
on the general characteristics of the distribution amplitudes in order to retain 
a generality for the other psuedoscalar and vector mesons, we do not refine 
these constraints further using specific data for the pion. It is therefore 
possible to draw up a set of general constraints on the (symmetric) part of 
our model DAs: 

• 1 < A < 1.2 with < a 2 < 0.2 for 0+ 2 : this is based on the observation 
[41] that DAs of mesons with higher mass tend to become narrower, 
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• 1 < A < Min(1.2, A max ) for 2 , with A max given in the table above, 

• b = 2, lacking further data. 

3.5.5 Extension to 4>k 

As mentioned before, we have to distinguish between those mesons for which 
the a Q dd vanish due to G-parity, such as the n, and the strange mesons, for 
which the odd moments are in general non-zero and induce an antisymmetric 
part into the DA. 

We can now consider the extension of our new model DAs to the leading- 
twist DA of the Kaon, (f>K- This differs from <j) n by the contribution of odd 
Gegenbauer moments. We can construct a model for this antisymmetric 
part in an analogous way as before, where we introduce A asym , c, d in place 
of A, a, b. Since there is almost no information at all about the form of the 
antisymmetric DAs we set d = 2 immediately. We have 

#(«) = ^ j\t{-\nt) c - 1 (/(2u-l,Vf)-/(2u-l,-Vf)) 

= ^ jT'cftC-ln*) 6 - 1 (f(2u-l,iVi)-f(2u-l,-iVi)) 

(3.52) 

The models give a value of the first odd moment as a\ = (2/3) c . If we wish 
to redefine these expressions to account for arbitrary a±, we can write 

= ai (3/2) c j>t(u) (3.53) 

An example of such a model is shown in Figure 3.4. The scale evolution 
proceeds as for the symmetric part. There is very little information about 
the di moment for the Kaon - not even the sign is known reliably [46]. From a 
heuristic standpoint it is expected a\ > 0. This is because the DA is expected 
to be skewed towards larger values of u, since u denotes the momentum of 
the s-quark (the heavier quark) in the meson. The first calculation using 
QCD sum rules was from Chernyak and Zhitnitsky [41], who confirmed this 
intuition with a central value of 

af (lGeV) = 0.17 (3.54) 

This result was countered by the work of [99], who claimed to find a sign 
error in the CZ method as well as adding higher order radiative corrections; 
this gave a value of a\ of 

af(lGeV) = -0.18 (3.55) 
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Figure 3.4: Models for the antisymmetric contributions to the twist-2 DA for 



It was reported in [46] that this re-analysis of the CZ sum rule calculation 
gave a problem with the q 2 dependence of the form factors. In light of these 
thoughts we use the original result from [41] and use the evolved value of 
ai = 0.15 at 1.2GeV. We can then combine the asymmetric and symmetric 
parts of the wavefunction to give a total DA for the K. We can write 



where A is the sum of the even-valued moments from the symmetric contri- 
butions to the DA, and we defined A 885 ™^ 



cti = 0.15. ip c is shown as function of u for c £ {1, 2, 3, 4}. Left: 
Right: ip~ . Like the symmetric models (j>~, ip~ is non-analytic at 
u = 1/2 for c < 3. 





I du 

o 



I du 

o 



i>t{u) 



ij c (u) 



ai (3/2) c C(c,3/2) 



oi(3/4) c [C(c,3/4)-C(c,5/4)] 



3.6 Numerical evolution of model DA 



Using the evolution equation for the meson DA 




(3.56) 



the kernel is given, to leading order in a s , as 



V(x,y,fi) = C F 



2tt 



V (x, y) + 0(a 2 s ) 



with, for pseudoscalar mesons, 



V (x,y) = V BL (x,y) - 5(x -y) dzV BL (z,y) 
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V BL (x,y) = L-E.(l + -LS)e(x-y) + -(l + —)e(y-x) 
1-y \ x-yj y \ y-xj 

Equation (3.56) is not well suited to the numerical evolution of 0, as it 
requires the use of an explicit formula for a s , which is itself the leading order 
solution of a renormalisation group equation. In order to exactly reproduce 
the correct LO scaling behaviour of the Gegenbauer moments numerically, 
we can rewrite (3.56) in terms of a differential equation in a s . To LO this is 
written as 

~7^r a s~^ — 4>( x > a s)= / dyV (x,y)(j)(y,a s ) (3.57) 
2C F oa s J 

which can be solved iteratively using Euler's method: 

4>{x, a s - Aa s ) = <j>(x, a s ) + — — / dyV (x,y)<j>(y,a s ) 

Po a s Jo 

We have checked that for Aa s = 0.01 and with 21 mesh points in x we 
correctly reproduce the known scaling behaviour of the truncated conformal 
expansion to within one per mille. We can show this for example, using the 
standard conformal expansion truncated at a 4 

4>(u, a s ) = Qu(l — u) 1 + a2{cts)C^ 2 {2u — 1) + a i {a s )C^ 2 {2u — 1) 

using the values for the first two moments from Ball and Zwicky [46] given at 
H = 1.2GeV =>> a L s ° = 0.6 as a 2 (0.6) = 0.115, a 4 (0.6) = -0.015. Figure 3.5 
shows the results for evolving this wavefunction up to a scale of 5.3 (af = 
0.26) using the numerical method, and for the full analytical evolution using 
the formulae in Appendix C. As we can see clearly, both methods give 
equivalent results. The change of our model DAs with the scale \i is shown 
for an example model (<j>f 3 ) in Figure 3.6. 
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Figure 3.5: Numeric (dashed red line) and analytic evolution (dotted blue line) 
of sample DA from 1.2GeV to 5.3GeV; un-evolved DA shown for 
comparison. 




Figure 3.6: Evolution of model DAs. Left: c/)J 3 (1.2) and Right: 0g 3 (1.2) from 
scale [i = 1.2GeV (solid curves) to (j, = 4.8GeV (dashed red line). 
We also show the asymptotic curve for comparison (dot-dash blue 
line). 



3.7 Application to B —> nn 

It is at this point that we can consider the effects of the leading-twist DA on 
the non-leptonic B decays calculated within QCD factorisation. As shown in 
Chapter 2, the matrix elements determining the B — > tttt decay amplitudes 
are written as convolution integrals involving the leading twist DA. We dis- 
cussed in Section 2.4 how the recent experimental data from the non-leptonic 
B decays (specifically B — > nix) imply the failure of QCD factorisation to 
explain the branching ratios and CP asymmetry measurements. This sec- 
tion investigates the effect of a non-standard DA on the predictions of QCD 
factorisation. 

In this study, we use the full QCD factorisation formulae for the factorisable 
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contributions (including the calculated l/m& corrections) and also include the 
model-dependent corrections for the hard-scattering terms (parameterised by 
the complex number Xh) which we include at the "default level" of Xh = 
\n(m n / 'Ah) ~ 2.4 [2]. We do not however include the non-factorisable weak 
annihilation terms. We use the input parameters as quoted in [3] except for 
that of the leading-twist <f> n , which we replace by our model parameter and 
the B — > 7i transition form factor which we take from light-cone sum rules 
as found in [46]. 

We begin by considering the time- dependent CP asymmetry in B — > 7r + 7r~, 
which is defined as 

A%> = ~ ^Tl ~ Z° - ^l] = sin (Ami) + CW cos (Am*) 

(3.58) 

where the mixing-induced and direct asymmetries (S wn and C n7V respectively) 
depend on the unitarity triangle angles (3 and 7 via 

2 Im\ n7T 1 — \X n7T \ . s 

™ ~ 1 1 ix — 12 ^ ~~ 1 , i\ — 12 

with 

x " = '""SI < 3 - 60 > 

In practise, we calculate the CP asymmetry via its parameterisation in 
Wolfenstein parameters, so that the CKM factors read 

±i7 _ p±ifj 2if3 _ (l-p) 2 -^ 2 -2^(l-p) 

~VFT? " F#T? (3 - 61) 

Although the penguin-to-tree ratio P/T is highly suppressed, it is not negli- 
gible and can be expressed via QCD factorisation in terms of CKM phases 
and pure strong interaction parameters. 

P re^ 



1 \ (aKmr) + ^(tttt)) + (a c w (inr) + r*a c 8 (inr)) 



v/p 2 + fj 2 J ai(7T7r) + (o^tttt) + rJa^TTTr)) + (a^(nn) + rJa^TTTr)) 

(3.62) 

P/T is then given by a ratio of polynomials in the Gegenbauer polynomials 
a n with complex coefficients. This enables us to write the time-dependent 
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CP asymmetry as 



s = 2y [p 2 + rf -r 2 - p(l - r 2 ) + (p 2 + - 1) r cos 0] g ^ 
[(1 — p 2 ) + ff 2 } [p 2 + f/ 2 + 2rp cos 0] 



In the limit where P/T is zero this reduces to S nn = sin 2a. 

The variation of the CP asymmetry in terms of A is shown in Figure 3.7. 
The left-hand figure shows the variation within the physical region for A, and 
shows that both S+_ and the direct asymmetry C + _ are largely independent 
of A. The right-hand figure shows the level of enhancement required to 
approach a 10% change in the value of the asymmetry, and that there is 
an increased effect for larger value of a. The effect is only significant for 
unrealistically large values of A. The current experimental results [101,102] 

S W7V = -0.30 ±0.17 ±0.03 (BaBar) 
S ww = -0.67 ±0.16 ±0.06 (Belle) 

can only be accommodated using very extreme values of A. For example, 
taking a model with same-sign fall-off, to reproduce the BaBar result with 
a = 2 would require a value of A > 10 to be within the la band and A = 20 
to approach the central value. This also produces unphysical values for the 
Gegenbauer moments, a 2 (2.2 GeV) = 7.7 and a 4 (2.2 GeV) = 0.5. Obtaining 
a value of S n7T = —0.67 is not possible for values of A > 1, which is required 
to keep a>2 positive. 
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Figure 3.7: CP asymmetries in B — > 7r + 7r - decays in QCD factorisation. Left: 

In physical region of A, (green) and (red). Right: as 

a function of A for a = 2, 3 ... 6, shown for models with same-sign 
fall-off (dashed) and alternating fall-off (solid lines). Curves converge 
as a increases. 



We now move on to study the effect of the model DAs on the B — » im 
branching fractions. We find that the effect is significantly more pronounced 
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on the branching ratios than for the CP asymmetries, as there is no longer the 
near cancellation of terms that occurs in the calculation of Acp- Considering 
first the decay of B — > tt + ?t~, the central value from QCD factorisation (for 
the asymptotic DA) is found to be 

BR(5 -> tt+tO = 5.5 x 10" 6 |0.25e i l5 ° + e" i7 | 2 (3.64) 

where 7 = 60° ± 7° [103], and the explicit dependence of the branching ratio 
on the Gegenbauer moments is again a polynomial in a n . Figure 3.8 shows 
the variation of the branching ratios for the model DAs compared with the 
current experimental data. These graphs show that large increases away 
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Figure 3.8: Left: Ratio of CP-averaged branching ratios for B — * tttt decays as 
compared with experimental result as a function of A, for model DA 
4>2 2'i B° — > tt + tt~ (blue) and B° — > 7t°7T° (green). Right: Example of 
T(B — > tt + vt _ ) (units 10 -6 ) as a function of A for a = 2,3. .. 5 for 
same-sign fall-off (dashed lines) and alternating (solid lines). Curves 
converge for increasing a. 

from the asymptotic value are possible by increasing A within its physical 
range. However, we see that in the case of 7r°7r° (where unlike 7t + ?t~, the 
enhancement is toward the experimental value) it is still not possible to 
reach the la error bound within the allowed range of A. 

This study does not aim to be a full error analysis and the graphs presented 
here are illustrative of the discrepancies that exist in the tttt system for QCD 
factorisation predictions. These results show that the lack of agreement 
between the predictions and the experimental results is not a result of the 
uncertainties from the leading-twist distribution amplitude, not even with 
our model DA with parameters well outside the allowed range. An alternative 
explanation of these discrepancies, namely the over-neglect of suppressed 
non-factorisable corrections, is discussed next in Chapter 4. 
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Non-factorisable corrections to 
charmless B-decays 

It does not do to leave a live dragon out of your calculations 

"The Hobbit", J.R.R. Tolkien 

This chapter is devoted to a study of the non-factorisable corrections to 
non-leptonic charmless decays of the 5-meson. This subset of decays is a 
crucial testing ground for QCD factorisation (QCDF). As we have access to 
increasingly accurate measurements in the flavour sector, it is essential to 
exploit the data to better understand the limitations and potential "pitfalls" 
with QCDF. Indeed, before we enter the era of LHC physics we ideally need 
to have complete grasp of the theoretical framework of B decays, as it is 
expected this sector will provide complementary, indirect evidence for new 
physics searches. As we discussed in Chapter 2, QCDF is only marginally 
consistent with the data on charmless B-decays, and investigation is war- 
ranted to study the source of this discrepancy. 

It is known that the factorisable contributions to exclusive B decays can be 
estimated using the formalism of QCD factorisation which is exact in the 
heavy quark limit, m& — > oo. Contributions that cannot be treated in this 
framework, i.e. the "non-factorisable" corrections, vanish in this limit, and 
are in general treated as unknown hadronic parameters. 

Our study specifically scrutinises the B — > tttt system in the context of these 
unknown non-factorisable corrections. A previous analysis of non-factorisable 
corrections was presented in [50], but we consider a different approach to the 
problem, and specifically quantify our comparison of the theory - data agree- 
ment. We begin by removing the model dependence from the QCDF pre- 
dictions and including "generic" non-factorisable (NF) contributions which 
aim to replicate both the known and unknown NF corrections. Our aim is 
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to determine the general size and nature of the NF corrections and test what 
exactly is required to reconcile the predictions with the experimental data. 

Our study analyses three separate scenarios for the inclusion of non-factorisable 
effects; in each scenario the sizes and phases of the NF contributions are var- 
ied. We also test a scenario which includes an enhanced "charming penguin" 
contribution, which we discuss in the following sections. We analyse the 
level of NF corrections required to bring each of the branching fractions and 
CP asymmetries within the B — > tttt system in line with their experimen- 
tal measurements. Finally, we give a short discussion on the implication of 
our results on the predictions for B — > irK, on which analysis is currently 
underway [104]. 

4.1 Non-factorisable effects in B — > tttt 

The continued improvement of experimental data on the charmless B decays, 
most notably for B — > mi, has culminated in a large enough number of results 
(with sufficient accuracy) to fully constrain all the QCD parameters relevant 
for describing this system. There are now six available measurements for the 
three B — > nn channels, for the branching ratios BR(£> — > 7r + 7r~) [105, 106], 
BR(B -> tt+tt ) [105,107] and BR(B -> [107,108], and for the three CP 

asymmetry measurements S+-, C + _ [101,102] and Coo [107,108]. The CP 
asymmetries are defined in equation (3.59). We use the most recent data (as 
of September 2005) from BaBar and Belle, as well as the combined results 
from the Heavy Flavour Averaging Group (HFAG) 1 . Although this system 
is well known experimentally, there are still some discrepancies between the 
two 5-factory experiments, most notably in the measurement of B — > 7r°7r°, 
and increasingly in the measurement of B — > 7r + 7r~. These are shown for 
illustration in Figure 4.1 and the CP asymmetries for 7r + 7r~ in Figure 4.2. 

As discussed in Section 2.4, there are discrepancies between these results and 
the QCDF predictions, leading to the coined U B — > titt problem" . In Chapter 
3 we examined the effect of the uncertainty in the hadronic input parame- 
ters - specifically that of the light-cone distribution amplitude of the pion 
- on the B — > nn observables. We determined that the resummed models 
can significantly affect the branching ratios and can account for both mod- 
erate and large deviations from the asymptotic DA. However, these models 
cannot cause enough enhancement to bring agreement between theory and 
data, within the physical range of the model input parameters. It is logical 
therefore, to examine the next large source of uncertainty in the theoretical 
predictions: non-factorisable corrections. 

1 http://www.slac.stanford.edu/xorg/hfag/index.html 
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Figure 4.1: Plot of experimental results for B — * ttit branching fractions in units 
of 10~ 6 : From left to right B -» 7r°vr (Yellow), and B -> vr+vr- 
(Magenta). HFAG error bounds are shown for comparison: BR(i? — > 
7r 7r°) = (1.45±0.29) x 1(T 6 and BR(B -> vr+Tr") = (5.0±0.4) x 10~ 6 . 
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Figure 4.2: Plot of experimental results for B — > 7T7r CP asymmetries from HFAG 
(Red), Belle (Blue) and BaBar (Green). The horizontal axis shows 
S n7T and the vertical axis C nn . 



The large quantity of experimental data makes this system ideal to test 
our hypothesis that non-factorisable contributions can have a large effect on 
predictions of charmless B decays. Annihilation contributions to B — > im 
are expected to be small in QCDF, and this is supported by attempts of 
direct calculation of these topologies via QCD sum rules [109]. If this is 
the case, then the enhanced NF contributions would be predominantly from 
hard-gluon exchange, as parameterised by the non-perturbative parameter 
Xh- This is supported to some degree by the results of factorisation fits 
to charmless data [110]. This study fits both and X& to all the ex- 
perimental data on charmless decays and finds best-fit values for Xh well 
outside the expected range. In our analysis we do not distinguish between 
the different sources of NF correction, instead categorising them according 
to their contribution to the B — > 7T7r isospin amplitudes. Using the decom- 



position presented in Section 2.3.3, we find that additional contributions will 
be required mainly for the AI = 1/2 amplitude. 

4.2 Charming penguins 

The charming penguin is a non-perturbative O (AQco/n^b) correction from 
enhanced penguin diagrams containing charm loops (Figure 4.3). It was first 
introduced by Ciuchini et al. in 1997 [111,112], who added long-distance 
contributions to decay amplitudes to improve agreement of fits to experi- 
ment. The charming penguin originates from a non-perturbative penguin 
contraction of the leading operators in the effective Hamiltonian, namely 
those with 0(1) Wilson coefficients Q\ and Q2, i.e., when the c and c an- 
nihilate. It is expected to give large contributions in some decay channels, 
notably B — > K + n~ and B + — > K°tc + , where the factorised amplitudes are 
colour or Cabbibo suppressed with respect to the penguins. 



The evidence and support for a charming penguin contribution has fluctu- 
ated since their introduction in [111]. The original parameterisation was in 
terms of a diagrammatic deconstruction into renormalisation group invariant 
quantities or topologies [113]. For example the leading (tree) contributions 
are described by emission topologies Ei and E 2 which are scheme indepen- 
dent combinations of Wilson coefficients and matrix elements related to the 
current-current operators Qi and Q2, as shown in Figure 4.4. 

The charm and GIM penguin followed a topology of the type in Figure 4.5. 
The charming penguin has only the charm quark in the loop; the GIM pen- 
guins are those where the diagrams always appear in the combination of u — c. 

This description was later updated to include further diagrams (annihilation 
and penguin contractions) in addition to the original diagram, to give two 
phenomenological parameters P\ and pp IM for the charming and GIM pen- 
guins respectively. Fits of the charming penguins to Kit and ittt decays gave 
promising results [114]. However, when QCD factorisation was introduced 
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Figure 4.3: Example of a charming penguin. 
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Figure 4.4: Emission diagrams. Left: colour allowed (DE) topology; Right: 
colour suppressed (CE) topology. The dashed line represents the 
four-fermion operator. 



Figure 4.5: Example of charming and GIM penguin topologies (CP). 



and the non-factorisable corrections were shown to be calculable in pertur- 
bation theory, the charming penguin fell out of favour. It should be noted 
that since the charm penguins are expected to have a size equivalent to a 
AQCD/fft'b correction they are not in fact in conflict with the QCDF results in 
the rrib — > oo limit. If we call into question the size of the power-suppressed 
corrections to the AQ CD /mfc expansion and expect - and in fact demand - 
some of the non-factorisable contributions to be large, the charming penguin 
can provide significant additional contributions at this order. 

The current QCD factorisation predictions (specifically for the B — > tttt 
and B — > ttK decays) do not satisfactorily explain the current experimental 
data. Both the B — > tttt and B — > ttK decays may require significant levels of 
power-suppressed corrections. It is therefore quite possible that the charming 
penguin can provide additional enhancement, enough to bring the predictions 
into agreement with the current experimental results. 

Charming penguins are not Cabbibo enhanced in B — > tttt decays, unlike 
the case for B — > ttK. It is naively expected that, along with the other 
O (AQco/^ib) corrections, the charming penguin contributions should be 
small. However, as we show in our study, the predictions for the B — > tttt 
branching fractions cannot be brought into agreement with the experimen- 
tal data without significant non-factorisable corrections. Since the charming 
penguin is of the same order as the O (Aqcd / ' m b) effects, there is no reason 
why a sizable contribution should not exist in these decays. 



B 




fifi 



The charming penguin parameter Pi, and the GIM-penguin parameter p^ 1M 
were expressed as a complex number fitted to experiment by Ciuchini et al., 
and added to the factorisable amplitudes via 

Pi = a\A^ + P 1 P? 1M = (at - al)A^ + P? m 

where the notation is from QCD factorisation, as found in Appendix B. 

There is some disagreement over the inclusion of charming penguins within 
the factorisation framework. The parameter Pi contains not only the charm- 
ing penguin contribution but additionally the annihilation and penguin con- 
tractions of the penguin operators. Ciuchini et al. suggest that all chirally 
suppressed terms should be dropped and replaced with this term (these cor- 
respond to dropping the coefficients ag from QCD factorisation). The propo- 
nents of QCD factorisation however, argue that the "charming penguin" con- 
tribution is completely described by the penguin annihilation contributions 
and non-perturbative corrections to the coefficient at , i.e those contained 
in dg. They argue that a large non-perturbative penguin enhancement is 
implausible, and that any modification required to bring theory predictions 
into agreement with data should be attributed to weak annihilation rather 
than charming penguins [2]. 

As our study does not include any of the model-dependent calculations of 
the C(AQ CD /m b ) terms suggested by BBNS, to obtain a charming penguin 
enhancement we need additional contributions. We do not replace all of 
the chirally- suppressed terms as suggested by Ciuchini et al., and instead 
take the charm penguin to be in addition to the QCDF penguin diagrams 
(i.e. keep the calculable part of cig, but remove any dependence on terms 
proportional to Xh)- We can then allow for the possibility of the charming 
penguin explicitly as an enhancement of the coefficient at = at + r x a%. 



4.3 B —> nn Analysis 

We use the expansion of the B — > ixix decay amplitudes in terms of the two 
isospin amplitudes Ay 2 an d A 3 / 2 . We construct three different scenarios for 
inclusion of non-factorisable corrections and enhancements on top of the fac- 
torisable amplitude as calculated within QCD factorisation. We include only 
the (fully) factorisable parts from QCDF, i.e. we do not include the hard- 
scattering or annihilation contributions incorporated in the model-dependent 
parameters X H and X A (Section 2.3.2). NF corrections are then added ac- 
cording to the prescriptions of our three scenarios, so we can determine the 
size and nature of corrections needed to reconcile with the experimental data. 

In each scenario we allow the size and phases of the NF contributions to each 
isospin amplitude to vary independently. The full set of B — > mr branching 
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ratios and CP asymmetries are calculated for each of the combinations of 
NF contributions we consider. We take a defined set of input parameters 
which are also varied within their allowed ranges in order to further optimise 
the agreement between theory and experiment. Our input parameters are 
a combination of CKM phases, the B — > 7r form factor and A from the 
parameterisation of the pion wavefunction using the models developed in 
Chapter 3. 





Input ranges 


A 


1.1< A < 1.3 




0.26 ± 0.04 


R u 


0.399 ± 0.08 


7 


60.3°± 6.8° 



We consider primarily the data from the HFAG which averages over all of the 
current data available for each channel. However, we consider the data from 
BaBar and Belle separately where it is required, if the discrepancy between 
the two experiments is significant. Each of our scenarios is constructed with a 
"base" of the fully calculated factorisable part of the B — > tttt isospin ampli- 
tudes, as found using from QCDF. Different combinations of non-factorisable 
corrections and charming penguins are then added. 

We add a percentage non-factorisable contribution to each of the isospin 
amplitudes individually. The amplitude A 3 / 2 is tree dominated with a tiny 
contribution from electroweak penguins; the expressions for the isospin am- 
plitudes given in equation (2.37) show that the amplitude A^, 2 contains only 
contributions from electroweak penguins. We can therefore safely assume 
that there will be no significant non-factorisable contribution to this partic- 
ular amplitude. This is supported by the good agreement between the QCD 
factorisation prediction and the measurement for the decay B — > 7r + 7r°. 

The non-factorisable corrections are added then as 

A 3/2 = A F 3/2 + \* u \A u 3 f 2 \N 3/2 e id ^ 

A 1/2 = A( /2 + X:\A^ 2 \N^ /2 e id - + K\Alf 2 \N c 1/2 e^ (4.1) 

where the notation is intuitive: the superscript F denotes the pure factoris- 
able amplitudes, the N^i e ldAI is the size of the non-factorisable contribution 
(0 < Nai < 1), with arbitrary phase d^i. The level of non-factorisable cor- 
rections is expected to be of the order of A/to& ~ 10%, where A measures the 
kinetic energy of the b quark within the B meson. As a conservative estimate, 
we allow corrections of up to 20% as our "normal" expected level, i.e taking 
each N < 0.2. As an indicator for the enhanced corrections we allow double 
this estimate with up to 40% non-factorisable corrections. Anything beyond 



this limit would be approaching the point where the factorisation scheme 
itself would begin to break down. The QCD factorisation scheme is based 
upon the calculation of order a s corrections, so any terms suppressed by l/m& 
should not be much larger than a s rrib, otherwise it would not make sense to 
include small calculable terms and to neglect large incalculable terms. 

We also consider the possibility of additional contributions from charming 
penguin diagrams, which are not accounted for in the calculation of the fac- 
torisable amplitudes. The charming penguin contribution arises purely in the 
isospin amplitude A\ , 2 , and as such, can be parameterised as an enhancement 
to the factorisation coefficient a c A 



where R is some complex factor R — \R\e lip with \R\ > 1. This will have an 
effect on both the decay amplitudes for BR(B — > it + it~) and BR(B — > 7r°7r°), 
but there is no effect on BR(B — > tt + tt°) which has no contribution from Ay 2 . 

The scenarios are summarised as follows: 

• Scenario I: QCD factorisation + expected level of non-factorisable cor- 
rections at < 20% 

• Scenario II: QCD factorisation + enhanced non-factorisable corrections 
at < 40% 

• Scenario III: QCD factorisation + expected non-factorisable corrections 
at < 20% and charming penguin contribution 

In order to quantify the "quality" of each of these scenarios we construct a 
X 2 statistic, where the form factor F^^iO) = F n and the phases of the non- 
factorisable contributions (NF) are left as fitting parameters. The absolute 
size of the NF correction and the CKM parameters R u and 7 are taken as 
fixed, and the fitting is performed separately over differing values of the pion 
distribution amplitude parameter A. We define 



where yi are the theory predictions, and Xi are the data values with their 
associated error <jj, calculated by combining statistical and systematic er- 
rors in quadrature. The index % runs over the 6 available experimental 
measurements of the CP-averaged branching ratios and CP asymmetries: 
{r + ~, r 00 , r +0 , S+-, C + _, C o}. We then numerically minimise the x 2 while 
scanning over the three phases and applying the constraint F n < 0.3. We 
perform this separately for the BaBar, Belle and HFAG averaged data sets. 



«4 = a\ + r x a% 
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The goodness-of-fit then follows the x 2 probability distribution with 2 de- 
grees of freedom (equal to the number of measurements minus the number of 
fitted parameters). We express the goodness-of-fit as x 2 / n d, where the closer 
the value to 1, the better the fit. 

4.4 Results 

We now present and discuss the results of the analysis for the B — > -kit system 
for each of the specified scenarios in turn. 

4.4.1 Scenario I 

In this scenario we calculate the CP-averaged branching ratios and CP asym- 
metries for B — > 7T7T decays with the expected level of non-factorisable cor- 
rections (< 20%). We find that it is very difficult to reconcile the predictions 
for all of the observables with their experimental measurements within the 
experimental error bounds. We varied all of the input parameters within 
their allowed ranges and found a set of optimum inputs (shown in Table 4.1) 
which gave the best agreement between theory and experiment; the results 
were approximately the same for all data sets. The two CKM parameters op- 
timised at their input values but the DA parameter A required its maximum 
allowed contribution to best optimise the agreement. For the form factor we 
used the best-fit value as calculated from the minimisation of the x 2 function 
on the full set of B — > nn data. 



A 




R u 


7 


1.3 


0.22 


0.399 


60.3° 



Table 4.1: O ptimum input parameter; R u , 7 are the default central values (Ru = 
Vp 2 + V 2 ) for 20% non-factorisable corrections. 

Since the optimum scenario is clearly obtained for the maximum possible NF 
correction, we specified the size of the contributions as iV = 0.2 for each am- 
plitude and scanned over all possible phase combinations. We then calculated 
the ratio of the theory over experimental results for the branching fractions 
and the CP asymmetries. The most illustrative results are of B — > 7r + 7r~ 
against B — > 7r°7r° and for S wn against C wn . Considering first the branching 
ratios, the value of BR(B — > ir + ir°) is quite well reproduced, but the ir + ir~ 
and 7r°7r° branching fractions show considerable discrepancy. This is illus- 
trated for the HFAG result in Figure 4.6 and also separately for the data sets 
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from BaBar and Belle in Figure 4.7 and Figure 4.8 respectively, highlighting 
the large impact of the differing measurements of BK(B — > 7r + 7r~). 

The ellipses shown in Figures 4.6 to 4.8, are the 2a and 3a error ellipses. 
These are calculated from the experimental results under the assumption of 
no correlation between the branching ratios, and using the measured cor- 
relations for the CP asymmetries [101, 102]. These curves correspond to 
confidence regions (confidence level across 2 degrees of freedom) of ~ 84% 
and ~ 95% respectively. 
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Figure 4.6: Theory/Experimental ratio using HFAG data for against T 00 
(left) and CP asymmetry S nn /C nn (right). Error ellipses shown at 
2a and 3a. 
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Figure 4.7: As Figure 4.6 using only BaBar data set. 
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Figure 4.8: As Figure 4.6 using Belle data set. 



These figures clearly show that with a 20% non-factorisable correction, it is 
not possible to reproduce the experimental data within the 2a error bound, 
and not even to 3a using the Belle data alone. Using the HFAG data and 
maximal contribution of A = 1.3, the x 2 minimisation yields a value of 
F n = 0.22 with x 2 / r of 4.41/2. For the BaBar data alone we also obtain 
F n = 0.22 but with a much improved fit and x 2 / r = 0.850/2. 
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We observe that in increasing the parameter A from the asymptotic value 
(A = 1) to its physically allowed maximum (A = 1.3) we find a steadily 
decreasing x 2 value, and correspondingly increasing best-fit value of F n . The 
distribution of the x 2 values as a function of the form factor is plotted for 
20% NF corrections, for the full range of possible phases in Figure 4.9. This 
shows clearly the minimum \ 2 m agreement with the numerical minimisation 
requiring F w ~ 0.22. 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 


20 40 60 80 100 
Figure 4.9: Dependence of minimum x 2 on F w for A = 1 and 20% NF correction. 



4.4.2 Scenario II 

In a direct extension of Scenario I, we repeat our study with a further en- 
hanced non-factorisable contribution beyond that of the expected levels. We 
again use the factorisable isospin amplitudes calculated within QCDF, but 
now add up to 40% non-factorisable corrections. As may be expected, this 
reproduces the experimental data quite well. We find the optimum param- 
eter set to be similar to that from Scenario I; the agreement is not greatly 
improved from changing the CKM input, but it is improved by the maximi- 
sation of A and use of the best-fit form factor. For a 40% NF correction, 
the x 2 minimisation gives a value of F n = 0.23 at x 2 / r = 2.0/2 when us- 
ing the maximum correction of A = 1.3. As in Scenario I, we see that the 
X 2 improves as we increase A, and the best-fit form factor increases corre- 
spondingly, around the value of F 77 ~ 0.23. This is illustrated in Figure 4.10. 

We find that there is a range of allowed values of from 0.21 < F w < 0.23 
which give good agreement with the data. We illustrate the best possible 
agreement (for both a 30% and 40% non-factorisable correction) using the 
BaBar data set in Figure 4.11. 

The results using the HFAG data are shown in Figure 4.12. In this case the 
experimental results are further away from the theory predictions, primarily 
due to the Belle measurement of BR(B — > 7r°7r°). We see that the enhanced 
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Figure 4.10: Dependence of minimum x 2 on F n with increasing A for 
correction. 
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Figure 4.11: The 55(5 vr+vr"), 55(5 vr 7r°) agreement for the best- 
fit case for 30% NF correction (left) and 40% (right), both using 
BaBar data. 
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Figure 4.12: The 572(5 -> vr+vr"), 55(5 -> 7r°vr ) agreement for the best-fit 
case for 30% NF correction (left) and 40% (right), using HFAG 
data. 



levels of NF corrections of around 30% can bring the theory values close to 
the experimental averages, at least to within the 3a error, although larger 
corrections are needed to approach 2a agreement. The optimum for A is 
again found to be maximal, A = 1.3, and for the form factor F n = 0.22 and 
F 71 " = 0.23 for the 30% and 40% cases respectively. 
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4.4.3 Scenario III 



Our third scenario involves adding contributions from charming penguins in 
addition to the expected non-factorisable correction (< 20%) to the QCDF 
base. As discussed above, we expect this to impact on the values of the 
branching ratios for B — > and B — > n + ii . For illustration, Figure 4.13 
shows the possible enhancement of these branching ratios for various sizes of 
charming penguin contribution. 




Figure 4.13: Dependence of BR(B -> vr+vr") (left) and BR(B -> 7r°7r°) (right) 
on charming penguin parameter R, in units of 10~ 6 ; curves shown 
for constant if of (long-dashed), 7r/4 (short-dashed), tt/2 (dotted) 
and 7r (solid). The experimental result and la error from HFAG is 
shown for comparison. 

This is presented for the expected 20% non-factorisable correction and our 
"best-fit" set of parameters (A, F n and NF phases for each branching ratio), 
as found in the analysis with no charming penguin contribution (Scenario I). 

The figure shows that it should indeed be possible to improve the agree- 
ment between the theory and experimental results for both branching ratios 
simultaneously, with a level of 20% non-factorisable corrections. 

To better gauge and quantify the level of enhancement that is needed, we 
again compare the agreement of the ratio of B — > 7r + 7r~ and B — > 7r°7r° 
branching fractions, and the B — > 7r + 7r~ CP asymmetry. The results are 
shown below in Figure 4.14 and Figure 4.15 for sample values of \R\ = 
1, 1.5,2, scanning over possible phases < ip < 2tt and using the combined 
HFAG dataset. 

Similarly to our previous scenarios, we find that the input parameters opti- 
mise for the maximal contribution A = 1.3 and with F n = 0.22. The figures 
show that considerable enhancement is possible, even for small values of \R\, 
and that it is possible to reproduce the data very well. It is clear however 
that there must be restrictions on the size and the phase of the parameter R, 
in order not to overshoot the experimental bounds - especially in the case 
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Figure 4.14: Theory/Experimental ratio using HFAG data for against T 00 
with charming penguin contribution of | R\ = 1 (magenta), \R\ = 1.5 
(yellow) and \R\ = 2 (Black), scanning over all possible phases <p. 
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Figure 4.15: As above for CP asymmetry S^/C- 



of the CP asymmetry. Any contribution of magnitude \R\ > 2 would be 
difficult to reconcile with the data. 

As can be seen in the figures, using any magnitude of enhancement between 
\R\ = 1 and \R\ = 2 can reproduce the data extremely well for both the 
branching ratios and for the CP asymmetry; however, the required phase does 
not agree for both. The best-fit for the branching ratios requires ip ~ 120°, 
and for the CP asymmetries ip ~ 300°. It is the value of C^, that dic- 
tates the best-fit value for the CP plot, as it is very sensitive to the phase 
of the additional charming penguin contributions. Using our x 2 minimisa- 
tion technique we found the best-fit to all six experimental results, requiring 
a charming penguin contribution of R = 2e* 260 °, resulting in a (minimum) 
X 2 = 0.44; this is shown below in Figure 4.16. 
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Figure 4.16: Best fit agreement from x 2 minimisation to full data set with charm- 
ing penguin contribution of R = 2e* 260 , using HFAG data with 
A = 1.3, F n = 0.22. 



4.5 Discussion and comments 

In this chapter we have presented an analysis of the non-factorisable correc- 
tions to B — > 7T7T and shown how the experimental data can be used to guide 
and constrain the levels of contribution needed to gain agreement between 
data and theory. Before discussing the full implications of our results, it is 
useful to compare the size of the model-dependent QCD factorisation pa- 
rameterisation of the non-factorisable effects to our generic non-factorisable 
contributions. Consider first the corrections arising from the hard-spectator 
diagrams, where the divergences from the chirally enhanced twist-3 contri- 
butions to kernel T t n are parameterised by complex parameter Xjj: 



HAim) 



m B \ B F£(0) 



SrlAX H 



with Ni = (i = 6, 8), -1 (i = 5, 7) or 1 otherwise. These contributions 
enter into each of the factorisation coefficients af^mr). We can extract the 
contributions proportional to Xh in both of the isospin amplitudes Ai/% and 
A3/2 and obtain 

/ ap n r rn B C F 7ia s ((x h ) 

{A 1/2 )hs - B^— — 1 l 2 x H 

(AP \ R mB C F™s{Vh) r ia y 

K A 3/2) hs - -B^- ^2 L3/2 r x AX H 



with B W7T = i^f B fl, and the Caj are combinations of Wilson coefficients 
given as 

Ci/ 2 = 5 pu (4C 2 - 2d) + 6C3 - 3C 8 + 3C 9 - 3C W 
C 3/2 = 5 pu (2C 2 + 2d) + 3C 8 + 3C 9 + 3C 10 

Comparing this contribution to the expression for our generic non-factorisable 
correction (4.1), we can estimate what size of correction corresponds to the 
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"default" estimate of X H ~ \n(mb/A h ) ~ 2.4 [2]. For A = 1 and equal size 
contributions from Ny 2 and N^ 2 , we estimate that using the BBNS model 
for X H equates to a non-factorisable correction of ~ 8%. 

We can perform a similar decomposition for the other main source of non- 
factorisable corrections - the annihilation contributions - parameterised by 
the parameter X A in the BBNS model. There is no correction to the A 3 / 2 
isospin amplitude (and correspondingly no annihilation contribution to the 
B — > 7T+7T decay), and the correction to Ay 2 can be expressed 

2C 



with Cj again combinations of Wilson coefficients 

Ci = S pu Ci + C3 + 2C4 + 2Cq — ^Cg + |Cio + \C% 
C 3 = C 5 + N C C 6 - \C 7 - \N C C % 

and the amplitudes A\ and A$, calculated with asymptotic distribution am- 
plitudes and assuming SU(3) flavour symmetry, are 



A\ = 2na s 



9(X A -A + —)+rlX 2 A 



3 

A{ = l27ca s rl(2X 2 A -X A ) 

We can see the explicit power-suppression of these terms relative to the 
leading gluon exchange contributions via the fact that the weak annihilation 
terms are proportional to f B rather than fg Trig / \b- By comparing these 
expressions with the isospin expansions of (4.1) and using the default value 
of X A ~ \n(m b /A h ) ~ 2.4, we obtain a corresponding size for the generic 
non-factorisable correction of ~ 5%. 

From these estimates it becomes very clear how this size of contribution can- 
not reproduce the experimental data. This would correspond to an overall NF 
contribution less than the moderate (20%) correction we tested in Scenario 
I, which in itself could not reconcile the predictions with the experimental 
data. 

It was the marginal agreement of the QCDF predictions with the B — > -kit 
data that originally motivated us to test the consequences of enhancing the 
non-factorisable corrections to a moderate (20%) and extreme (40%) level. 
Our results show clearly that even for extremal values of our input param- 
eters, the best-fit scenarios do not provide reasonable agreement between 
theory and experiment until at least 30% contributions are used. In the 
case of the combined HFAG data, a 40% contribution is required to obtain 
a reasonable 2a agreement. This also highlights the problems with the dis- 
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crepancies between the BaBar and Belle measurements. 

Our final scenario incorporated the possibility of the charming penguin as 
an additional non-factorisable effect. We found that significant enhancement 
of the branching ratios is possible, but that a considerable contribution with 
\R\ > 2 is not supported by the data. There are also a number of questions 
which must be addressed if we are to consider this contribution as a realistic 
contender for a significant, neglected l/m b correction. Firstly, the required 
best-fit phase of the charming penguin parameter R does not naturally agree 
for both the branching ratios and the CP asymmetries simultaneously. The 
best-fit agreement (shown in Figure 4.16) is however reasonable, although 
yet more NF corrections would be needed to gain a la agreement with the 
data. 

Secondly, we must ensure that the size of the contribution suggested by our 
results is appropriate with respect to the leading QCD factorisation terms. 
To check this, we can determine the relative size of the charming penguin 
contribution to the relevant isospin amplitude, A 1 / 2 , compared with the fac- 
torisable contribution. For the maximum allowed charming penguin contri- 
bution \R\ =2, the ratio of the charming penguin contribution to the QCDF 
contribution is ~ 0.3, and for a moderate contribution of \R\ = 1.5, we find 
a ratio of ~ 0.15. We can compare this to the relative size of our generic 
non-factorisable corrections, for example with equal phases and | = 0.2, 
which provides a correction 15% the size of the purely factorisable terms. 

It would seem that the enhancement from the charming penguin must be ap- 
proached with caution, as even a moderate increase, \R\ = 1.5, gives an over- 
all contribution equal to the size of the non-factorisable corrections coming 
from all of the annihilation and hard-spectator topologies. This suggests the 
implementation of the charming penguins in this way can provide a relative 
size contribution too large for them to be realistic additional NF contribu- 
tions. We should then increase the constraint on the value of R, restricting 
the maximum contribution to be \R\ ~ 1.2. The theoretical agreement to 
the B — > 7T7T observables would still receive significant improvement from the 
additional contributions, although not enough to bring both the branching 
ratios and CP asymmetries into 1o agreement using only a 20% generic NF 
contribution. The best-fit agreement to the HFAG data with this constraint 
is shown in Figure 4.17. 

Finally we can consider the application of this work to other charmless de- 
cays, as the techniques and results that we determined for B — > tttt can be 
applied to B — > irK [104]. An isospin decomposition can be performed in an 
analogous way to B — > nit; there are now six independent isospin amplitudes 

AP AP AP 

3/2,1 1/2,1 1/2,0 

with p — u,c. The factorisable contribution can again be determined from 
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Figure 4.17: Best fit agreement from x 2 minimisation to full data set with charm- 
ing penguin contribution of R = 1.2e* 269 , using HFAG data with 
A = 1.3, F n = 0.22. 



QCDF, and the non-factorisable contributions defined as corrections to the 
above amplitudes. These corrections can be related to the non-factorisable 
corrections to B — > tttt via SU(3) flavour symmetry. 

It is possible in principle to use the best-fit data found from the B — > Tin 
case for A(tv) and F^^ n . The separate txK input parameters, A(K) and af 
(from the <pK distribution amplitude), F^ K , and the CKM factors R u and 
7, can be varied within their allowed ranges. 

There is experimental information available for four ttK branching ratios 
(tt°K + , 7t + K°, tt°K°, it~K + ), four direct CP asymmetries C + o,Co + ,Cqo and 
C_ + and one indirect CP asymmetry for B — > n°K . Using this information, 
the level of non-factorisable corrections needed to bring the predictions for the 
ttK system into agreement with the data can be tested and constrained. The 
possibility of charming penguins also has more significance in this system, as 
they are expected to have a significant contribution to the decays B — > K + n~ 
and B+ -> K°tt+. 
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Chapter 5 

Exclusive radiative B-decays 



It is wisdom to recognise necessity when all other courses have 
been weighed, though as folly it may appear to those who cling to 
false hope... 

"The Lord of the Rings", J.R.R. Tolkien 

Radiative penguin transitions such as b — > S7 or b — > are examples of 
flavour changing neutral currents. They are rare decays, arising only at 1- 
loop in the Standard Model and are a valuable test of predictions of flavour 
physics, enabling us to test the detailed structure of this sector at the level 
of radiative corrections. 

There has been much work and speculation regarding the phenomenology of 
FCNC processes, as they are ideal candidates for indirect searches for new 
physics phenomena e.g. as reviewed in [115,116]. The process b — > sj allows 
a large parameter space for new physics and has undergone scrutiny by a 
number of authors [117-121]; this process is seen as ideal to constrain the 
possible new physics effects, assisted by available experimental data on both 
the inclusive process B — > X s 7 and the exclusive channel B — > K*^. The 
b — > g?7 process is suppressed with respect to b — > S7 by the CKM factor 
\ — VtbVtdi an d the exclusive modes have previously had only experimental 
upper bounds from the B factories. In June of this year however, Belle 
announced new measurements of B — > p7 and B — > [122]. These new data 
for b — > d decays and updated results for the b — > s processes present an ideal 
opportunity to perform an analysis of the radiative B decays, allowing us to 
both test the predictions of the Standard Model and examine the sensitivity 
of b — > (d, s)7 to new physics. 

In this chapter, we consider the exclusive channels B — > (p, u)^ and B — > 
-fT*7 within the Standard Model, based on the formulation of B — > V'j decays 
in QCD factorisation from Bosch and Buchalla [4]. We construct ratios of 



these branching fractions in order to limit the theoretical uncertainties in 
the predictions. We also examine the effect of the uncertainty from the 
distribution amplitudes of the vector mesons - the p, uj or K* - by using 
our new resummed DA model as developed in Chapter 3. By testing the 
branching fractions against our DA parameter A, we can determine how 
much of an effect the LCDA has upon the final predictions for these decay 
channels. Using the ratio of the B — > (p, uj)^ and B — > K*^ branching 
fractions we show how the CKM ratio | Vtd/V^ s | can be extracted, and obtain 
an estimate of this quantity using the QCDF predictions. 

We then go on to motivate an analysis of B — > V'y decays within generic su- 
persymmetric models with the so-called mass insertion approximation [123]. 
We then determine the constraints on such a generic MSSM model for both 
b — *■ g?7 and b — > S7 transitions, using all the available experimental data 
from the B factory experiments. 

5.1 B —> Vy in QCD factorisation 

A model-independent framework for the analysis of the radiative decays has 
been presented in [4]; we summarise the important points and calculational 
formulae here. The method presents an expansion in the heavy-quark limit, 
giving expressions valid at leading order in a Aqcd / f^b expansion. There are 
contributions from power-suppressed annihilation topologies, particularly rel- 
evant for the B — > p7 decays, where they are numerically enhanced by large 
Wilson coefficients. We implement these corrections following the method 
of [124]. 

The process, shown in Figure 5.1, is dominated by the electromagnetic pen- 
guin operator Qi^. The annihilation diagram shown is important in B — > fry, 
specifically for the charged decays. 
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The decay amplitudes are expressed following a similar procedure as in QCDF 
for decays into two mesons, i.e. in terms of Wilson coefficients and hadronic 



SI 



matrix elements, which are then calculated using a factorisation formula. The 
matrix elements for the radiative processes are found using a factorisation 
formula valid up to corrections of order (Aqcd/^6) 



(V-r{e)\Qi\B) 



Jo 



e (5.1) 



where e is the photon polarisation four-vector. For the matrix element of the 
operator Q 7 , the formula (5.1) becomes trivial: The kernel T? is a purely 
kinematic function and the spectator term is absent, so the only dependence 
on non-perturbative input parameters that remain is from the B — > V tran- 
sition form factor. At the leading order, this is the only contribution to the 
B — > V7 decay. At the next-to-leading order we gain additional contributions 
from the other AB = 1 operators, dominated by Q\ and Q$- 

The matrix elements of the additional operators can be expressed in terms 
of the calculable matrix element for (Q7). There are two contributions, one 
from the hard-scattering kernel T/ for vertex corrections, and the other from 
T t H which includes the hard-spectator corrections. These are given by 

The type-I contributions depend on Gj(s p ), which are simple functions of 
s p = nip /ml (for p — u,c) given explicitly in [4]. The spectator contributions 
are incorporated into the functions Hi, and are more interesting than the 
type-I functions: 



ffi(s p ) = 



2tt 2 f B ft n 



T^iMil f 1 dvh(l-V,S p )$ ± (v) 

Jo s Jo 



3N c F v m 2 B Jo 
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The hard-scattering function h(u, z) is given as 



h(u, z) = 



4z 



1 - 



u—4z 



u 



1 + 



u—4z 



2 

u 



(5.2) 



and is real for values of u < Az, and complex for u > Az. For z — > 
(for example the up-quark contribution z u ~ 0), the function is regular and 
becomes 

2 

fr( u , 0) = — 
u 
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The distribution amplitude for the 5-meson is parameterised as before using 
the first inverse moment 



1 d* * B ® = — 

£ 

These type-II contributions also include a dependence on the light-cone dis- 
tribution amplitude of the vector meson of the B — > Vj decay, specifically, for 
the transversely polarised vector mesons, <p±. The vector meson is described 
at leading twist, by two distribution amplitudes <f>\\ and <p± corresponding 
to longitudinally and transversely polarised mesons. The vector meson in 
B — > V'-y decays is to leading power transversely polarised, and hence can be 
described by <p±, as any contributions from 0y are suppressed. 

As with the pseudoscalar mesons, the distribution amplitude <p± can be 
treated using a conformal expansion that is often truncated at the moment 
a-i. In our analysis, we wish to examine the effects of a non-standard DA 
on the radiative decays, so we use our resummed model DA in place of this 
truncated expansion. For definiteness we use the model 0^ 3 (A), where A = 1 
reduces to the asymptotic DA. This model applies to the mesons symmet- 
ric in the fractional momentum u, i.e. it, p, uj. Where we consider the K* 
channels, there is a non-zero anti-symmetric part to the meson wavefunction 
which we can also model using our resummed DA models. This is constructed 
in an analogous way to the symmetric DA, as discussed in Chapter 3. We 
use the model tp^(ai), which allows for the inclusion of an arbitrary value of 
the first anti-symmetric Gegenbauer moment af . We take the value of af 
from our discussion on <px unless otherwise stated. The value of A for the 
K* wavefunction is the sum of the symmetric and antisymmetric parts, i.e. 

^tot,+ — a _|_ A asym ' + 

with A as y m -+ = - ai (3/2) c C(c, 3/2). Hereafter, we will refer to the total A tot 
as A and imply the inclusion of the anti-symmetric part. 

Finally we combine all of the contributions to the B — > V7 decay into a 
factorisation coefficient, written as [4] 

o?(v 7 ) = + [cr\m b ) Gl ( Sp ) + c r s (m b ) G S 
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CfVfc)#i(sp) + ^ )e V)#8 (5.3) 



The Wilson coefficients C[ ^ and C^ cS as indicated, are taken at leading 
order. 

The decay amplitudes are then expressed as follows, keeping both the charm 
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and up-quark contribution 



A(B -> Vj) 
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(V 7 |07|B> 



(5.4) 



giving the final B — > V7 branching ratios as 



GpCteTn 



mi 



32tt 4 



ml 



where q 



for V = and q = d for V = p; Cy 



1 for K\ 



(5.5) 
and 



cy = 1/V2 for p°. The CP-conjugate branching ratios are obtained by 



replacing 
averaged. 



X^"; in our analysis we take all branching ratios as CP- 



The annihilation contributions can be included in the factorisation coefficient 
dy via 



a: 



a 7 + <nn 



The dominant annihilation contributions are suppressed by one power of 
AQCD/^bj but can still be calculated within QCD factorisation, (a proof to 
0(a s ) can be found in [124]). These can give a sizable contribution to the 
radiative decays, such as B — > pj, where the contributions originating from 
operators Q\ and Q 2 are enhanced by large Wilson coefficients C\ and C 2 . 
These have small impact on the B — ^ K* decays due to CKM suppression, 
however this decay can receive a large contribution from the operator Q e 
[125]. The annihilation contributions are also sensitive to the flavour of 
the light quark of the 5-meson, and so differentiate between decays of the 
B° (B°) and the B ± , and additionally are dependent on the distribution 
amplitude (fi±. Since this is clearly important to our analysis, we describe 
the contributions a v ann in more detail in Appendix B. 

There is no indication for large power corrections beyond these calculable 
annihilation terms [4,126]. 



5.2 Analysis and comparison with Belle data 

For channels that have low statistics such as the rare radiative processes, the 
experimental measurements often combine different individual channels in 
order to improve the significance, and allow a signal to be discovered. Belle 
is to date the only experiment to have released an actual measurement of the 
b — > g?7 processes [122], where they combine the data from three channels 



84 



p°7, p + 7, ^7 . The average is defined as [127] 



BR[ J B^(p,cu) 7 ] 



BR(B 



+ -> p+ 7 ) + ^1 rBR( J B° -> p° 7 ) + BR(/7° -> 



We can also define an isospin-averaged branching ratio BR(B — > i^*7) which 
is used frequently 



BR(B -> iT 7 ) 



I {br (5 - 



iT+ 7 ) + ^BR(£° -> iT 7 ) |> (5.6) 



We can predict the branching fractions or CP asymmetries for the B — > 
decays directly from these averages using the QCD factorisation formulae of 
Section 5.1. We can also find the ratio of these averaged branching fractions 
which is fitted directly from the experimental data. We define, using the 
CP-averaged branching fractions [127, 128] 



BR [B^(p,w)i\ 



BR(B 



V, 



is 



2 C 1 " m l P ,u>)l m lf 



:i 



m 2 K ./m 2 B ) 3 



Col 



a c 7 (K*>y) 



[1 + AR] 
(5.7) 



C denotes the form factor ratio T 1 p (0)/T 1 K *(0), which characterises the 577(3) 
breaking in the transition form factors. The form factor 71(0) is one of seven 
independent form factors which describe all B — > V decays, and is the only 
one relevant for the radiative B — > V'-f decay. It is defined in conjunction 
with two other form factors T 2 3 via 



1 

cv 



(V(p)\q<T lu ,q v (l + >y 5 )b\B( PB ))=ie 



+ T 2 (q 2 ){e; 



m, 



m 



s^p^T^q 2 ) 
v) -e*q)(p B +p),,} 



+ T 3 (q 2 )(e*q){q„- 



m B — ITLy 



The parameter ( has been estimated using several methods, the most recent 
value is from Ball and Zwicky using light-cone sum rules [129] who find 
= 1.25 ±0.18. 



AT? is given within QCD factorisation by [128] 



AR = 2Re 5a 



R\ - R u cos 7 
l-2.Ru cos7 + i?2 



(5.8) 



written in terms of R u = \J p 2 + r] 2 and the CKM angle 7. 5 a is given by 



5a = 



a "((p^)7) -a£((p,q;)7) 



(5.9) 



where the factors a? implicitly include the annihilation contributions, which 
contribute only at order If the annihilation contributions to 5a 

were neglected then 5 = 0(a s ). To obtain expression (5.9) we make use of 
the identity: 



where the second term in the bracket in (5.10) is neglected for B — > K*^ as 
it provides a contribution of less than 1%. This results in the correction 5a 
depending only on the factorisation coefficients from (p, uj)'~f. 

The advantage of calculating the ratio of branching fractions in this way is 
to reduce the error from theoretical uncertainties. Firstly, from the theory 
parameters that cancel in taking the ratio, including parameters such as nib 
which have a sizable associated uncertainty. Secondly, the quantity AR is 
a small correction with central values close to zero, primarily due to the 
smallness of the CKM factor /, multiplying Re(<5a). The ranges of 7 and 
R u (as suggested in Chapter 4) imply —0.1 < / < 0.03. The contribution 
to AR specifically from the neutral decays is also small due to an accidental 
cancellation between the 0(a s ) and annihilation effects in 5a. The two main 
sources of uncertainty remaining in the expression i? th are from ( and AR. 

We use our prediction for _R th in combination with the experimental results 
to constrain the possibility of new physics in the radiative decays, or in its 
absence, to determine the CKM factor |Vtd/Ks|- We can use R t h in the 
following ways - either to extract the value of |V i( i/V is | assuming AR is fully 
known, or to extract information about AR using the best-fitted value for 
the CKM ratio from another source. 

We begin by analysing the individual branching ratios for B — > p7 and B — > 
K*j before going on to study the ratio R t h- We summarise the experimental 
results for the radiative decays in Table 5.1 [11, 122, 130, 131]. 



Using the formulae collected above we can obtain a theoretical estimate for 
the combined and individual branching ratios for B — > p'j. With a conser- 
vative estimate of the uncertainty on our values, and using the asymptotic 
distribution amplitude of the rho and omega mesons (A(p) = A(cu) = 1), 



we find BR(B ± -> p± 7 ) = (1.52 ± 0.45) x lO" 6 and BR{B° -> p° 7 ) = 



A c fly -|- X u a 



-W 7 1- 




(5.10) 



5.2.1 B^(p,u)-f 



Decay channel BaBar Belle Average 



B fSxp (B^>K*i) 40.6 ±2.6 43.0 ± 2.5 42.0 ± 1.7 

B exp (B° -> K*°-f) 39.2 ±2.0 ±2.4 40.1 ±2.1 ± 1.7 40.1 ±2.0 

B exp (B+ -> K* + -f) 38.7 ±2.8 ±2.6 42.5 ± 3.1 ± 2.4 40.3 ± 2.6 
B cxp (B -> (p, w ) 7 ) < 1.9 l-34l°inio 

i? exp [(pu;) 7 /ir 7 ] 0.032 ± 0.008±g-ggi 



Table 5.1: The CP-averaged branching fractions for exclusive B — > K*7 and 
-B — > p7, in units of 10~ 6 (with exception of R exp ) and using the 
combined branching fractions defined in the text. 



(0.72 ± 0.22) x 10 6 . We also find, for the combined weighted averages 

BR(B -> (p, u))i) = (1.60 ± 0.48) x 10~ 6 

which is above, but within the error bounds of the current experimental 
measurements. We take the values of the form factors as F p = 0.267 ± 0.021 
and F w = 0.242 ± 0.022 [58]. The error on the branching ratio is estimated 
by taking into account the dominant uncertainties on the input parameters: 
F p , \b, rrib, m c , /# and the CKM inputs. Adding these in quadrature 
leads to an error estimate of O(30%) for the B — > p7 channels. 

The distribution amplitude for the rho meson required for the calculation 
of the hard-spectator and annihilation contributions to the process is that 
for the transversely polarised mesons 0j_(-u,p). The leading twist DA (both 
the transverse and longitudinal) is often expressed as a truncated conformal 
expansion, where there is a prediction for from QCD sum rules [69, 70] of 
a,2 {Pi lGeV) = 0.2±0.1. In this analysis we replace the conformal expression 
of 4>±(p) with our resummed DA model. The suggested region for A was 
examined in [58] , and is based on the indication from all available calculations 
that a 2 > 0. This demands, as for <f>(n), that A > 1. Using the model DA 
with a = 3 suggests A(p) = A(uj) = 1.15 ± 0.10. In light of this we use a 
"physical range" for A of 1 < A < 1.4. 

We test the effect of our non-standard DA by plotting the dependence of the 
b — > g?7 branching ratios on the parameter A, as is shown in Figure 5.2. 

This figure shows that within this range the dependence on A appears rel- 
atively small. There is a downward trend in the p + 7 case and conversely 
an upward trend in the p°7 case. These conspire to cancel in the combined 
branching fraction, leaving B — > (p, 07)7 free of any significant dependence 
on A and of uncertainty from the distribution amplitude. The actual change 
in the value of the branching ratio over the full (physical) range of A is of 
the order of 3% for (p° 7 ) and 1% for the combined ratio. 
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Figure 5.2: Branching fractions for B — > p~f decays from QCD factorisation using 
our model DA </>^~ 3 , in units of 10 -6 and in possible physical region 
of A: B + — ► p + 7 (green), i? — > /5°7 (red) and the combined -B — > 
(p, w)7 (blue). 



5.2.2 5 -». i^* 7 

The theoretical estimate for the B — > K*^ decays can be found in an anal- 
ogous way to B — > p7, using the definition of the weighted average from 
equation (5.6). We obtain an estimate for the A"* 7 branching ratio as 
BR(_B — > K*j) = (5.49 ± 1.64) x 10~ 5 , which is somewhat higher than, 
but within the error bounds of, the current experimental measurement. The 
considerable decrease between this value and that presented in [4] is predom- 
inantly due to the updated value of the form factor F K * = 0.333 ± 0.028 [58] 
from its previous value of F K = 0.38 ± 0.06 [64]. 

We begin by considering the effect of the K* distribution amplitude, by 
testing the dependence of the branching fractions on our non-standard DA. 
Unlike the p, we have the additional complication for the K* meson of an 
anti-symmetric part to the distribution amplitude. The moments for the 
truncated conformal expression are largely unknown, although there have 
been predictions from QCD sum rules [58,99]. This gives us: 

ai(K*, lGeV) = 0.10 ± 0.07 4(iT, lGeV) = 0.13 ± 0.08 

We can examine the dependence on both the symmetric part of the K* DA 
and on the leading anti-symmetric moment af\ The value of the symmetric 
contribution to A for the K* is estimated at A(K*) = 1.12 ±0.10 - we 
consider the wider range of 1 < A (A"*) < 1.4, in order to compare directly 
with the physical range used in the pj case. The dependence of the B — > K*^ 
branching ratios on A with a constant value of af = 0.13 is summarised in 
Figure 5.3. Although this graph suggests that the K* branching ratios have 
a larger dependence that the B — > /ry decays, the percentage change across 
the full range of A is in fact the same for both sets of decays. As before, 



we have a 1% change for the combined K* decay and a 3% change for K*°. 
This effect suggests that the discrepancies between the experimental values 
and the theory predictions can not be attributed to the uncertainty from the 
DA - increasing the contribution of the higher order moments of the DA 
by increasing A serves only to increase the prediction of the averaged K* 
branching ratio. 




Figure 5.3: Branching fractions (in units of 10 -6 ) for B — ► K*j decays in QCD 
Factorisation using our model DA with antisymmetric part t/j^, in the 
possible physical region of A and constant a\ = 0.13; B + — ► K* + j 
(green), B — ► K*°j (red) and the combined B — ► K*j (blue). 

Figure 5.4 shows the dependence of the B —> K* branching ratio on both 
A and ai, and shows clearly that a higher value of the moment a± will lead 
to lower values of the branching ratio for all values of A (if*), although the 
change is not considerable. 
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Figure 5.4: Density plot of BR(i? — ► K*j) (with increasing values indicated by 
lighter shading) for variation in A (horizontal axis) and a\ (vertical 
axis). Contours shown at 5.5, 5.52 . . . 5.58 x 10 -5 . 



5.2.3 R[(p,uj)^/K*j] and the extraction of \V t d/V ts 



This section proceeds to build on the results from the individual branching 
fractions, by calculating the ratio R t h \{p,uj)^ / K*^] as introduced above. As 
discussed, this is important for the extraction of information on the hadronic 
quantity AR and the value of the CKM ratio |V^/V^|. 

We begin by using equation (5.7), and the calculation of AR as in equations 
(5.8) and (5.9) to find the ratio of the (p, u)'~f and K*^ branching ratios. We 
use the asymptotic distribution amplitudes for <fi±(p), (f>±(uj) and 4>±(K*), 
and obtain an estimate of R t h as 



R 



tb 



BR(£?^(p,^) 7 ) 
BR(B -> K*i) 



0.030 ±0.008 



which is in rather good agreement with the experimental measurement. The 
theoretical uncertainty remaining in i? t h is dominated by the form factor 
ratio C, and the value of the quantity AR. We can use the ratio to extract 
a value of |Vtd/Vt s |. Using C^ 1 = 1.25 and our estimate of AR from QCD 
factorisation as AR = 0.006, we obtain: 



\V td /V t 



ts 



0.195 ±0.051 



which is well within the SM range of \V td /V ts \ = 0.197 ±0.013 [132]. We plot 
the dependence of R t h against |Vtd/Vt s | in Figure 5.5. 




Figure 5.5: Plot of Rth [(p, uj)j/K*j] (within quoted ±1<7 errors) against 
I^W^isI with current Belle experimental result and la errors shown 
(horizontal band in magenta) and the SM best-fit for |Vjd/Vj s | (ver- 
tical band). 

The constraints on the unitarity fit are continually improving, so we are 
approaching the stage of having consistent measurements for many of the 
CKM parameters. We can then in principle use the measured (fitted) value 
for |V^/^4s| in order to examine the parametric uncertainties entering into 
-Rth- 
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We expect the dependence on the distribution amplitudes that enter into R t h 
to be small, even though there is dependence on A(p), A(K*) and af*. We 
approximate by taking A (a;) = A(p). Increasing the contribution from A(p) 
acts to increase the ratio, and increasing A(K*) acts to decrease it. Since 
these effects are of the same magnitude, the overall effect is negligible. 

The dependence of -Rth on the form factor ratio is much greater and is illus- 
trated in Figure 5.6. This shows the sensitivity to changes in the value of 
(~ 1 especially any movement towards lower values. 
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Figure 5.6: Plot of [(p, uj) r y/K*'y\ (within quoted ±1<7 errors) against form 
factor ratio £ _1 = (0)/T^(0) (vertical band) with current Belle 
experimental result and la errors shown (horizontal band in ma- 
genta) . 

Our prediction for i? th \(p,uS) r y/K* r y\ is close to the central value determined 
experimentally, showing clearly how well the uncertainties can be reduced by 
calculating the combined ratio. The parametric uncertainty in the individual 
branching ratios is dominated by the dependence on the form factors which 
enter quadratically in the expressions for BR(I? — > V7). The ratio R t h 
reduces the sensitivities to F p and F K * by limiting the impact of the long- 
distance hadronic physics to the SU(3) breaking in the form factor ratio (. 
The other main uncertainty is from the QCDF prediction for AR; this is 
largely due to Ab, as it is this that determines the strength of the weak 
annihilation contributions. This is especially relevant for the charged modes 
£^p± 7 . 

It is worth questioning however, why the individual branching ratios have 
high central values with respect to their experimental determinations, most 
notably, the B — > K*j branching ratio. From the calculation of inclusive 
B — > X s 7, which is in good agreement with the data, it is justifiable to 
conclude that the short distance physics is compatible with the Standard 
Model expectations [126]. This suggests two lines of reasoning: firstly, that 
the uncertainties which are removed when calculating the ratio Rth have a 
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considerable effect on the branching ratio. We illustrate the sensitivity to the 
form factor F K in the K*^ branching ratio in Figure 5.7. This shows how 
the experimental measurement leans towards a lower value for this result and 
how the reduction to the updated value from [58] leads to such a significant 
drop in the theory estimate of BK(B — > K*^). 
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Figure 5.7: Plot of BR(S — > If* 7) in units of 10~5 (and within conservative 
estimated errors) against form factor F K (0) (vertical band) with 
current experimental average (horizontal band in magenta). 

A second line of reasoning is to consider if there are any other factors con- 
tributing to the over-estimation of the branching ratios, which conspire to 
cancel in the combined ratio i? t h- This could include errors from the cal- 
culation of AR, large power corrections which have been neglected, or new 
physics in the underlying b — > (d, s)7 transitions. 

We now go on to consider the possibility of new physics in the B — > Vj 
decays. The errors on the experimental data are reducing to a level where it 
will be possible to impose considerable constraints on new physics models. As 
we discussed in the opening motivation of this chapter, the radiative decays 
(especially those involving the b — > s transition) have long been seen as an 
excellent candidate for searches for new physics. We proceed to perform 
an analysis in a generic minimal supersymmetric model as an example of 
the most popular extension of the Standard Model. We show how effective 
constraints can be imposed on the parameter space for possible new physics 
effects in the radiative decays using the bounds from the experimental results. 
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5.3 Introduction and motivation for super- 
symmetry 



Supersymmetry (SUSY) is the most popular low-energy extension of the 
standard model (see for example [133,134]). Beyond the elegance of the 
symmetry itself, which unifies fundamental fermions and scalars, supersym- 
metric models have many virtues which can alleviate some of the issues within 
the Standard Model. 

The main benefit that motivates many SUSY models is the resolution of 
the "hierarchy problem". In the Standard Model (SM) there are a number 
of corrections which lead to infinite contributions to a particular process - 
these are removed by renormalisation by introducing some cut-off energy A. 
Quadratic divergences appear in the contribution of heavy fermions to the 
Higgs self energy. These can be again be removed by renormalisation, but 
large corrections (~ mj) remain. Hence fine-tuning is required to ensure the 
small Higgs mass required in order to produce a vacuum expectation value 
of the correct size to give the observed W and Z boson masses. If we are 
envisioning that our theory contains the large scales of grand unification, 
namely the Planck scale (mpL ~ 10 19 GeV), then we have to ask how we can 
arrange for the Higgs mass term (and the corresponding electroweak scale) 
to be so much smaller then the underlying mass scale mpL- Supersymmetry 
requires that each scalar particle has a fermionic "superpartner" (and vice 
versa) of the same mass so that all the quadratic divergences of the scalar 
mass terms will automatically vanish. These cancellations occur at every 
order in perturbation theory protecting the light scalar masses. 

Supersymmetry is also a natural bridge to incorporate gravity, as implement- 
ing SUSY as a local gauge symmetry demands local co-ordinate covariance i.e 
linking to general relativity. Supersymmetry is also an integral requirement 
at some energy scale in most string theories [135]. 

If SUSY were to be manifest then all particles and their superpartners would 
be mass degenerate. Since to date, no state which can be identified as a 
superpartner has been observed, then we must have broken supersymmetry 
at our low energy scales. The mechanism for the SUSY breaking is unknown, 
but as we discuss below it has a large impact on phenomenology. We can 
give rough bounds for the masses of the superpartners; all of the masses 
must be over lOOGeV (otherwise they should have been observed at LEP), 
but at least some masses have to be less than 1 TeV, in order to resolve the 
hierarchy problem. 
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5.3.1 Low energy MSSM 



The Minimal Supersymmetric Standard Model (MSSM) is the simplest and 
most popular SUSY extension to the SM. The particle content is a set of 
supermultiplets, constructed from the SM states and their superpartners. 
The superpartners of the left and right handed quarks are denoted as left 
and right handed squarks even though the squarks are scalar and have no 
chiral structure. In summary we have the following: 
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The interactions among the scalar particles are governed by a superpotential 
W, where the dimensionless numbers yi contain the Yukawa couplings of the 
SM fermion fields. 

W = Uy u QH u + Dy d QH d + Ey e LH d + fiH u H d (5.11) 

The gauge bosons of the Standard Model also get a fermionic superpartner 
generically called a gaugino. The gluino is the superpartner of the gluon, 
and hence carries colour charge. The remaining gauginos corresponding to 
superpartners of the W, Z and photon, combined with the higgsinos of the 
same quantum numbers all mix together to give four new physical states 
- two charged and two neutral fermions labelled charginos and neutralinos 
respectively. 

Since SUSY cannot be manifest, we need to find a mechanism to break the 
symmetry. It is in principle possible to achieve spontaneous symmetry break- 
ing of global SUSY, however it is very difficult to build a realistic model. The 
problems with the different spontaneous SUSY breaking scenarios include, 
for example, violation of mass sum rules between quarks and squarks or the 
need for additional field content beyond the MSSM. It is instead expected 
that SUSY is softly broken by dimensionful parameters such as squark mass 
terms. These explicit breaking terms are put in by hand, and can be thought 
of as originating from dynamics beyond the MSSM operating at some very 
high energy scale and "communicated" to the low energy sector via some (un- 
known) mechanism. The soft breaking Lagrangian contains a host of mass 
terms for scalar particles and gauginos, bi- and tri-linear scalar interactions, 
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giving a total of 105 new physical parameters in its most general form. In 
general scenarios there is no reason why these parameters should be flavour- 
blind, but often in particular models many simplifying assumptions are made 
about the structure of the flavour sector. We consider the implications of this 
in the next section. 



5.3.2 Implications for the flavour sector 

Given the most general form of supersymmetry breaking as we have discussed 
above, there are many new parameters and contributions in the flavour sec- 
tor which could impact on the phenomenology of flavour changing and CP 
violating processes. There is potential for many supersymmetric contribu- 
tions to these processes, allowing them to arise at rates much higher than 
experimentally observed. This is however coupled with the very stringent 
constraints on any new physics, for example, from the fitting of the unitarity 
triangle - which combines all available experimental data on flavour and CP 
violation. The success with which this fit determines the CKM parameters 
to be in agreement with the Standard Model expectations shows where new 
large sources of flavour changing neutral currents or CP violation are not 
favoured in a possible new physics model. 

There are a large number of different possible scenarios for SUSY breaking 
[136], including the simplifying assumption of flavour universality at the high 
energy scale. This universality removes all of the terms which could lead to 
"dangerous" FCNC and CP-violating effects which are already excluded by 
experiment. For example, in the minimal supergravity (mSUGRA) model, 
which imposes gravity-mediated SUSY breaking, the canonical scenario has 
universal scalar masses 

However, even in this picture evolving the dynamics from the high energy 
( "GUT" scale) to the electroweak scale induces flavour non-universality due 
to the contribution from the Yukawa couplings [137], although they are ex- 
pected to be small. 

Alternatively, in other SUSY breaking scenarios, it is possible for the dynam- 
ics of the hidden sector to be communicated to the low energy theory in a 
flavour specific way - even before radiative corrections or the RGE evolution 
are included. This will allow flavour changing neutral currents which are in 
general not GIM-suppressed or CP invariant. This clearly will produce a 
very rich phenomenology with many new sources of flavour and CP physics. 



5.3.3 Mass insertion approximation 



To test the impact of a particular SUSY model on the flavour sector we 
would, ideally, diagonalise the squark mass matrices so we could calculate 
diagrams with the squark mass eigenstates. Without knowledge of the full 
mass matrices however, this is a rather difficult problem, as the matrices 
are related to the parameters of the soft SUSY breaking Lagrangian in a 
non-trivial manner. We instead utilise the fact that the flavour changing 
effects originate from the off-diagonal elements of the squark (sfermion) mass 
matrices. If the off-diagonal elements are small relative to the diagonal entries 
we can use the mass insertion approximation [123]. This enables us to simply 
compute ratios of the off-diagonal entries over the diagonal elements of the 
squark mass matrices and compare with the bounds we can derive from 
experimental data. 

In order to allow the most general method of SUSY breaking and flavour 
structure, we can consider a "generic" R-parity conserving MSSM. We pa- 
rameterise the FCNC and left-right squark mixing by taking the squark 
mass matrices as both flavour universal and real. This mass matrix is writ- 
ten with the off-diagonal elements denoted (Aab)^, where AB stands for 
LL, RR, LR, RL and denotes the helicity of the squarks. For example, m 2 LL 
can be written 

/ (m| L ) n (A LL )i2 (A LL ) 13 \ 
™ 2 ll = (All) 2 i (m| L ) 22 (A LL ) 23 

V (A LL ) 3 i (A LL ) 32 (m| L ) 33 / 

We can then constrain the perturbations from universal scalar masses by 
normalising the off-diagonal elements by a common mass, normally taken as 
the average scalar mass m 1 -. The mass insertions are defined as 

Ksh = {A 4^ (5-12) 
m f 

For the b — > d transition the relevant insertions are (<5l B )i3 and for the 
b — > s transitions (5^b) 23 . The universal scalar mass is often taken as an 
approximate mass scale, such as the average squark mass ~ 500 GeV. 

Different realisations of the MSSM will be characterised by different sets of 
soft breaking terms which in turn induce different sets and combinations 
of mass insertions. We classify the effect of these different MSSM models 
on 5fj in terms of the helicities AB = LL, RR, LR, RL. We consider for 
the most part cases where a single helicity mass insertion dominates, but 
also test some cases where two of the mass insertions are sizable (a double 
mass insertion). This allows us to find the leading contributions to various 



processes by considering diagrams with a single (or double) mass insertion 
and express the phenomenological bounds in a model-independent way. 



5.4 Gluino contributions in generic MSSM 



We add supersymmetric contributions to the Standard Model expressions 
using an effective Hamiltonian induced by gluino exchange as first derived 
in [138] using the mass insertion approximation. Gluino exchange is the 
dominant contribution and we do not consider charged Higgs, chargino or 
neutralino exchanges. This is a good approximation in a preliminary analysis 
of this nature as these contributions are proportional to the weak coupling 
rather than a s as for the gluino contributions, and so can be safely neglected. 

Two classes of diagrams contribute to the AB = 1 processes, namely penguin 
and box diagrams, which give both additional contributions to the existing 
Standard Model operators Qi, and induce new operators Qi obtained from 
the Qi under the exchange L <-> R. The gluino contributions which interfere 
with the SM contributions (i.e those from operators Qi) are combined at the 
level of the Wilson coefficients. The contributions from the remaining new 
operators Qi are added separately since they do not interfere with the SM 
amplitude. 

The Wilson coefficients for the SUSY operators act to modify the initial 
conditions of the SM operators at the scale /j, = M w . For the radiative 
processes, we only need to consider the relevant contributions, namely to 

Ci° )cS and Cf eff : 
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C 8 



(^J^ + ^^M^x) 

+ ^)l R ^X~1 Mi{x) ~ 3M2{x] 



(5.13) 



where the Mi{x) are functions of x — m^/rn^ obtained from the calculation 
of the gluino penguins [138], (m, is the average squark mass and rrig is the 
gluino mass). In all of the following analysis we use equal values for the 
average squark and gluino masses, i.e. x — 1. The co-efficients for Q 7 and 
Qs are obtained from these under the exchange L <-> R. 

In considering the SUSY contributions to B — > U7 the Standard Model con- 
tribution is added to the gluino contribution and we determine the possible 
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constraints on the flavour- violating sources in the squark sector. The bounds 
we determine on the 5f 3 and 8$ 3 insertions are only indicative as they are ex- 
tracted ignoring the error of the theory calculation. We also indicate the 
bound on the parameter space when the la theory error is included. 

5.5 New physics in b —> d transitions 

Although the constraints on the flavour sector are quite stringent there is still 
room for new physics to affect the flavour sector. Using the new results for the 
B — > p7 decays we can constrain the possibility of new physics in a generic 
super symmetric model using the mass insertion approximation introduced 
above. Constraints on Sf 3 derived from the B d — B d mass difference AM d 
and the CP asymmetry in B — > J/^K S were determined in [139]. This study 
found that with NLO evolution, for a single mass insertion and trig ~ rrig, 
the limits on 5f 3 are of the order of 10 _1 for the LL or RR insertions, and 
~ 1(T 2 for LR, RL. 

We expect the radiative b — > d decays to show a similar hierarchy, as the LR 
and RL insertions are those most effectively constrained by measurements of 
B — > p7 as they enter into the Wilson coefficients with the factor (m~ g /m b ). 
The LL and LR insertions contribute to the same operator that is responsible 
for the B — > U7 decay in the SM, and so yield a SUSY contribution which 
interferes with the SM. As a consequence, the rates tend to be larger than 
the RR (and RL) cases since these insertions do not add to the leading SM 
amplitude. The consequences of this in light of our high central values for 
the SM amplitudes, are that the resulting parameter spaces for the RR and 
RL insertions will be more constrained than the LL and LR ones. 

With these considerations in mind, we use the gluino contributions to AB = 
1 transitions and the constraints from the newly measured radiative processes 
to constrain the 5f 3 insertions. We can then determine if the measurement 
of the B — > p7 decay can put stringent limits on the possible SUSY models 
that may have impact on the b — > d sector. 

We consider single mass insertions for each helicity combination, and dou- 
ble mass insertions for a left-right symmetric case: (Sf 3 )iL = ($13) rr and 
($i3)rl = (Si 3 )lr as motivated in a number of SUSY models e.g. left-right 
super symmetric models. 

5.5.1 Constraints on 5f 3 

Using the full SM + SUSY contribution to the B — > pj processes, we can 
determine the allowed parameter space in the Sf 3 plane using the available 



experimental results. We use both the B — ► (p, lu) 7 branching fraction 
and the ratio R [(p,uj)'y/K*'y] as these were fitted separately from the Belle 
data set [122], so we can apply the constraints concurrently We begin by 
considering a single mass insertion with an average squark mass of rriq = 
500GeV. For constraints based on the ratio R t h we take all insertions to 
B — > K*j to be zero in the first instance. 

The absolute value of the B — > (p, uj)^ branching ratio is plotted against the 
absolute value of the mass insertion (for each single mass insertion dominat- 
ing in turn) in Figure 5.8. 
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Figure 5.8: Left: Dependence of BR(i3 — ► (p, (jj)'y) with gluino contributions on 
the single mass insertion |(<^ 3 )ll| (solid line) and \(5f 3 )RR\ (dashed 
line) with rriq = m g = 500GeV. Experimental la bounds shown for 
comparison. Right: same for \{5i 3 )lr\ (solid) and \(Sf 3 )jiL\ (dashed). 

The difference between the LL and RR cases is due to the way the gluino 
contribution interferes with the Standard Model. For the case where the 
(^13) LL insertion dominates, the supersymmetric contribution is to the same 
operator responsible for the B — > p7 decay. This interference in general 
causes the decay to proceed at a larger rate than that of the RR case. For 
the RR case where no interference with the Standard Model is present, the 
results will be symmetric around (Sf 3 )nR = 0, and similarly for (5f 3 )nL- 

The allowed values for the LR and RL insertions are expected to be much 
smaller than for the LL or RR insertions as they enter into the Wilson 
coefficients with factor {m g /m h ). Hence (df 3 )LR and (5f 3 )jiL will De the most 
effectively constrained. The high central value of SM amplitude for BR(5 — ► 
(p, uj)j) will lead to a minimal allowed area of parameter space for the RR 
and RL insertions. If the error on the SM amplitude is taken into account 
the constraints the parameter space would clearly be less stringent. 

The full constraints on the parameter space for the four single mass insertions 
are shown in Figure 5.9, where we extract the bounds on the insertions using 
both the central value of the SM amplitude for B —>■ (p, 07)7 and taking into 
account the la error. 
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Figure 5.9: Allowed regions in Re(5f 3 )ij - Im(5f 3 )ij parameter space for ij = 
LL, RR, LR, RL. We use rriq = rrig = 500 GeV, and constraints 
from Belle measurement of BR(£? — * (p, u) 7) and additionally 
R [(p, uf)^/K*^f\. The overlay in magenta indicates allowed regions 
using central value of SM amplitude. 
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The allowed areas in the parameter space indicate regions where the insertion 
can produce a supersymmetric contribution which is within the experimental 
la bounds on BR(£> — ► (p, These figures show that using the data 

imposes quite stringent restrictions on the parameter space, the most impact 
coming from the low error on the R exp [(p, u)^f/K*^\ measurement. 

We can also consider the effect of a double mass insertion, for the cases where 
{813) ll = {813) rr or {Sf 3 ) L R = {8f 3 )RL- This sort of scenario is motivated by 
a number of different SUSY models, such as the left-right supersymmetric 
model or minimal supergravity (mSUGRA) which demands 5ll = Srr = 
and Slr = <5rz- The results for both cases are shown in figure 5.10. 
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Figure 5.10: Allowed regions in Re(5f 3 )ij - Im(df 3 )ij parameter space for double 
mass insertions LL = RR and LR = RL. We use rriq = rrig = 500 
GeV, and constraints from Belle measurement of BR(.B — > (p, uj) 7) 
and R[(p,u)y/K*-y). 

We see a similar pattern as with the single mass insertion, of a much greater 
constraint on the helicity-changing insertions than the helicity conserving. 
The LL = RR case for example, allows ample space for a sizable mass 
insertion. 



5.6 New physics in b —> s transitions 

We now perform a similar analysis with the B — ► K* / y decays. Although 
there are strong constraints on new physics contributions to s — > d and 
b — > d transitions, this is not the case with b —>■ s transitions. The unitarity 
fit for example strongly constrains the s — > d and b — > d transitions, but 
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not the b —>■ s transitions as they do not affect the fit directly (unless they 
interfere with the amplitude for B s mixing and reduce the expected value of 
AM S below its experimental lower bound). In the radiative processes such 
as b —>■ sj, there can be many new contributions with SUSY particles in the 
penguin loop. The main contribution is shown diagrammatically in Figure 
5.11. The cross on the internal squark line represents a flavour changing 
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Figure 5.11: Main supersymmetric contribution to b — > sj. 

mass insertion. The size of this will be determined by the nature of the 
supersymmetric breaking - our ignorance of this means that the insertion 
could in principle be as large as any of the other terms. 



5.6.1 Constraints on 

There have been previous studies on constraining new physics in b — > 57 de- 
cays such as [117], which concentrates on using the bounds from the inclusive 
process B — > X s j. These studies have shown that there is still ample room in 
the parameter space for SUSY contributions to various different observables. 
The constraints that can be determined from the branching ratio and CP 
asymmetry in B — > X s 7 mainly affect the helicity-changing contributions 
i.e. the SUSY insertions (S^lr an d (S^rl- We can therefore use the mea- 
surement of our exclusive decay B — > K*j, and see if additional constraints 
(specifically in the LL and RR sectors) can be determined. 

The graphs in Figure 5.12 show the constraints on the Re(Sf 3 )ij - Im(5f 3 )ij 
plane once the bounds from the B — > K*^y branching ratio have been applied. 

These s show that using the central values for the Standard Model ampli- 
tude for the K*^ decay gives extremely tight constraints on the ReiS^ij 
- Im(52 3 )ij parameter space, leaving virtually no room for any significant 
contribution from a single mass insertion. Also, since there is no available 
parameter space for the RR and RL single insertions, no reduction of the 
constraints is obtained by considering the double mass insertions. Including 
the la theory error on the SM contribution opens up the parameter space a 
little, as we see in Figure 5.12, specifically with the LL and RR insertions. 
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Figure 5.12: Allowed regions in i?e(#j 3 )ij - Imft^ij parameter space for double 
mass insertions LL = RR and LR = RL. We use rriq = rrig = 500 
GeV, and constraints from Belle measurement of BR(£> — > K*~f). 



Overall however, these results do not support the conclusion that there is 
any new physics in the exclusive K*j decays. 
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Chapter 6 

Conclusions and outlook 



Few can foresee whither their road will lead them, till they come 
to its end 

"The Lord of the Rings", J.R.R.Tolkien 

This thesis has examined the non-leptonic decays both within QCD factori- 
sation and beyond it, challenging the assumptions and limitations of the 
method. Our study was three-fold: we analysed the treatment of the light- 
cone distribution amplitudes for the light mesons and introduced a new model 
described by simple physical parameters. We then scrutinised the charmless 
non-leptonic decays in the context of their incalculable non-factorisable con- 
tributions - which we consider the main limitation of the QCD factorisation 
method. Finally, we analysed the new results on the exclusive B — > V'-f 
within QCDF and generic supersymmetric models. 

Our main new result is the development of a new model for the symmetric 
part of the leading-twist light-cone distribution amplitude of light mesons. 
This was done with particular reference to the pion, although the results 
are applicable to all light mesons. The models depend on three parameters, 
two which control the fall-off behaviour of the Gegenbauer moments a n in 
n, and a third A, which parameterises the maximum possible impact of the 
higher-order moments on a factorised physical amplitude. This is given by 
the first inverse moment of the DA. The model is defined at a low energy 
scale (pL = 1.2GeV for the tc) and can be easily evolved up to higher scales 
by solving the evolution equation numerically. In contrast to the conformal 
expansion, which in computation becomes numerically unstable at high or- 
ders, our model DAs are very well behaved numerically and can be easily 
implemented in computer codes. We discussed how the evolution equation 
can be solved at the leading order, but this could in principle also be solved 
to next-to-leading order. 



We have introduced a similar model for the antisymmetric part of the distri- 
bution amplitude, relevant for the K, and the K* mesons; these models are 
normalised to a±, the first Gegenbauer moment. For the ir DA, where experi- 
mental data exists, we have formulated constraints on the model parameters, 
likely to be valid for other meson DAs as well. The constraints on the allowed 
values of A are tighter for the models <p~ with an alternating sign fall-off of 
the Gegenbauer moments than for those with a equal sign behaviour <fi + . 

We feel that these models are better suited for estimating hadronic uncer- 
tainties of processes calculated within factorisation than using the standard 
truncated conformal expansion. This is true in particular for processes that 
involve convolutions with singular kernels ~ as the corresponding inte- 
gral, A, is the basic parameter of the model. 

Our models also should prove particularly useful for describing DAs of mesons 
other than the it which are also symmetric by virtue of G-parity, but for 
which no experimental or other reliable theoretical information is available 
- in particular the p, ui and 0. For these particles we argue that existing 
theoretical indications from local QCD sum rules [41, 99, 100] point to the 
DAs being narrower than the ir, which implies A < 1.2 at the scale 1.2GeV. 
These results also imply that a 2 is positive, hence constraining A > 1. For 
the parameter a, which controls the fall-off the a n with n, we have found that 
the perturbative contributions to QCD sum rules indicate it to be 3 [69,70], 
but that smaller values of a are not excluded unless one can rigorously prove 
that all leading-twist DAs must behave as ~ u(l — u) near the endpoints 
u — 0, 1 at all scales. Any further input, either from experiment or lattice 
would be welcome to help constrain the model, however, even without the 
constraints being taken literally our models provide a convenient way to test 
the impact of non-asymptotic DAs on physical amplitudes without resorting 
to the conformal expansion. 

We illustrated the impact of our models on the calculation of the B — > nn 
branching ratios and CP asymmetries and found that the branching ratios 
are more sensitive to the precise values of the model parameters than the 
CP-asymmetries. The observed discrepancy between experimental data and 
the predictions of QCD factorisation cannot however be explained by non- 
standard DAs, which indicates the presence of non-negligible non-factorisable 
contributions — a conclusion that agrees with the findings of other authors, 
e.g. [50,140]. 

In order to quantitatively justify this premise, we analysed in detail the non- 
factorisable Aq CD /m b corrections to B — > ixix decay amplitudes. We showed 
how the experimental information can be used to extract valuable informa- 
tion on the size and nature of these non-factorisable corrections. We aimed 
to move beyond the model-dependent treatment presented in [1], where non- 
factorisable contributions are identified from convolution integrals which have 
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endpoint divergences. The chirally enhanced power corrections which orig- 
inate from hard-spectator scattering diagrams and annihilation topologies, 
are each parameterised by an arbitrary complex number, the modulus of 
which is estimated from regularising the end-point divergences. We replaced 
these terms with generic non-factorisable corrections via a decomposition of 
the B — > 7T7T decays into isospin amplitudes A 1 / 2 and A 3 / 2 . Three scenarios 
were constructed using differing levels of non-factorisable corrections on top 
of the factorisable terms found via QCDF. 

Our analysis included the effect of a non-standard DA, to simulate greater 
impact from higher-order moments, as parameterised by the physical pa- 
rameter A. We found that the maximum contribution allowed within the 
physical range of A was required to give the best-fit to the data. Variation 
of the B — > 7i transition form factor also required an extremal (minimum) 
allowed value, as found from the \ 2 fitting across the full data set. Our 
combined analysis of all of the most up-to-date experimental data available 
on the B — > ixix system (as of September 2005), showed clearly that sizable 
non-factorisable effects are required to bring all of the B — > txtt results into 
agreement. Even for extremal values of our input parameters, the best-fit 
scenarios do not give a reasonable 2a agreement to the experimental data 
until at least 30% contributions are included. For the Belle data and com- 
bined HFAG results, a minimum level of 40% is needed to produce the same 
level of agreement. 

Our third scenario showed how the long-distance enhanced charming pen- 
guin diagrams can contribute constructively to the decay amplitudes. The 
size of the enhancement must be tightly constrained in order to keep within 
the error bounds of the data (particularly for CP asymmetries) and as tak- 
ing \R\ > 1.5 will make the charm penguin the dominant non-factorisable 
effect. The branching ratios and 7r + 7r~ CP asymmetry can be independently 
brought into excellent agreement, however we discussed how the phases of 
the charm correction are at odds for these two scenarios. The best-fit for the 
scenario with 20% NF correction in addition to the charming penguin gives a 
considerable improvement above the scenario with only the 20% correction. 
However, further NF corrections are still required to bring the theory into 
the experimental 2a error bound. Hence, we can conclude that a sizable NF 
correction is still required even in this scenario. 

It will be interesting to watch how any future experimental results on B — > rnr 
develop, particularly those on the branching ratios for B — > 7r + 7r~ and B — > 
7r°7r°. The current situation makes it difficult to make conclusive statements, 
as measurements from BaBar and Belle currently disagree by approximately 
2a for both the 7r + 7r~ and 7r°7r° channels. A notable rise in the central value 
for the 7r + 7r~ channel (as was the case for BaBar and the figures released 
in summer 2004 and 2005) would act favourably for a lower requirement of 
non-factorisable corrections. The converse however, could leave us looking 
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at approaching 50% non-factorisable corrections in order to obtain a good 
2(7 agreement with the HFAG data. 

The results of our analysis leave us with many more questions than when 
we started. When the experimental situation is resolved and the accuracy 
improves the question remains of the source of the disagreement with the 
B — > 7T7T data; whether this can be explained by significant non-factorisable 
effects within the Standard Model, or if it points to new physics beyond 
it. This also presents the more pertinent question of the validity of the 
factorisation scheme itself, if indeed such high levels of power-suppressed 
corrections are required to satisfy the data. The next step in this puzzle will 
be to analyse B — > irK decays, by fitting a set of generic non factorisable 
contributions to the experimental data [104]. 

Finally, we took a step sideways and analysed the recently released data on 
B — > (p, u;)7 decays. We examined the Standard Model predictions within 
QCD factorisation using the technique first presented in [4] using updated 
non perturbative input, e.g. form factors [58], and testing the effect of our 
non standard DA introduced in Chapter 3. We found predictions for B — > fry 
in agreement with that presented in [4] 

BR(B° -> p° 7 ) = (0.72±0.22)xl0~ 6 BR(B ± -> pS) = (1.52±0.45) x 10~ 6 

and estimated the combined branching ratio as 

BR(B -> (p, ui)i) = (1.60 ± 0.48) x 10~ 6 

which is within the error bounds of the new Belle measurement. We found 
that the variation of the branching ratio with the non standard DA was 
minimal, and that there was less than a 3% effect across the full physical 
range of A. We found the same results for an analogous analysis of B — > K*^ 
where we estimated a central value for the branching ratio of 

BR(B -> K*j) = (5.49 ± 1.64) x 10~ 5 

This differs from the central value presented in [4] predominantly due to the 
drop in the value of the B — > K* form factor. We then discussed how these 
two branching ratios can be combined into a ratio R[(p, / 'K*j] which can 
be used to extract the value of the CKM ratio IV^/V^I- We obtained a value 
of 

R[(p, cc07/iT 7 ] = 0.030 ± 0.008 

in excellent agreement with the measured result. This led to an extracted 
value of | Vtrf/Vts| = 0.195 ± 0.051. This shows how well the ratio can reduce 
the theoretical uncertainty. 

Since previously only an upper bound for the B — > (p, cj) 7 has been available, 
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this was an opportune moment to examine how much room for new physics 
remains in the b — > d'j decays. We examined these decays and B — > K*j in 
the context of a generic minimal supersymmetric model, using the dominant 
gluino contributions as derived using the mass insertion approximation. The 
contributions from charged Higgs, neutralino and chargino contributions are 
sub-leading and can be safely neglected in an indicative analysis of this type. 

We found constraints on the parameter space for the mass insertion 5f 3 
which is relevant for the b — > d transitions. We showed how the allowed re- 
gions are very stringently constrained, primarily due to the constraint from 
R[(p, u)^/K*^\. Taking into account the la theory error opens up more 
of the parameter space, but still leaving only the LL and RR insertions as 
realistic possibilities. We examined the parameter space for the 5^ inser- 
tion using the bounds from the B — > K*^ decay. We used the central value 
of the SM amplitude and found virtually no parameter space available for 
a significant mass insertion. This is to be expected, considering the good 
agreement of theory with the measurement for the inclusive decay B — > X s 7, 
which implies that all of the short- distance physics in the b — > 57 transition 
is under good control. 

As we rapidly approach the new era of particle physics heralded by the LHC, 
precision measurements in the flavour sector will become an essential com- 
ponent of testing the Standard Model expectations, and to uncover potential 
sources of new physics. We need therefore complete control over the long 
distance physics to ensure that the theory errors do not remain significantly 
larger than the experimental errors as the measurements continue to im- 
prove. I believe that this work presents an unbiased examination of methods 
to tackle and understand the uncertainties and limitations of the QCD fac- 
torisation formalism and points to a direction in which further study could 
go. 
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The Road goes ever on and on, 
Down from the door where it began. 
Now far ahead the Road has gone 
And I must follow if I can. 

Pursuing it with eager feet, 
Until it meets some larger way, 
Where many paths and errands meet. 
And whither then? I cannot say. 

J.R.R. Tolkien 
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Appendix A 



Wilson coefficients for 1 



A.l 6x6 operator basis 

We give the explicit formulae for the Wilson coefficients in the naive dimen- 
sional regularization scheme. We consider the 6x6 matrix which is made 
up of the operators Qi...Q 6 . Recalling that we can expand the Wilson 
coefficients in powers of a s , we have 

C iM = cfV) + ^C« W + ... (A.l) 

The evolution of the effective Wilson coefficients is governed by the anoma- 
lous dimension matrix (ADM). The expansion of the ADM including only 
gluonic corrections (no electroweak penguins) is, to order a s , 



2 



iin 



The leading order matrix in this expansion is given (using the operators in 
Section 1.2.3) 



7 



(0) 



{ -2 


6 


6 


-2 



























V 















2 
9 


2 
3 


2 
9 


2 

3 


22 
9 


22 
3 


4 
9 


4 
3 


-¥ 


-2 + f 


2/ 
9 


2/ 

3 








2 


-6 


2/ 
9 


2/ 
3 


2/ 
9 


-16 + f 



\ 



/ 



(A.3) 



The number of quark colours has been set to N c = 3, and the number of 
active flavours is denoted by f. At next-to-leading order we have the following 
matrix governing the mixing of current-current and QCD-penguin operators 
among each other: 
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(i) 





21 
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2/ 
9 


7,2/ 

2 ~ 3 


79 
9 


7 
3 


65 
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7 
3 




\ 




7 
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r 3 


21 2/ 
2 9 


202 
243 


1354 
81 


1192 
243 


904 
81 
















5911 | 71/ 
486 ' 9 


5983 , / 
162 ' 3 


2384 71/ 
243 9 


1808 
81 


/ 
3 














379 , 56/ 
18 ^ 243 


91 , 808/ 
6 81 


130 502/ 
9 243 


_14 + 
3 T 


646/ 
81 














61/ 
9 


11/ 

3 


71 1 61/ 

3 ~ 9 


-99 + 


11/ 

3 




V 










682/ 
243 


106/ 
81 


225 , 1676/ 
2 ^ 243 


1343 , 
6 ~ l ~ 


1348/ 
81 


/ 



(A.4) 

At /xo = Mw the matching conditions can be split into the leading order and 
NLO corrections, Cf\M w ) and C^\M W ) respectively. These are given by 



C?\M W ) 



1 for i = 2 
otherwise 



(A.5) 



and 



C?\M W ) 



for i — 1 



n 

2 

-ii for % = 2 

6 

-\E Q (x t ) fori = 3, 5 
\%{x t ) for % = 4, 6 



(A.6) 



1 1 1 



with Eq(x) being the Inami-Lim function 
Eq(x) = Eq(x)-1 



2, . . x(18- llx-x 2 ) x 2 (lh-lQx + Ax 2 ) , , 
3 + 12(1 -*)» + 6(1-*)' 

(A.7) 



and 



* = w < A - 8 > 



A. 2 10 x 10 operator basis, including electroweak 
penguins 

The expansion of the ADM including gluonic and photonic corrections is 

7 = + 4^ + fe) ^ + 4^ (A ' 9) 

where we include the leading QED contributions from but do not in- 
clude the two-loop QCD-QED anomalous dimension matrix 7^ . The other 
matrices are given below, again in the NDR scheme. 

/ -2 6 000 

6 — 2 — - — - 00 

00^ 22 4qq 

6-f-2 + f ^ f 

000 2 -6 00 

00 f -16 + f00 

000 02 -6 

q q -2(u-d/2) 2( u-d/2) -2{u-d/2) 2(u-d/2) q 
9 3 9 3 

0- -- - --00 

UUg g g 3 UU 

p. p. -2{u-d/2) 2( u-d/2) -2(u-d/2) 2(u-d/2) ^ ^ 
\ UU 9 3 9 3 



































































-2 


6 




6 


-2 


/ 



119 



(1) 
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16 N 
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27 ~ t ~ 
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—16 N 
27 


■ 16(u-d/2 ) 
27 
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81 



















10793 
486 


5623 
162 


6020 
243 


2852 
81 
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18 


6 
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81 


243 


81 
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2735 
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1835 








209 
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V 
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602 


1504 
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41 
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81 
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81 
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iV is the number of colours, /, the number of active flavours and u and d 
denote the number of up- and down- type flavours respectively (u + d = /). 
We have also specified / = 5 for presentation of 7^. 
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The additional initial conditions for the Wilson coefficients of order a e are: 





35 
18 



Ct{M w ) = ^*{ 



2 1 



^ (2B (x t ) + C (x t )) 
f(4C (x t ) + D (x t j) 



| (4C (x t ) + D (x t ) + jr (10B (x t ) - 4C (x t ))) 



with 



B (a; t ) 
C (x t ) 



+ 



x t 



4(x t - 1) 4(x t - l) 2 
x t (x t - 6) x t (2 + 3x t ) 



8(x t - 1) 8(x t - l) 2 



\nx t 
In a;j 



for i = 1,4-6,8,10 
for % = 2 
for « = 3 
for i = 7 

for % = 9 
(A.10) 



(A.11) 



~, , 4, ,-19x3 + 25x1 x?(5x? -2x t -6), 4 

D (a; t ) = — InxA — *- + ty \ — — J -\nx t - - 

° v * ; 9 V 36 (x t = 1 3 18(x t - 1 4 9 
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Appendix B 

Additional formulae from QCD 
factorisation 



B.l Decay amplitudes for B —> tttt 

In terms of the QCD factorisation coefficients defined in Section (2.3.1), the 
B — > 7T7T decay amplitudes are written as follows: 

-V2A(B~ -> tTtt ) = [A^ai + |Af (-a 7 + a 9 + a p 10 + r£a?)] A n7T 
-A{B^^-) = [AiV + AWCoS + ^ + AWr^og + al)]^ 

- ttV) = [Afa 2 - Af « - |< ) - A^« - X) + 1^ - «r)] 4nr 

(B.l) 

with aj = aj(7T7r), A^ = V p bV* d for p = u,c and 

A„r = <^(m| - (B.2) 
The charge conjugate decays are obtained by replacing \ p — > A*. 



B.2 Annihilation contributions to B — > V7 

The contributions to the annihilation amplitudes contain two different factors 
depending if the photon emission if from the light quark in the 5-meson, or 
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from one of the constituent quarks of the vector meson. The former give the 
factor b v , and the latter d v , which are expressed as [124]: 

bV = 2tt 2 f B m v f v dV = 4tt 2 f B f$ C 1 <j) v {v) ^ 
F v m B m b \ B ' F v m B m h J v 

The integral over the distribution amplitude in the factor d v will give a 
dependence on A (and ai where relevant) for all mesons. We include these 
components in the decay amplitudes via 



a" — > a" + a 
a c 7 — >■ a c 7 + a 



u 

ann 



where 



<Cm(A) = Qd [-o 2 b p + a A b" + 2a 6 d"(v)] 

<nn(A) = Qd M P + 2a 6 d»(v)} 

<n>~7) = Qu [aibP + a A b p + (Q s /Q u + 1) a 6 d»] 
<nn(P"7) = Q u Kfc p + (g s /g u + l)a 6 d»] 

Cn(^*°7) = Q^^V a 6 (d*>) + 
< nn (7r 7 ) = g d [a 4 6^+a 6 (rf x *(t;) + ^*(t;))] 

< nn (K*- 7 ) = Q u [tuft** + a 4 b K * + Q s /Q u a 6 d K *(v) + a 6 d K \v)] 
a c ann {K^!) = Q u [a 4 b K * + Q s /Q u a 6 d K * (v) + a 6 d K * (v)] 

The a,- are combinations of Wilson coefficients: 



0-1,2 


= c ^ + k 






a 6 
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Appendix C 



Analytic evolution of light-cone 
distribution amplitudes 



The pion distribution amplitude is by nature a non-perturbative quantity, 
but has evolution governed by perturbative QCD. The evolution is described 
by the ER-BL equation given in (3.22); we can recast this as 

<p n (u,fx 2 )= dvU(u,v,fj, 2 ,fil)<l> v (v,fj,l), (C.l) 
Jo 

where the operator U(u, v, /z, /Xq) represents the solution to an evolution equa- 
tion equivalent to (3.22) and describes evolution from the scale /Zq to some 
scale fi 2 . This presents a general solution to NLO [87] of 



\u,v,fi 2 ,fi 2 ) = ^2E n {n 2 ,n 2 ) 



n=0 



( 2\ 00 

CT(2u-l) + ^ £ d kn (»\»t)C 3 k /2 (2u 



k=n+2 



+ 0(ai 



u(l — u) 



Cl'\2v - 1). 



(C 



At leading order, the mixing of the moments of the distribution amplitude are 
triangular, so that the LO kernal is diagonal with respect to the Gegenbauer 
polynomials. This implies that only the diagonal terms of the (triangular) 
anomalous dimension matrix 7^ = will appear and are completely en- 
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coded in the function E%° (/J? , nl) 
£^V,/^) =exp 

with 



as(p2) ia . a 7n(a a ) 









.«-(A*o). 



7i 0) /2&0 



7i 0) = 



-3 - 



o "+ 1 

- +4V 

(n + l)(ra + 2) ^« 



(C.3) 



(C.4) 



The anomalous dimensions are identical to those known from deep inelastic 
scattering [87,141]. 

At next-to-leading order, the conformal operators mix under renormalisation, 
so that both the diagonal and non-diagonal terms will contribute to the 
evolution. The equivalent function to (C.3), E^ LO , is written 





lk 0) /2b0 - 







1 + 



1 + A *M) 



w(n) 



(C.5) 



with 



win) 



i , (0), 
Tn Op - 7n 0i 

2b bi 



(C.6) 



The off-diagonal contributions to U(u,v, fi 2 , (/%) denoted d nk in the full solu- 
tion (C.2) given by 



(C.7) 



where 

c fcn = (2n + 3) 



- 7 i 0) - 2&o + 8C F A kn 2C F [A kn - ^(k + 2) + 
2(k-n){k + n + 3) (n + l)(ra + 2) J 



_ (0) 

7fc 7n 

(0) (0) n , 

7fc - 7n - 26 



«s(/-t 2 ) 



-l+(7i 0) -7i 0) )/2 6 o" 



fc + n + 4 



k — n 



+ 2ip(k - n) - 4>(k + 2) - -0(1) 



(C.S 
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and with the function ip(z) defined as tp(z) = dhl J} z ^ ■ The NLO anomalous 
dimensions are given [141, 142] as 



7 i 1} = -(Cf ~ \C F C A ) llGS^n 



2n + 1 
n 2 (n + 1) 



+ 16 



2S 1 (n) 



n(n + 1) 



+ QAS(n) + 2AS 2 (n) - 3 - SS'^n) - 8 \ f - .„ - 16(-l) n 



[S 2 (n) - 5£(±n)] 
2n 2 + 2n + 1 1 



.43 _ 4 
6 



536 I o 

9 n 2 (n+l) 2 



2 ' n 3 (n + l) 3 v ' n 3 (n+l) 3 

2n + 1 



- 165*1 (n)S 2 (n) + S 2 (n) 



_52 _|_ u 

3 n(ra + l) 



151n 4 + 263n 3 + 97n 2 + 3n + 9\ C F n f 



9n 3 (n+ l) 3 



{-m Sl ( n ) + fS 2 (n) + l 



+ 16 



lln 2 + 5n - 3 



9n 2 (n+ l) 2 
where the n-dependent functions are 



Si(n) 



S[{\n) 



n 1 

E 1 



= 



2 

1H 



(CIO) 



and Cf = §, CU = 3. 

The distribution amplitude to NLO will then be of the form 



<^K/i 2 ) = Qu(l-u)J2 \aM)E^ LO (f, 2 ^l)C 3 J 2 (2u-l) 

n=0 ^ 

^ E an(^)^ LO (^,^)4n(/i 2 ,^)Cf(2 M -l) 

(C.ll) 



+ 



fc=n+2 



(C.9) 



The Gegenbauer moments a n contain all of the important non-perturbative 
information. They can be independently evolved from their input values at 
some scale /i 2 , to any scale jj 2 using the (NLO) evolution functions discussed 
above. Generalising to include the dependence on all higher-order Gegeg- 
nbauer moments we can write the following as an alternative to equation 



1 1 Q 



(C.ll), using the shorthand a n ($) = a°: 

oo 

Mu, fi 2 ) = 6u(l -u)J2 a n (fi 2 )C 3 J 2 (2u - 1) (C.12) 



n=0 

where 

a n (v 2 ) ee a° n EZ">(v?,fi) + ^J2 d ^^l) E ^ 2 ^iyk (C13) 

fc=0 

for example: 

a 2 (/i 2 ) = a°-E 2 + 7—^20 

47T 

a 4 (/i 2 ) = a°E 4 + -r-(dw + d 42 E 2 al) 

47T 

a 6 (/i 2 ) = ag-^e + t^(^60 + d§ 2 E 2 a\ + d^E^ 

47T 

When the evolution from //jj — >■ /i 2 crosses any heavy flavour thresholds these 
expressions will be modified. Taking /i 2 , at a low scale ~ lGeV 2 where the 
number of active flavours rif is three, the above expressions apply for all 
/i 2 < m 2 . Above this value there are two distinct regions: (i) with rif — 4 
when m 2 < fi 2 < m 2 and (ii) with rif — 5 for /i 2 > m 2 . These can be 
summarised in a similar expression to (C.13) 



k=0 



(C.14) 

where the evolution function E^ 1 ^ takes on a different form in each region 
above the charm and beauty thresholds. 

Region (i) 

E^\^,»l)\ m = E^^ 2 ,m 2 c )E^(mlti) 



(C.15) 
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gion (ii) 



+4 3 2k 2 ,^)^ 5) (^ 2 ,^)4 4) K,^) 

(C.16) 
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Appendix D 

Summary of input parameters 



Masses 


m c 






M w 


1.3GeV 


4.8GeV 


174.3GeV 


80.4GeV 


m n 


m p 




m K * 


140MeV 


770MeV 


782MeV 


894MeV 



CKM Parameters and couplings 


V 


P 


v cb 




0.347 


0.196 


0.0415 


-0.2258 


7 




a 


a s (M z ) 


60.3° ±6.8° 


23.4° 


1/137 


0.1187 



Decay constants [58] 


fs 


u 




fk* 


200 ± 30MeV 


133MeV 


217MeV 


170MeV 


f P 


f 1 - 
J P 


feu 


f 1 - 


205MeV 


160MeV 


195MeV 


145MeV 




Form factors [46, 58] 








F^ 


0.258 ±0.031 


0.333 ±0.028 


0.267 ±0.021 


0.242 ± 0.022 



B-meson parameters 


m B 


\b 


t b + 




5.28GeV 


350 ± 150MeV 


1.65ps 


1.53ps 
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